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ABSTRACT 


The Interaction of an underexpanded hydrogen jet coaxially injected 
into supersonic flow is Investigated experimentally. Experimental re- 
sults are d scussed and analyzed* Comparisons are made between the 
experimental results and theoretical predictions computed using an ana- 
lytical technique. Changes to improve the theory are suggested. 
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SUMMARY 


An experimental data base for the injectlor, mixing, and combustion 
of an underexpanded hydrogen jet in a supersonic test stream has been 
obtained. Experimental pitot pressure data have been compared with theo- 
retical predictions. 

The experimental tests were conducted with both air and nitrogen as 
test media which led to reacting and nonreactin^ flov's, respectively. 
Tests were conducted in a free- jet and in a di ted mode. Theoretical 
values were computed using two different viscosity models and a wii.e 
range of Prandtl number (0.7 to 1.4) with a Lewis number of 1. 

The comparison of the experimental and theoretical data Indicattis 
that the theory is inadequate for predicting the flow field resulting 
from the injection of an underexpanded (hydrogen) jet into supersonic 
flow. Suggestions are made for improving the theory. 
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CHAPTER 1 


INTRODUCTION 

The hydrogen fueled supersonic combustion ramjet (scramjet) engine 
is envisioned as the prime candidate to fill the propulsion requirements 
for future hypersonic aircraft. However, feasible scram jet engines face 
problems in several technological areas. (Status evaluations of the 
scramjet concept may be found in references 1, 2, 3, and 4.) Three such 
areas are of concern in this work* These are the injection, mixing, and 
combustion of hydrogen* Note that the last two are directly related to 
the first by the following sequence: injection controls mixing and mix- 

ing controls combustion. As a result, fuel injection holds an important 
position in the total scramjet problem. Thus, it is not surprising that 
numerous fuel 5' iection schemes have been investigated in both cold and 
hot supersonic flows. Simplicity in flow field modeling has made paral- 
lel coaxial injection the scheme most widely investigated (references 5, 
6, and 7 present investigations of this type). 

These previous investigations of coaxial Injection were limited to 
cases where injector exit pressure matched the test stream static pres- 
sure. These matched pressure cases were selected primarily because the 
theory available was designed to handle them* 

On the other hand, recent theory (see references 8 and 9) is design- 
ed to handle tl:a more complex underexpanded (jet pressure greater than 
the test stream static pressure) injection. The significance of such a 
theory becomes apparent when one notes tnat any practical scramjet engine 
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is likely to use hydrogen injection by an underexpanded jet. 

In fact, all scrarojet engines must be capable of operating with 
underexpanded injection, although this may not be the primary type of 
injection. However, a search of the literature indicates hat there 
was very little information on an underexpanded hydrogen jet coaxially 
injected into supersonic flow. Particular information, such as data on 
the unierexpansion (exit) shock wave’s affect on the hydrogen mixing and 
combustion, is completely lacking. The present investigation was there- 
fore undertaken to experimentally determine some of the fundamental 
characteristics of the mixing and combustion of an underexpanded hydrogen 
jet in supersonic flow. In addition, the theory of reference 8 was test- 
ed by comparing experimental data with theoretical data computed using 
the computer program (reference 9) based on tb.- theory of reference 8. 



CHAPTER II 


APPARATUS AND INSTRUMENTATION 
Facility and Test Conditions 

The experimental portion of this work was conducted in the Langley 

II* Inch Ceramic*Heated Tunnel* This facility, described in reference 10, 

has a bed of zirconia pebbles which is heated by the combustion products 

from a propane burner* The products from the burner are passed through 

the bed until the desired stagnation temperature is reached* The hot 

test gas is obtained by passing the test medium (air or nitrogen) through 

the heated pebbles* In this manner, test gas total teirperatures up to 

2530 K (maximum usage temperature of the zirconia pebbles) can be fur* 

2 

nished with a maximum stagnation pressure of 4 MN/m * 

For the purpose of the present tests, the facility was fitted with 
the Mach 2 test stream nozzle which is a scaled version of one gWen in 
reference 11* This axisymmetric nozzle was constructed of stainless* 
steel and cooled by about 6 kg/sec of water* The facility was operated 
in two modes, a free*jet mode and a du'^ted mode* In the ducted mode, 
the ducting around the supersonic flow formed a circular combustor* A 
schematic of the facility (in the ducted mode configuration) is given in 
figure 1* The free*Jet mode configuration is obtained by removing the 
constant area duct whic* extends from plane A*A to plane B*B of figure 1. 
In each configuration, the exit plane of the Mach 2 hydrogen injector 
nozzle was 0*3175 cm downstream of the exit plane of the test stream 
nozzle* Tests were conducted with both air and nitrogen as test media. 
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Figure 1.- Schematic of the coaxial supersonic combustion apparatus. 




6 


and for all tests the total temperature of the test stream was 2167 K, 

with a nozzle exit (static) temperature of 1338 K* This temperature was 

high enough to give ignition without a pilot flame or ignitor* The 

stagnation pressure ranged from 0.759 to 0.8S8 MN/m , which gave rise to 

test gas flow rates of 1.23 to 1.39 kg/ sec and nozzle exit (static) 

2 

pressures of 0.099 to 0.112 MN/m . 

A summary of the test conditions is presented in Table 1. 

Hydrogen Injector 

The hydrogen injector, which was mounted coaxial with the main noz- 
zle, is a 0.953 cm (3/8 in) stainless-steel tube with a 5^ conical nozzle 
at the exit. This nozzle, with a 0.635 cm exit diameter and 0.488 cm 
throat diameter, gives a nominal exit Mach number of 2. The injector 
exit lip thickness is 0.159 cm. 

This injector lip of finite thickness introduces the problem of 
wake effects in the base region of the injector. However, it is considered 
to be a good compromise between the ideal and technically practical noz- 
zle. Ideally, for ease of analysis, the injector should have an infi- 
nately thin lip, and parallel flow at its exit. Unfortunately, the 
contoured nozzle ueeded to fulfill these ideal conditions cannot be 
built and a compromise must be sought. If the requirement of parallel 
flow is dropped, the infinitely thin lip can be approached by at least 
two designs* One is the boattail conical type nozzle given in figure 2. 
This design produces two undesirable results. First, the boattail causes 
the test flow to expand to a lower pressure, and second, the expansion 
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Table 1 


TEST CONDITIONS 

X 

^'0 

MN/m2 

^00 

MN/m^ 

!i 

P« 

Re X 1 (T * 

m 

kg/sec 

s 

kg/ sec 

Test 

Type 


0.789 

0.103 

2.03 

2.446 

1.302 

0.015 

A-FJ 

19 

0.794 

0.103 

2.031 

2.460 

1.307 

0.015 

A-FJ 

30 

0.802 

0.104 

2.189 

2.485 

1.320 

0.016 

A-FJ 

40 

0.80 

0.104 

2.138 

2.478 

1.316 

0.016 

A-FJ 

56 

0.817 

0.106 

2.022 

2.531 

1.345 

. 

0.015 

A-FJ 

80 

0.791 

0.103 

2.136 

2.451 

1.302 

0.016 



0.792 

0.103 

2.292 

2.452 

1.281 

0.017 

N-FJ 

19 

0.794 

0.103 

2.269 

2.460 

1.285 

0.017 

N-FJ 

30 

0.806 

0.105 

2.356 

2.498 

1.305 

0.018 

N-FJ 

40 

0.782 

0.102 

2.123 

2.425 

1.267 

0.016 

N-FJ 

56 

0.842 

0.11 

2.054 

2.609 

1.363 

0 016 

N-FJ 

80 

0.779 

0.101 

2.205 

2.413 

1.260 

0.0 x 6 

N-FJ 

30 

0.817 

0.106 

2.101 

2.530 

1.343 

0.016 

A-D 

40 

0.812 

0.106 

1.952 

2.517 

1.336 

0.015 

A-D 

96 

0.811 

0.106 

1.966 

2.511 

1.334 

0.015 

A-D 

144 

0.819 



2.538 



A-D 

30 

0.759 

0.099 

2 285 

2.352 

1.229 

0.016 

N-D 

40 

0.833 


2.035 

2.582 

1.349 

0.016 

N-D 

96 

0.781 

0.102 

’..996 

2.422 

1.265 

0.015 

N-D 

144 

0.858 

0.112 

2.05 

2.657 

1.388 

0.016 

N-D 


A - Air FJ - Free Jet N - Nitrogen D - Ducted — » Duct Length 

^ X 
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Figure 2.- Half cross section of a nozzle having a very thin lip end boattall* 
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turns the flow (near the Injector) so that it is no longer parallel with 
the rest of the test stream. The other design is the one of figure 3, 
where the boattail has been eliminated. Unfortunately, this design 
suffers from the increased chance of separation of boundary layer on the 
injector. Such separation of boundary layer would be caused by inter- 
action with the exit shock and the jet flow* If a nozzle of this design 
is cut off (see figure 3), the resulting nozzle has a finite lip thick- 
ness with a base region. Although the wake effects of this region can- 
not be computed close to the base, the probability of boundary layer 
separation is reduced. This result is obtained from the fact that the 
boundary layer can bleed into the wake and the compression effects of 
the divergent flow are eased. It was felt that the exit thickness 
(0.159 cm) of the injector chosen was sufficient to prevent separation 
but small enough to get a far field (several solution for the wake 

region. It is also pointed out that the experimental data of reference 
12 indicates that the jet spreads better when injected from a blunt body 
of this type. Increased spreading (mixing) suggests better burning. 

This design was therefore adopted for the present investigations. 

The cooling needed to protect the injector during each test is pro- 
vided by the injectant (hydrogen). In the present tests, the hydrogen 
supplied at ambient temperature was heated to a total temperature of 


approximately 470 K as it cooled the injector before injection into the 


stream. With this 'total temperature, and stagnation pressures ranging 

2 

from 1.59 to 1.94 MN/m , the injector supplied hydrogen mass flow rates 
of 0.015 to 0.018 kg/sec. The resulting equivalence ratio, based on 
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Figure 3.- Half cross section of a nozzle with no boattall but with a very thin lij 
(Cutting the tip off this nozzle as indicated by the vertical line pro- 
duces an Injector of the type used in the present work.) 
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total flow In the test stream nozzle, varied from 0.381 to 0.467 and 

2 

the exit (static) pressures ranged from 0.203 to 0.248 MN/m . The in- 
jector exit pressure was therefore about 2 times the test stream static 
pressure for each test, and the injected hydrogen was thus underexpanded. 

circular Combustor 

In the ducted mode, constant area ducts of four different lengths 
(9.53, 12.70, 30*48, and 45.72 cm) were individually attached to the 
facility nozzle to form circular combustors. These combustors, con- 
structed of stainless-steel, are uncooled (heat sink) and have numerous 
pressure orifices for measuring static pressure. The orifices are ar- 
ranged in three rows (designated P, Q, and R in figure 4) that run 
axially along the duct with each row spaced 120^ apart. A schematic of 
the 12.7 cm combustor, accompanied by a table summarizing the orifice 
locations for all four ducts, is given in figure 4. 

Pitot Probes 

The pitot probes used in the present tests were of two different 
' jigns. One design is a modified version of a probe developed by the 
Applied Physics Laboratory of Johns Hopkins University and reported in 
reference 13. It has an outside diameter of 0.633 cm and a tip half- 
angle *>f 30^ (see figure 5 for details of probe tip). The other design 
is a slightly modified version of a probe described in reference 14. It 
has an outside tip of 0.914 cm and a tip half-angle of 20^ (see figure 
6 for details of probe). 



1.27 cm 
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Figure 4.- Half section of 12.7 cm length duct (» ken in the plane bisecting 

the pressure taps of Row P) with axial locations for all four ducts. 
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Figure 5.- Cross section of the modified Johns Hopkins' probe. 
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Figure 6.- Half*-section of the modified Eggers* probe. 



15 


Probes of both designs were water-cooled by a no return method. In 
this method, water is supplied through a single passage in the main body 
of the probe, sprayed against the rear of the probe tip, and then Inject- 
ed into the test stream at a location behind the pressure sensing region. 
Once in the test stream, the water is swept downstream over the probe 
body furnishing further cooling. 

Pitot-pressure profiles were obtained with a single moving probe 
which was driven perpendicularly across the flow field at a rate of 
approximately 0.3 cm/sec by a dc motor. Comparisons of pitot pressures 
taken at the same points with the probe moving and stationary indicated 
that response of the pressure transducer was» sufficient to give accurate 
measurements while moving. In addition, probes of either design gave the 
same results for identical test conditions. 

Photographs and Shadowgraphs 

Data obtained in the form of photographic records were of two types: 
black and white movies, and shadowgraphs. The movies were taken at a 
frame rate that varied from 20 to 64 frames/sec. They were used to check 
the pitot probe alinement and vibration. Tlie shadowgraphs were taken at 
a constant frame rate of 24 frames/sec. They were used to define the 
flow quality and are q^j^ite valuable for analyzing the flow field. 

Photographic records of bo^h types were obtained on 16 mm black and 
white mrvie film with an AS^t number of 400 (Lin number of 27). The total 
photographic records will not be included in this work. However, an exam- 
ple of the shadowgraphs are given in figure 7. 



It 



(a) Test stream and jet* 



(b) low-up of wave structure* 


Figure 7*-- Shadowgraph of the test stream and jet in t!*e free- jet 
mode with comboatlon. 
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Example of Pitot Test Data 

Although the majority of the test data Is to be presented In Chapter 
IV, the exit pitot surveys are introduced here to provide a feel for the 
experimental data. In figures 8 and 9, radial exit pitot pressure pro- 
files for the free-jet reacting and nonreacting cases are given respec- 
tively. Both profiles have the same general shape, however their peak 
(centerline) values are not equal. The nonreacting peak value is less 
than the reacting, since it is taken at an axial location slightly down- 
stream of the axial location of the reacting case. The solid line Df 
both figures is a straight line connection of adjacent .ata points in- 
tended as a guide to the data trend. 

Since both cases have the same shape, only one discussion will be 
offered. This discussion uses the letters common to both of these fig- 
ures, and the flow schematic of figure 10. The pitot pressure varies 
racially in the following manner. The pressure decrease in going from 
points a to b is due partly to the radial travel across the conical 
jet flow field, and partly to an expansion fan from the injector lip. 

Both processes result in higher Mach numbers, and thus lower pitot pres- 
sures. Ttu small peak at c is the resu.lt of the shock wave which 
terminates the expansion fan. The decrease in pressure from c to d 
is due to the shock wave indicated at c and the fact that d is in 
the base region of the Inject^'r. The shock ind cated at c is a curved 
shock which extends from the injector lip at the exit to the centerline 
at a slightly downstream location. Thus, much of the region c to d 
is behind the curved shock, whose strength varies from a minimum near the 
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Figure 8.- Pitot profile at exit of Injector (x/r 
with air test medium. 
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Figure 9.- Pitot profile at exit of injector (x/r 
with nitrogen test medium. 
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Test stream mixing boundary 


Jet /Transmitted portion 
mixing / of intersecting shock 
boundary , 



.Expansion fan with 
intersecting shock wave 


Survey location 


Figure 10.- A schematic of the free-jet flow field with 
various prominent features at the survey 
location labeled. 
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injector lip to a maximum at the centerline. The strength variation in 
this region produces a radial pitot pressure profile which varies in the 
same direction (minimum to maximum) » whereas the radial Mach number pro- 
file varies in the opposite direction. In idealized flow, d would be 
the location of a slip line separating the test stream and jet flow. In 
the present work, the radial region near d is probably a mixing bound- 
ary. Point e is the underexpansion or exit shock wave which extends 
from the injector's outer lip to the test stream boundary, where it is 
reflected as an expansion fan. Therefore, the region d to e is 
similar but in opposite sense to the region c to d. The region from 
e to f is the test stream without any interaction* The dip from f 
to g is an indication of the free-jet test stream interacting with the 
ambient air. 

It may be surmised from the above discussion that the flow field 
resulting from the underexpanded injection of hydrogen into supersonic 
flow is quite complex. As a consequence, the theoretical treatment by 
necessity must be rather sophisticated. The theory used for comparison 
in this work is that of reference 8, and is outlined in the next chapter. 



CHAPTER III 


THEORY 

General Governing Equations 

The basic governing equations are the well knowi. ••viscous-lnviscid” 
equations used in higher order boundary layer and viscous flow field 
analysis with tne finite rate chemistry terms included. These “viscous- 
inviscid“ equations are supplemented by the Rankine-Hugoniot and Prandtl- 
Meyer relations to facilitate the computation of shock and expansion 
conditions respectively. The basic equations are given in Appendix A 
along with a limited discussion of how they are applied. The reader 
interested in a more thorough delineation of the equations and the nume. 
ical application may consult references 15, 16, and 17. 

Viscosity Models 

The program as published in reference 9 had a turbulent eddy vis- 
.sity model referred to as the “Ferri-Kleinsteii" model. This model, 
which v^as developed in references 18 and 19, has viscosity variation in 
the axial direction only. However, it was felt that Eggers* viscosity 
model (see reference 20), which varies both axially and radially, may be 
Fire accurate. Thus, it was decided that the program would be run with 
both models individually incorporated. 

Ferri-Kleinstein Model 

In this model, the turbulent eddy viscosity undergoes an axial 
variation from the jet exit to the end of the potential core. The length 
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of the potential core is defined as: the distance x from the Jet exit 
to the downstream location where the mass fraction of hydrogen on the 
centerline becomes less than 0.99. The viscosity is then assumed to be 
constant for all locations downstream of the potential core length x* 
The viscosity is computed, for stream locations {x/r^) less than y 
with the nondimens ional equation, 

g - ((P<1)^, - (P<l)^„) <«/rj + K 3 ( 1 ) 

where, « 7.5 x 10“^ and = 100. 

For stream locations equal to or greater than X 

„-KiRe(<p,)^^- (P,)_j^)x+K3 (2) 

and since y is constant downstream of the length x« ation 2 is 
executed once* The resulting value of y is stored for a 1 future down 
stream calculations. 


Eggers' Model 

There are two viscosity models generally referred to as Eggers* 
model, thus one must be careful to specify the model intended. The two 
models, which are similar in mathematical structure, are called Z-differ 
ence and kinematic Z-difference models by Eggers (reference 6). In the 
Z-difference model (see reference 20) the absolute viscosity varies axi^ 
ally only and is computed using the nondimens ional equation, 

y - KZ (pq)^^ (3) 

In the kinematic Z-difference model, the kinematic viscosity varies 
axially and is computed using the nondimens ional equation. 
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e = KZ (q)^j^ (4) 

The absolute viscosity is obtained by multiplying the kinematic 
viscosity (of equation 4) by the local density which varies radially. 
Thus, the absolute viscosity varies both axially and radially, and is 
computed with the equation. 


U * plocal KZ (q)^^ 

In all three equations (3-5), the empirical constant K has a value 
of 0.01. The quantity Z is defined as the radial distance between the 
points where the local velocities are and U2 as given by the 

equations, 

U = + 0.95 (Up. - U^) (6) 

1 a a 

and. 


0.5 (U^,_ 


where equal the stream velocity external to the jet. 

It is the model computed by use of eqiiatior. 5 that is referred to 
as the Eggers' model in this work. 



CHAPTER IV 


RESULTS AND DISCUSSION 

The experimental data and theoretical predictions of the present 
study are presented in dimensionless form. All pressures are nondlmen- 
sior.alized by dividing by the test stream stagnation pressure (P ). 
Similarly, dimensionless coordinates and lengths are obtained by division 
by the hydrogen jet radius (r^) at the exit of the injector. It is also 
noted that all theoretical calculat'<''’\s were performed with a Lewis 
number of 1. 


Free- jet Data 

Radial pitot pressure surveys were taken at several axial stations 
for the free- jet mode and at the end of the ducts when operating in the 
ducted mode. The pitot pressure data (surveys) for each mode of operation 
can be subdivided into reacting and nonreacting cases. In the reacting 
cases, the test stream is air, and in the nonreacting cases the test 
stream is nitrogen. 

The data for the free- jet reacting cases given in figure 11 are 
typical and will be discussed. The pitot surveys were made at axial 
locations of 1, 19, 30, 40, 56, and 80. The data for the axial 

location equal one were previously presented in figure 8, and will 

not be covered here. The prominent features, such as high jet centerline 
pressure bounded by jet mixing boundaries, present at the ° ^ 

location extend downstream. In fact, the high centerline pressure is 
present for the 30, and 40 locations. However, the mixing 
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Figure 11.- Experimental and theoretical pitot profiles at various 
axial locations for the reacting free-jet mode (Ferri- 
Kleinstein viscosity model)* 
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Figure 11.- Concluded 



31 


region has engulfed the centerline at the x/r^ * 56 location and the 
centerline is not discernible. The regions of no interaction, previously 
discussed in the section on the pitot sample, have become tenuous at the 
x/r^ = 56 location. This demise of these regions is attributed to the 
fact that free- jet tr-^t mixing boundary spreads inward to meet the jet- 
test-stream mixing region which spreads outward. 

Other details of the data ere given in the following discussion, in 
which theoretical predictions are compared with the data. 

The theoretical calculations at the free-jet test stream boundary 
were not expected to agree with the experimental data, since the program 
does not have the necessary theory fot handling the test stream mixing 
boundary. The program tak<.s a constant pressure boundary approach which 
is sufficient for mathematical consistence, but improper for actual 
boundary conditions. This approach does not affect the accuracy of the 
calculations performed for the reg'*on inside the test stream mixing 
boundary since this region is supersonic. Thus, the boundary disturbances 
cannot be transmitted to the internal region of interest, and the calcu- 
lations should be in agreement with the experimental data. Unfortunately, 
an actual comparison of the theoretical calculations and the experimental 
data does not show such agreement. In figure 11, for example, there is a 
comparison of the experimental pitot pressure data to theoretical calcu- 
lations perfonsea with the Ferri-Kleinsteln viscosity model. The test 
stream is air, and ‘the theoretical data are for Prandtl numbers of 0.7 and 
1.4. As expected, there is no agreement in the region of the test stream 
mixing boundary. For axial locations x/r^ ■ 19 and 30 where there is a 
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region of test stream not affected by the mixing boundary or jet inter- 
action, the agreement is excellent. (At x/r^ = 19 these regions extend 
from y/r^ ■ -7.5 to -4 and from y/r^ = 4 to 9.) This agreement indi- 
cates that the constant pressure boundary approach does not affect the 
accurary of the program for the region internal to the test stream mixing 
boundary. However, the only other semblance of agreement is at the 
centerline region y/r^ > i 1 and that is not complete. For example, 
the centerline differences between the experimental and theoretical 
values Pr * 0.7 are given in Table 2. 

Table 2 


x/rj 

% bxfference @ g 

19 

21 

30 

0 

40 

13 

56 

56 

80 

^ 


This erratic agreement on the centerline suggests that the analytical 
technique does not handle the wave structure Internal to the Jet (see 
figure 7). 

Theoretical Jet mixing (spreading) effects, as indicated by pitot 
pressure, are much too large at all the axial lo. at ions (this may be 
observed by comparing the theoretical and experimental widths of the 
region of interaction in figure 11). 
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In figure 12, the theoretical data computed using the Eggers* vis- 
cosity model (and Prandtl number of 0.7, 1., and 1.4) are compared to the 
same experimental data given in figure 11. As can be seen, the agreement 
with this viscosity model is about the same as that of the Ferri- 
Kleinstein model. Likewise, the discussion of figure 11 is in general 
true of figure 12. 

The nonreacting free- jet case, resulting from the use of nitrogen 
as the test medium, is presented in figure 13. The theoretical results 
obtained with each of the viscosity models are so clo'''' together that 
only the theoretical results obtained with Eggers* model will be pre- 
sented. In this figure, only the theoretical results obtained with a 
Prandtl number of 1 are offered since this gives the best agreement. At 
an axial location of ^ theoretical and experimental results 

have the same general shapes. The numerical agreement, however, is quite 
poor in the near centerline region y/r^ = ± 1.5. In addition, the shape 
agreement is short lived and iisappecrs by the time an axial location of 
40 is reached. For all values of experimental data have 

a minimum at the centerline and the theoretical have a maximum. That is, 
the theoretical data exhibit a valley in the near centerline region. 
These results indicate that the theoretical near centerline Mach numbers 
are too low, thus pij* ’ucing pitot pressures which are too high. This 
contrary behavior of the theoretical predictions is probably due to im- 
proper handling of the jet wave structure. 

The expected disagreement for the test stream mixing boundary is 
also present. Furthermore, the region of no test stream interaction (for 
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Figure 12.- Experimental and theoretical pitot profiles at various 
axial locations for the reacting free-Jet (Eggers' 
viscosity model). 











Figure 12.- Continued 
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Figure 13.- Continued. 
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Figure 13.- Concluded. 
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example at 1 7. to ± 2.) gives excellent agreement between the 

experimental and theoretical values of pitot pressures* 


Ducted Data (Circular Combustors) 

The experimental data for the ducted case are presented in figures 
14, 15, and 16* Figure 14 presents the pitot surveys made at the exit 
of the four ducts using air test medium (reacting case). Figure 15 gives 
similar data for the nonreacting case (nitrogen test medium) • Figure 16 
gives the static pressures measured along the various ducts. The static 
pressure measurements and exit pitot piofiles were made simultaneously 
for each duct. The same technique was used for both air and nitrogen 
test gas. 

The program was unable to calculate the flow field for even the 
shortest (length * 30 r^) duct, therefore a comparison between the 
experimental and theoretical data cannot be made. The inability of the 
program to compute the flow field for the ducted cases stemmed from the 
fact that the underexp - sion shock wave, which reflects from the duct 
wall, is unable to traverse the region of test stream jet interaction. 

The flow angle computed for the jet and its Interaction with the test 
stream ere inconsistent with the shock wave and the rest of the flow 
field. The reflected shock is not suspect since it is fully compatible 
with the portion of the test stream which has not interacted with the Jet. 


An Evaluation of the Analytical Tool 
The utility of the analytical tool as applied here appears quite 
limited with either of the twc viscosity models employed. This Is not 
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Figure 14.“ Lateral pitot pressure distributions at the exits 
of the four ducts. (Test medium, Air) 
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Figure 15.- Lateral pitot pressure distributions at the exits 
of the four ducts. (Test medium, Nitrogen) 






Figure 16.- Static pressures along the duct wall. 
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to Imply that the Inadequacy of the program is due to the viscosity models 
employed. In fact, the similar results obtained with the two different 
models indicate that the fault lies in handling the internal (jet) flow, 
and its interaction with the test stream. This hypothesis is further 
supported by the fact that the spreading as indicated by theoretical 
pitot pressure is excessive for the reacting case (see figures 11 and 12). 
Similarly, pitot pressure was consistently overpredicted by the program 
for the nonreacting case, particularly in the near centerline region 
which is highly dependent on the jet flow. 

Another failing of the program is that it did not detect the fact 
that the expansion fan is terminated by an intersecting shock front. A 
discussion of why there is an intersecting shock terminating the fan can 
be found in reference 21, and one may refer to the discussion of figures 
7 and 8 for the location of the shock. The program has incorporated in 
it a subprogram which checks the entire flow field at one axial location 
before each downstream marching step is taken, to see if a pressure 
gradient exists. If a pressure gradient of sufficient strength is en- 
countered, it inserts an embedded shock wave. Thus, it must be concluded 
that the computed jet flow did not produce the pressure gradient necessary 
for shock wave insertion. 

It should be noted that all of the shortcomings of the analytical 
approach, detected by the present work, are associated with the divergent 
internal (conical jet) flow and its interaction with the test stream. 

Thus, there is the possibility that it can be successfully applied to the 
case of underexpanded jets whose flow divergence is small at the injector 
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exit. 

It is felc that the analytical predictions can be improved by re- 
placing the present method of handling the jet flow with a characteristic 
expansion network. This network, wnich would require considerable effort 
to implement, would be terminated by an intersecting shock front. The 
intersecting shock may require little additional effort since the sub- 
program previously mentioned may insert the required shock once the proper 
jet flow is computed by means of an expansion network. An improved pro- 
gram, which correctly handles the jet flow by an expansion network, would 
be very us* ful and probably effective in analyzing underexpanded jets. 



CHAPTER V 


CONCLUDING REMARKS 

One of the major technological problems facing the scramjet engine 
concept is the ability to successfully predict the flow field resulting 
from the injection, mi:cing» and combustion of hydrogen fuel. Such pre~ 
dictions are necessarv for good design of major components of the 
engine (i.e. fuel injectors, tht combustor, and the exit nozzle). Of 
particular importance here is the ability to predict the flow field re- 
sulting from underexpanded injection of hydrogen. More fundamental, 
however, is the need for experimental data on an underexpanded H 2 jet 
in a supersonic flow. 

The present work has accomplished the task of furnishing a small 
data base on the coaxial injection of an underexpanded jet into 

supersonic fl.jw* The data obtained are for a Mach 2 test stream of air 
or nitrogen, and a Mach 2 hydrogen jet whose exit pressure is approxi- 
mately twice the test stream static pressure. Since the air or nitrogen 
test stream has a static temperature of 1338 K, data with and without 
combustion is provided. In addition, the facility was operated in a 
free- jet mode and in a ducted mode furnishing data for four different 
duct lengths. The free-jet data consist of radial pitot profiles at 
various axial locations. The ducted data consist of radial pitot pro- 
files at the* duct -ends and static pressures measured along the duct walls. 
In addition, the present work tested the utility of an analyticax tech- 
nique designed to predict the flow field resulting from the injection of 
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an underexpanded jet into supersonic flow. The theory is tested by com- 
paring experimental data with theoretical predictions. The theoretical 
calculations, which cover a wide range of Prandtl number (0.7 to 1*4), 
were unable to correctly predict the experimental results. 
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APPENDIX A 
THEORY 


General Governing Equations 

The basic governing equations are the well known "viscous-inviscid** 
equations used in higher order boundary layer and viscous flow field 
analysis (see references 15, 16, and 17) with the finite rate chemistry 
terms included. These equations are evolved from the full Navier-Stokes 
equations by assuming that the transport effects depend only on gradients 
normal to the stireamlines. (The normal momentum equations are kept in 
the inviscid form.) 

These equations, written in nond imensional form for an intrinsic 
coordinate system (with s along the streamlines and n normal to the 
streamlines), are as follows for axisymmetric flow. 

Gloval Continuity: 

+ pq ^ sin 0 = 0 (Al) 

ds ^ dn y 

S-Momentum : 

I? ^ 


where. 


N-Momentum: 


Re 



Ji 


(cos 0) 9q 1 
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2 30 

P‘1 3^ 


+ 


3n 




(A3) 
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Energy: 


‘“1 8 - ’ ■ ^2 - ^ “l "l 


(AA) 


where. 


1 [l- 


Re (y - 1) -- 

' *00 ' oo 


E^|IjC fi (y„ - 1) „ |l’ 

Pr 9n pj dn “ ' oo ^ \9n 


U C (cos 0) 

P ^ + 

y Pr 3n 



Species Conservation: 


9«i 

IT = ^3^ ”i 


(A5) 


where , 




(hLJL >il 

U Le cos 0 9«i I 

Re 

[an 

V pr 9n / 

y Pr 9n J 


State: 


W^pT 

P = -^1 — (A6) 

Y M^W 

' O 


where , 


w “ 



-1 


For supersonic flow fields, the above equations (A1 to A6) have a 
dual mathematical nature (see reference 15). That is, they exhibit 
features of both hyperbolic and parabolic systems. The analytical tool 
of reference 9, therefore, uses a numerical scheme employing a character- 
istic network in conjunction with a boundary layer type network to yield 
a coupled solution. This scheme is thoroughly discussed in reference 17, 
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and will not be fully covered here. However, the following description 
of the approach used by the scheme is offered. 

Essentially, the approach finds a characteristic solution which 
feels the effects of diffusion and finite rate chemistry. This is 
done by creating the diffusive and chemistry terms as forcing functions 
in the ’V.ompatabllity relation** along characteristics. Treating these 
terms as forcing functions r iults in the characteristic directions of 
the viscous system being exactly those of the inviscid system. Namely, 
the frozen Mach line (C±) 

^ = tan (0 ± (A7) 

and thus the streamlines are defined by the equation 

^ » t<in 0 (A8) 

clx 

The compatibility relation can be shown to be (see reference 17 for 
an excellent derivation) 


sin cos 

yIp 


dp ± d ) + 


sin 0 


(Yf - 1) 


Y P 
*00^ 


S2 + 


iL 

pq 



W. 




- > 

- 1) PI 

I sin 
I cos ± y^) 



Z 


W, 



dx = 0 


(A9) 


The program of reference 9 Is designed Co analyze the mixing and 
combustion of an underexpanded H2 jet; therefore, it is apparent that 
the equations previously presented are not sufficient. Since Che jet is 
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underexpanded, it has an exit pressure greater than the test stream 
static pressure and must expand Into the test stream. Tne expanded jet, 
however, is seen by the test stream as an obstruction and an exit shock 
wave is generated. In addition, embedded shocks caused by combustion 
compression are possible downstream. The equations required to perform 
the expansion and shock calculations are also incorporated into the 
program (reference 9)* 

The expansion was assumed to be isentropic, two dimensional, and 
invlscid in the limit of vanishing radial distance with respect to the 
injector lip. These assumptions allowed the use of the following 
isentropic relations (Prandtl-Meyer expansion) near the injector’s lip. 


i. 

State P/p^ = constant 

(Alt) 

2. 

2 

Energy h + 1/2 V ■ constant 

(All) 

3. 

Momentum ^ + 1/2 d(V^) = o 

P 

(A12) 

4. 

Compatibility — d(inP) ± 0 

Y cos p sin p 

(A13) 


In the case of the shock wave, it was assumed that the chemistry 
was frozen across the. shock and that it was two dimensional (the 2-D 
shock is an exact solution for the conical shock if there is no angle 
attack). Thus, the following Rankine-Hugoniot relations were incor- 


porated 

into the p''ogram. They are: 


1. 

Continuity * P 2 U 2 

(A14) 

2. 

2 2 

Normal Momentum Pj^ + Pj^(U^) “ P 2 P2^^^2^ 

(A15) 

3. 

Tangential Momentum V • V 

'1 2 

(A16) 

4. 

Energy H ■ h + (l/2)v2 - constant 

(A17) 

where, h 

- 1 «j^h^(T) 
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5* State p = p(p>T,<^^) (A18) 

Exit (Underexpansion) Shock 

As previously stated, when the jet expands into the test stream an 
exit ^u^^d^rex;'ansion) shock wave is generated. The idealized flow re- 
sulting from such an interaction is depicted in figure Al. 

Although it can easily be deduced that the pressures (p’s), and 
flow angles (0’s) are aqual on either side of the slip line separating 
the regions 1 and 2 of this figure, it is not p ssible to calculate them 
by a direct ethod. Fortunately, the downstream conditions can be cal- 
culated bv the iterative process that follows. A shock angle is chosen 
(an angle slightly larger than sin ^(1/M^) is a gov^d choice) and the 
downstream properties (P^, 0^, eti.) are computed. The jet is then 

expanded from its exit pressure to the pressure P 2 * Fp* flow 

angle ©2 associated with this pressure does not equal -he flow angle 
0^ downstream of the shock wave, a new shock angle is selected and the 
above procedure is repeated until convergence is obtained. The stream 
properties (p, 0) for which convrrgeno'^^ is obtained are the properties 
existing across the slip line of figure Al. 
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Figure Al.- A schematic of the flow field resulting from the Interaction 
of the test stream and the underexpanded jet^ 
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APPENDIX B 

MODIFIED COMPUTER PROGRAM 

The original program of reference 9 has been streamlined and mod- 
ified to the extent that it is not readily recognized as essentially 
the same program. Numerically, both old and new versions give the 
same marhematical results for the cases they are both able to handle. 
(The original program was not able to handle shock waves which ran from 
the outer boundary toward the centerline, and various other subtleties.) 
The version given here has the Ferri- Kleins tein viscosity model as did 
the or^'ginal of reference 9. 
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Asqrr 

PROGRAM CHARI INPUT* OUTPUT* TAPE7.TAP£SsINPUT,TAPC6=OUTPUT> 

C MAIN PROGRAM FOP CHARACTERISTICS WITH SHEAR 

CONMON/Ae/EPP*EPC*EPT 
COMNON/AC/IBOO*FIN 
COMHON/AI/GAR*GEM 
CONHOK/AX/JSUSL * JSU3U 
COMNON/OA/A(.P(7*$St *EHINF.UlNF 
COMMON/eC/IOCHEM 
COMHON/eC/XHASSISGI 

C0MH0N/CG/AUP.eUP*CUP*0TSPRI(55l*0UP*EUP«JC0NV*THPRIIS5l*VPRI 1551 

COMNON/CK/UTHOLE 171 

coHNON/oe/PETet«ii*is<u) 

CONNOK/DE/MM 

CCMMON/OP/TNI55I 

COMMON/EO/CPIN*RC 

C0HN0N/EF/EMf55l*GAMI55)*PI55l .TH<55>*V(55) 

C0MH0K/EG/EIN*PR*KLE 
COMHON/EF/GAMINF.HI (7I.PINF 
COMMON/FE/OEL 

COMHON/GE/RAQ.ROC.UIN*VISINF 
COMHON/GF/0£LV.CVISA.KOUNTO*VISA 
CCMNOK/GK/OEi.X 
COMHON/HJ/KOUNT *LL*NPT 
CONNON/HL/ACPKA .BETA 

COMNON/HP/ At. PNt 7.55)* CPNI 7*551 .CPXNI55I *EMN 155) *G AHNI5S I *HN (7 *551 * 
IL.PNI55) •GNI55 I *RhON( 55)*RN(55) .THNI55) *TNI55i *WN<55) ,XhUN(55 ) 
CCMM0N7 HP/B|£TANI<i)*IEHSE 0 
CCHM0N/CF/AtPe(7) ,PHI 155) 

C0PMCN/FC/M(55)*XI55) 

COHHOP/PC/ALPHN 17 ) * IFUEL, PRES 
CCMMON/PO/JCHEM.NSP.T 155) 

COPHON/QA/.'<lir*55)*Q«55)*RHO(55l *XMUI55) 

CO)*MON/QS/RHOPI?)*NOOTf 7.55) .M00TCI7) .HPI2I .XHUPI2) 
COHHON/RC/APO.AP1.AP2 

commok/rs/gphs.ps.ths.thsl.thso 

C0)«M0N/ST/I13*I»FGI.K .XFIRST.KKKQ.PSTAR 

CO) 1 MON/TS/OVlSe.OVISr*IFS.MHM.V ISO. vise 

COPMON/UV/IIll.IERR.IPRESS.IPRESU.ISUe 

C0HH0N/WT/0ACM7.55) . 0TCHI55) .DVISO.VISO 

COMNON/tfN/ICCNT.IENO.KT.tHePN.xePN 

COMHON/MV/NPTS.RE.XOP.XJ 

COMMON/XG/XOO 

COHNOP/VX/ABOOS .EP9ESS.CPRESS 
C0MN0K/Y2/aPRESU.CHEMFC*CPRESU.tMSUe.RTH.XSTEP 
CO)<NON/ZV7AQOO*eeCO.CeOO.EeOO.FeoC.G300.IAVE.TPUNCH.J800.KKKKK 
OIMENSIOA XSI7)*PITOTI55) 

DATA 1111/0/ 

begin inputting parameters 





r>0000000« 
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MRtTE(6«400t 

112 F0rtMATa5*5X*9E13.S) 

406 FORNATtlHl) 

J«0 TNC OINENSIONAL 
J*i AXISVRHETftlC 
SPECIES 1 IS H 
SPECIES 2 IS 0 
SPECIES 3 IS H20 
SPECIES 4 IS H2 
SPECIES S IS 02 
SPECIES 6 IS Oh 
SPECIES 7 IS N2 
MTNOLEID^l.OOO 
MTP0LE(2I«16. 
wTHOLEisisia.oie 
HTH0LE(4>32.916 
MTMOLE(S)s32.0 
MTHOLEI6M17.00A 
MTNOLEir|x20.014 
PAS«MTnCLet4)/16. 

J22*0 

XOOsO. 

IOG>0 

DEL>0. 

00 0220 1*1.4 
OEtANtltsO. 

BEI8m«0. 

ism*o 

0220 CONTINUE 
IFS«0 
NSP=T 
R0«1.902 

R00«R0*3. 007*32. 2/2. 20 S«100C. 

EPP*1.E-10 

EPTH«1.E-10 

EPQ«1.E«1C 

EPT«1.E-10 

I0CHEN«1 

HN*6 

EXXXsl.E'CO 
113*0 
XFtRST»-l 
KKKQ*10000 
JCONV«0 
1NPTSH«0 
363 CONTINUE 
CALL INOATA 

C •«••*•••••••••••••••««•«•**«**••••» hake initial SHEAR 

VISO*XVIS(XBP> 
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CFF=0. 

CALL SHCARKCFF.VISQI 
67 S9 CCNTIKUE 

CALL EMOCC 

W’S=»»IS0=!WISC=VISe*WISA=X¥IS<X8PI 

OVISO=OV1SC=OVISB=OFISA=0. 

OC 71B8 X=l,NPfS 
00 7J99 J=1,NSP 

7199 XSIJ>=HIKI«ALP(J*Kl/HTNOLE(Jl 

FUAIR=1.00E*tXSCll«2.«XSf4l«2.*XS(3)*XSI6n/ll6.*(XSC2l«XS(3) 

l«2.*XSI6)*XS<6)>«28.01<i>XS(7l) 

PHIfKIsFOAIR/. 029161 
7108 CONTINUE 

IFIKOUNT. :Q.KKKKK) GO TO 437 
IFIIIll.EC.il GO TO 407 
IFIKOUNT. EQ.KOUNTOI GO TO 4p7 
IFI(|KCUNT/LLI*LLI.NE.K0UNTIG0 to 179 
437 WRITE 16. 4G8IK0UNT 
4C8 FORMAT <7H1K0UNT= 151 
MRITEI6.5206I XII) 

52C6 F0PMATI75H X = E13.5) 

DO 5455 123=1.4 

IFIBETBII23) .EO.O.) GO TO 5485 
IFII23.LT.il 

IHRITE 16.84841 I23.BETBII23I 

8484 FORHATI/5X.20HEIieEODEO SMOCK TYPE . II .1CX.7H8ET A = .Ell. 3) 
1FII23.GE.3) 

1WRITEI6. 23311 I23.aET3II23) 

2331 FORNATI/14X.11HSMOCK TYPE .11. 10X.7HBETA s .Ell. 3) 

8485 CONTINUE 

V1SM*WISA*VISINF 
MRITEI6. 72221 VISM 

7222 F0RHATI/5X.11MVISC0STTY =El3.5,l5H ILB*SEC/FT**2I I 
MRITF|6.5eB7l 

52C7 F0RHATt/3X.3HPT..llX.lMV.l?X.lHQ,12X.lHT,12X.lHP.llX.2HTH.llX, 
12HEMHX3HRHO10X3PGAN9X5HPIT0TI 
no 73 I=1,NPTS 
P|II=Pai/PIN 

70 P1T0TI1I=PCII*PRES*H6BNII)*1.» *.5*EHII)**2I**CGANIII/(GAMIII-1.II 
a*l2.*GAHII)/(GAMn*l.)*EHm**2-CCGANIII-l.l/(GANIll*l.m**ll./l 
ei.-GAMlIIII 

WRITE 16.112) (I. VI I) .Otll.TII) .Pm.THIlI.EHlII.RHOlII.GANIII.PITOT 
HI) .I*1»NPTSI 
DO 71 I*1,NPTS 

71 PII)=PIII*P1N 
WPITEI6.160) 

16C F0RHATI/////3X.3HPT..8X. 6HALPI 1) .7X. 6HALPI2) .7X. 

26HALP I3).7X,6HALPI4) .7X.6HALPI5I .7X.6HALPI6).7X.6HALPI7) .9X.3HPHI 
2.11X.1HWI 
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MRITEf6«112MI. (ALPCJ.l), .PHim .M< II .Isi.NPTSI 

179 CONTIKUE 

IFtKOUNT.GE.KKKKKIOO TO 1S72 
IFdlll.EC.ll GO TO I57e 

ALPHAxl.O 

BETA^C.O 

CALL STFPfVISI 

IFdlll.EC.ll GO TO GC7 

IF(K0UNT.NE.KFIPST.0R.I13.NE.1I GC TO TOO 
CALL PUNCF 

00 331 I=JSUBL*JSUBU 
Y PRIdMY III 
THPRldl=THdl 
3E1 CONTINUE 
3C0 CONTINUE 

CALL CHENIFASI 
8203 CONTINUE 
ICONT=0 
IEN0=t 
K=1 
L=2 

887 IFIL.CE.JSUSL.AND.L.LE.JSU9UI GO TO 9C0 

888 K=L 
KT=L-1 

IFIL.EQ.NFTSIGOTO 612 

IFIL.EQ.JSU8U.AN0.IC0NT.CQ.il GO T-0 1622 

IFI8ETefll.EQ.O..ANO.8FT0l2l .EO .0 . . ANO. EE T6 (31 .EQ.O. . AND. 8ET8 <41 . 
AEQ.O.IGO TO 777 

IF(9ET8(ll.GT.O..OR,BCTei3l .GT. 0. I 777.776 
777 IF(K.EC.IS(2I.0R.K.EQ.IS(4I IGO TC 0236 
IF IK. NE. IS 111 -II 8231, 775 
8231 IF(K.Ne.IS(3l-ll02l4.773 
776 IFIK.NE.ISIII-IIGO TO 11 
775 NHPsl 
K=IS(1) 

KT=K 

GO TO 8232 

11 1F(K.NE.IS(3I-1IGC to 22 
773 MHN=3 
K»IS(3I 
KT=K 

GO TO 8232 

22 IFIK.NE.ISI2I IGO TO 33 
NNN=2 

GO TO 8232 

33 IFIK.NE.ISIAIIGO TO 44 
NHHs4 

GO TO 8232 

44 IF(K.Ea.IS(2l»1.0R.K.Ea.IS(4t«lia8888.8234 
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98888 K=K«1 
L=K 

KT*t-l 
GO TO 8239 
8232 L=K 
8230 IFS=1 
IPOIsl 
ALSVs ALPHA 
8ESV=BETA 

IF(aETam.EQ.:..ANO.BET8l2l.eQ.O..AN0.8ETR(3l.£Q.O..AN0.8ETe(9l. 
AEQ.O.IGC TO 722 

IFt8ETBIlt.GT.O..OR.8CTei3l .GT.C.IGO TO 772 
KTSAW=KT 
772 call CPOIHT 
THOE-THKIKI 
ALPHAS. S 
BETA=.5 

IFtaET8a».EQ.3..ANO.BET8f2I.EQ.a..AN0.8ETP(3l.£Q.O..AN0.8ET8<4l. 
AEQ.O.tCO TO 2199 

IFtGETblll .GT.0..CR.8ETB(3I .GT.C.IGO TU 2199 
K'-=KTSAV 

21% UAiL -?r.T9T 
I-t ^ ' 1 

IF.»F .LT.20IG0 TO 2195 
MHITE <6.9191 I 
MRITC<6«2196» 

2196 FCRHAT(99H ERROR IN CPOINT ITERATION FOR SHOCK IN CHARI 
STOP 

2195 ERTHO»AeS<THOE-THN(KII 
THOE^THNIKI 

IF<8ETe<ll.EQ.O..ANO.eET8<2l.EQ.>]..AN0.8ETe(3).EQ.O..ANO.BET6l9l . 
AEQ.O.IGO TO 771 

IF<8ET8<ll.GT.0.*0R.BET8<3l .GT.O. IGO TO 771 
KT=MSAW 

771 IFIERTHO.GT.EXXXIGO TO 2199 
ALPHAsALSV 
8ETA=BESW 

IF<K.EQ.IS<^n NNH=2 
IF<K.EQ.IS<9ll HHMs9 
IFS»2 

CALL HSHCCKTHHNI 
IFS=0 
K»K*1 
L«L*1 
GO TO 687 
8239 CONTINUE 
IPOI=l 
ALSVsALPHA 
BESVsBETA 
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CALL CPOINT 

KT=L-1 

TH0E*THNIK) 

ALPHAS. 5 
BETAS. 5 

26C1 CALL CPCINT 
KT»L-1 
IPOI*IPOI»l 

IFIIPOI.lt. 20)GO TO ?f>A2 
HRITEI6.B191) 

9191 FORMAT I1H1» 

HRITEI6.2197) K 

2197 F0RHATC53H ERROR IN STANDARD CPOINT ITERATION IN CHAR AT POINT 121 
STOP 

26C2 EPTHjsA8S(TH0E>THNIK) I 
THOEsTMNTK) 

IFIERThO.GT.EXXXI GO TO 2601 

ALPHAsALSV 

9ETA«e£SV 

C ••••»««•••••••••••»••••• INCREPEAT counters 00 NEXT C POINT 

9CC CONTINUE 
K*K*1 

IFtL.EO.NPTSI GO TO 7676 
L*L*1 

IF(ICCNT.EQ.l) GO TO 666 
GO TO 667 

C nozzle NALL calculation 

612 continue 
IPOI sl 
ALSW* ALPHA 
BESVsBETA 

call LPCINT CNPTS.l.l 

KsNPTS 

THOe»THN(K) 

IF( IPRESU.EQ.CI ThOEsPNIKI 

ALPHAS. 5 
BETAS. 5 

26C7 CALL LPCIAHNPTS.1.1 
KsNPTS 
IPOIsIPOI*! 

IFIIPOI.lt. 20)GO TO 2686 
MRITE(6. 91911 
XRITE (6*2196) 

2196 FOPNAKSIH ERROR IN NOZZLE WALL CALCULATION ITERATION IN CHAR) 

STOP 

2608 ERTH0sA6S(TH0E-THN(K) ) 

IF( I^RESU.EQ.O) eRTHO»A0S(l.-THOE7PN(<) ) 

THOEsTHNiK) 

1F( IPPESU.EQ.C) TMCE»?NCK) 
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tFIERlHO.GT.EXXXI GO TO 2G07 

ALPHAsAlSV 

BETAsBESV 

C COMPLETE FIRST POINT 
7676 CONTINUE 

IFtJSUBL.EQ.il GO TO 18CC 

CALL LPCIKTtl.O.I 

K»1 

IPOIsl 
ALSV* alpha 
BESV sBETA 
THOEsTHNtK) 

IFIIPPESS.EQ.C I TNDEsPNtKI 

ALPHAS. S 
BETAS. 5 

2600 CALL LPCINTll.O.) 

X*1 

IP0I*1P01»1 

IFtIPOI.LT.?OlGO TO 2610 
WRITE <6.91911 
HRITE<6.2199> 

2199 FORNATI39H ERROR IN FIRST POINT ITERATION IN CHARI 
STOP 

2610 ERTHOsABStTHOE-THNUI I 

IFtlPRESS.EO.C I £RTHDsAeS(l.-THDE7PN(Kll 

THOEsTHNIKI 

IFtIPPESS.EQ.O I TH0E=PNtX» 

IFIERTHD.CT.EXXXI GO TO 2609 
ALPHAsALSU 
BETAsBESV 

C SUBSONIC PRESSURE ITERATION 

1600 CONTINUE 

IFtlSUB.EC.CI GO TO 1622 
IFtICCNT.EQ.il GO TO 162? 

IFtllS.NE.l. OR. FOUNT. NE.KFIRST) GO TO 1777 

CALL OPOTHtTHS.JSUBUI 

CALL OPOTHtTHSU.jSUaU^ll 

CALL CPCTHtTHSL.JSUBU'll 

CALL THSSStTHSSI 

AUPsVtJSUBU) 

BUPsTANtTHtjSUeui I 

IFtIREGI.NE.Q.ANC.JSUau.EQ.J22l GO TO 6375 
CUPsTHS/CCStTHt JSU8UI l**3 

OUPstTHSS^S.’TANtTHtJSUBUll *THS*THSI/CCS<TH(JSU9UI l^’N 
EUPs.u.095 
GO TO 6376 
6375 CONTINUE 
CUPsQ. 

OUPs-1. 
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EUP=0. 

8376 CONTINUE 

00 381 I=jSUOL*JSUflU 
381 CALL CPCTHIOTSPRICII.I) 

1777 CONTINUE 

CALL SSCNICaOG) 

I«=’linG.EC.OI GO TO 16?2 
1111=1 
IPUNCH=1 
GO TO <•'57 
1672 CONTINUE 

IFdSue.EC.CI GO TO 359 
IMJCCNV.EQ.l) GO TO 36C 
lFCK'OUNT,NP.KKKQ*n GO TO 361 
IFlINPTSH.eO.CI nptsh=npts 
INPTSP=J 
REWINC 7 
:: TO 363 
36C NPTS=NPTSH 
iNPTSHxC 

359 OC 357 I=1*NPTS 
II = NPTS-M1 

IF (EMNfll) ,GT.E«SU3) GO TO 357 
IFIISUB.EQ.II GO TO 355 
ISUfl=l 

WRITE (6*364) 

354 FOPMATC37H1 SU3S0NIC REGION ENCCUNTEREO) 

EMST=EMSUP 

FHSUB=1.15 

AP1=0. 

IREGI=0 
GO TO 359 

355 K^IINI 

IF (IREGI.EQ*!) K= JSURU 
CO TO 358 

357 CONTINUE 
IREGI=2 
GO TO 361 

358 CONTINUE 

IFIJCCNV.EQ.O ) GO TO 1417 
JCCNV=0 

OC 1418 I=1*NPTS 
II=NPTS-1N1 

ir(EHN(II).GT.ENSTI GO TO 1418 
GO ro 1417 
1418 CONTINUE 
IREGI=2 
113=0 

EMSU0=EMST 
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XOCsXPPN 

ABOCsVNil) 

aecos=vMi» 

BeOOsTAMTHNd)) 

ceco=c. 

GO TO 361 
Hiir CONTINU6 
APCsPNUJ 

AP2= (PMK»*PN tl)>APl*IVN(K»-VN(l))»/ tVN(K>-VN(l> )**^ 

00 356 IsX.K 

PN<n=APO«APl*<fNII)-VNm » ♦AP2*<lfNm-YA<l H ••a 
9M0N < I» = GEM*WMn *PN< II / TNdl 
IFII.6C.U GO T- 355 
XJ1=1.»XJ 

YFUMs tYN<l)*»XJl-t*Mll)»»XJl»/XJl 

TeRK=(RHCMl) •QMl)*C0S«THNmi ♦PHCN < II I »0N ( H> *COS (THNIIim /2. 
XMA5S(I)*XHA5SII1I»TE9H*VFUN 
356 CCNfUUE 

0S=2.*DFtX/ <COS<IMI<I l♦COS<THN<KI)l 
PS*(“N<XI-P<<M/OS 

GPHS= <CAPK<<! •PN<KI*FHN<KI**2-GAH<KI*PlXI*fM<KI ••21 /OS 
ns=2 .*0ELX/<C0S(T*-N<1 II*C0S(TH<1M I 
361 CONTINUE 
C COMPUTE shear 

VISE*xviS<xnPfri 

CFF=0. 

CALL SHEAP2<CPF,VISE) 

(• RESET ALPHA AND BETA 

IF<IAVE.EQ.O» GO TO A306 
IF<BeTA.6T.C.CI&0 TO 83^6 
ALPhA=0.5 
BETA=C.5 
GO TO 8282 
8396 CONTINUE 

C ••••••••••••••••••••••••••••••••••••••••••• STEP TAKEN OUTPUT 

J22=JSueC 
call PSET 
00 1431 I=l<4 
1431 BET8m=eETAN<II 
K0UNT=KCLNT*1 
GO TO 6785 

1572 IFIIPUNCH.EQ.OICALL EXIT 
CALL PUNCH 
CALL EXIT 
END 

SUPROLTIKE ENBEO 

COHHON/AC/IBOO.PIN 

CCKMON/AL/GARtGEM 
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COt'N ONV A X/ J St PL . JSU>1U 
COt*MON/eA/ALP(7,55l »EHlNF,HINf 
COHHOK/8C/XMASS <551 
Cr>MM0N/CJ/CP<7,55» ,CPl(7» ,CPK<55I 
COPMOt/CK/MTPOLE <71 
C0PM0N/0B/BETe<4» ,IS<4» 

CCHMON/EC/CPIN,RC 

C0HM0K/EF/E*<<55»,GAM<55I,P(55I .TH<55I ,V<55» 
nCHMON/EF/GAMINF,Hl <71 ,RINP 
COPMON/GK/OELX 
CCphOK/HL/ALPHA .BETA 

CCPMON/HP/ALPN< 7,55I.CPN<7.55» ,CPXN<55» ,£MH<5P| .'■.aNN(55».HN<7,55» , 
IL»PN <551 ,CN<55» ,ReON<55J . =N < 55 » , TMN< 55 1 t TN < 55 » i MN < 551 ,XMUN<55 J 
COMMON/hF/B'’TAN<m,IFPBPO 
COMM OK/PC/ W< 551 ♦X<55l 
COMMON/PC/JCHE1,MSP,T<55I 
CCHMON/QA/M<7t55l ♦ 1 <5 5 1 tPHT <55 » ,XMU<55 I 
COMMON/PC/P <55 » 

C0f'M0t/U/0r!<55l ,CAlp <‘.551 ,OBO < 551 . OCPX <55 I ,00 At P < 7. 55 1 , 3T AU <55 1 , 
1TAU<55I 

C0MM0N/hW/NPTS,PF,X5P,XJ 
DATA EPPPES/1 .E-Ct/ 

IEMeEC=: 

AtPHA = l . 

PETA=C. 

AtPH=l. 

P€T=C. 

0ELX=1, 

DO 53C M=l,2 

IF<IS<MI ,KE.O> GC TO 50C 

IM=MPTS-2 

00 1 1=2, IM 

if<i.ge.jsuol.anc.i.lt.jsu5ui go to 1 
12=1 
11=12-1 
13=12*1 
I<.= 12*2 

TlC=V<IH-r <I2» 

Tll=r<I2l-r<I3) 

T12 = Y<I3I-Y<I<.I 

lF<T10.LT.l.e-C4.CP.ril.LT.l.E-C«..OP.Tl2.LT.l.E-0-,l GO TO I 
0Z=Y<I*1I -Y<I» 

IF<<M72I*2.E3.M| GO TO 200 

XP2 = XM1<ALPH ,BET , TH < I ♦! I , XMU ( I* 1 1 , 0 , , C . I 
XP1=XMKALPH ,BcT , TH < I »,XMU<I ),0,,C.I 
GO TO 201 
22P CONTINUF 

XP2=XM2<ALPH ,PET , TH < I *1 ) , XMU < I * 1 1 ,0 , , C . I 
XPl=y*'2<AtRH ,9ET ,TH<I »,XMU<I ),i).,S.I 
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201 0ZLAH=KP1-XP2 

IFIOZtAP.LT.l.e-lCI GO TO 1 

OI=nz/CZLAM 

IMOII 1,1,7 

7 IMOI,CT.10.*OELK» GO TO I 

pis=piiii*pmi 

P2S=Pa2»*P<I2l 

P7S=P<I3»*PTI3» 

P 4 S=p(i 4 i*pati 

T1=P<11I-P(I2I 

T2=P«I2)-FCI3I 

T3=PJ I3l-P<14) 

Ti»=PlS-P2S 

T5=P2S-P3S 

T6=P35-P45 

Tr=PlS*P(Ill-P2S*P< 121 
T6=P2S»P<12I-P3S*P(I3» 
T9=P3S»P<I3»-P4S*PII4» 


CAU 

SOLVE (T1 

• T2 

.T3 

.T:. ,T5 

.Tb .T7 

.18 


CALL 

SOLVE (TlO 

♦ Tn 

.T12 

.T4 ,TS 

.T6 .T7 

• T8 .T9 

.091 

CALL 

SCLVEITl 

*T? 

.T3 

tUO. Hi 

.T12.T7 

.Tfl ,T9 

• DC I 

CALL 

SOLVETTl 

.12 

.T3 

• Ti. »T5 

♦T6 .TlO 

♦ TU.T12 

• 001 


9=08/E 

C=OC/f 

0=00/E 

A = YtIH*Pai»*<-B*P«I II *<-C* 0 *P(IlM I 
TRES 1 ./ 3 . 

COS 1 ./ 27 . 

YST=A-C*e*TRE/D*2.*C *3*C0/0**2 

IFUST.LE.Ydl .OR.YST.Gf.Y(I*lM GO TO 1 

YSTP=e-C*C*TRk /0 

IFIYSTP.GE.EPPRESI GO TO 1 

ISCMI =!♦! 

IF<(H/ 2 I* 2 ,EQ.MI IS«MI=I 
¥P 5 =*P 1 
XP 6 *XP 2 
GO TO 501 
1 CONTINUE 
GO TO 500 
501 ISM*IS<NI 

eETB(NM(ATAN(XP 5 l*ATANIXP 6 ) 1 / 2 . 

MRITEI 6 , 506 IM 

506 FORNAT« 1 M 1 , 20 X, 19 HENBEOOEO SHOCK TYPE I 2 // 13 X, 2 MIS , 5 x, 4 mbETA ) 
HRITE( 6 , 5 C 0 I IS(N),BET 9 (HI 
5 CB FCRHAT( 10 X,I 5 ,E 11 . 3 I 
L =1 

IFnH/ 2 »* 2 .EQ.HI L = -l 
ISNM» ISP-L 
ISP«lSN*t 
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RAT= ( (YJISMI ♦ •( ! /2 .-Y ( ISP »> / ( Y ( I SM l-Y < ISP» » 


X 

-X 

(ISPI hPAT*(X 

(ISHI-X 

(ISP) ) 

Y 

<ISH)=Y 

CISPMPAT*(Y 

< ISHI-Y 

(ISP)) 

Q 

<tsh: =0 

(TSPMPAT*(Q 

(ISHI-Q 

(ISP)) 

P 

dsf') --P 

(ii>PI^PAT*(P 

(ISHI *P 

(ISP)) 

T 

<ISH>=T 

(ISP) ♦94T»( 1 

( ISHI-T 

(ISP)) 

TH 

<ISM = TH 

(ir'>)^PAT*(TM 

(ISHI-TH 

(ISP)) 

BO 

CISM) “BQ 

(ISP) ♦PAT*(BO 

(ISMl-80 

(ISP)) 

TAU 

CI?M)=TAU 

(ISP)^PAT^(TA J 

(ISMI-TAU 

(ISP)) 

C0O 

nsH) = nBQ 

( ISP) ♦PAT*(OeO 

USHI-08Q 

(ispir 

OCPX 


ISP) ♦PAT» (OCPX 

(ISHI-DCPX 

(ISP)) 

OTAU 

nsH» ' 

(ISP)^PAT*(0TAU 

(ISHI-OTAU 

(ISP) ) 

XMASS 

<ISH ' ' 

( ISO) ♦RAT*(XMASS 

1 ISMI-XHASS 

(ISP) ) 


CPK< I$K»=C. 

W<I3fM=C. 

CALL THER‘.'0<MISH|,«i .CPU 
00 IOC Kl-l.NSP 
J = KI 

ALP <KI,ISM|=ALP (KI.ISPI *OAT*(ALP (Kl.ISPI-ALP (KI.ISPI) 
DALP (KI.ISHI^OALP tKi.rSPI ♦RAT’IPA lP (K I • I SHI -OALP (KI«ISPII 
OOALPIKI, ISMI =OOALP<KI. ISP| ♦ RAT • < COALP < K * . I SHI -OCALP< KI, I SP 1 1 
HI J, ISMMU (Jl 
CP<J.1SHI=CPUJI 

WIISMI ^WUSHI*ALP(J.ISHI/HT‘10LeUI 
CPXC ISHI-CPxnSHI ALP U. ISHI •CPIJ.ISMI 
P N<KI,1SHI=H (KI.ISHI 
ALTN|t<I,ISMI=ALP(Kl,ISHI 
too CONTIEUt 

H(I3HI=l./MtlSMI 
R«ISHl=RC 'wnSHI 

GAH(ISMI=CPX(ISHJ/(CPX(ISM| -RIISHI 'CPIM 
RHO< tSHI=P(ISMI *W<ISHI»&EH/TaSH| 

Rr=l./RUSMI 

EHdSPI -C<ISPI*EHINF*SORT (GAR/GAH (ISHI*RI/T IISMI I 

XHuaSHI = ZHUIEM(ISMI I 

Q N(ISPI=Q (ISHI 

R N<ISPI=R (ISPI 

T N»rspi=T <rsMi 

p Naspi=p (ISHI 

M NdSHIsH (ISHI 

TH N(ISHI=TH (ISHI 

RHC NaSHI=PHO ISHI 

GAH N(ISHI=OAH(ISHI 

IF(JSLei.GT.IS(MI I JSU8L = JSU0LH 

lF»JSl;eU.CT.IS(HI I JSU8U=JSUBU*1 

00 101 XK=1,4 

IFEIS(KKI .6T.ISCHII I S (KK I = IS « KK I ♦ 1 
101 CONTINUE 
lEMBEO^l 


0MGDtALTA«EB 
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CALL HSHOCK(H) 

V ; CISMI 

itHAiS aSMHI r XHftSS as: 

W MISKNI 

P «ISH*»I = P 
C aSHPMQ 
T fISHHIrT 

TM IISPP)=TM mispmi 

£M 

RhC lISHWMPMONt iSMMf 
CPK i I SHPI -CPXN ( 

GAH(ISH*<I=GAHM ISMM) 

XMUl TSWMI =XHUN 
no 1313 Xlrlt^SP 

H N(KI*IS^M) 

CP (Kl.IS^'fll^CP 

fttPCKl, =ALPN (KT* ISP*«» 

1313 C0NTIMJ£ 
iPOo=r. 

P£ = 3. 

xep=o. 

CALL SHeALin,,n.) 

sen coktikue 

PElUPr 

ENO 

SuePOUTlNE SOLVE (All « Al3* Al 3 * A?l * A?2 t A2 3 ♦ A3 1 • A3^ • A33# OE T I 
0£T= All^f Ac2*A33-A3?^A23I-A1c*( A2i»A 33-A31»A23MA13*( A2l*A32*A22*a 
131) 

PETUPK 

END 

SUPROUTINE HShCCxIK) 

C0PH0K/AL/GAQ*0EW 

COPrtON/BA/ ^P(7,S5) #E**lNf ,»(INF 

C0PH0f/9e/SlA*S?A»S3AT 

C0MHnf4/pr/SAHP,PP.''6*RHCR,THP,WetXHU3,V^ 

COPHON/CA/WOOTN(7,S5) •XNfSS) 

C0PH0r:/Cj/CP(7*SS) *CP1(7) tCRXISSi 
COH«Cf/Oe/BETe«A) •ISCAI 
CO^^HOK/LPXYNf 5S) 

COP«CN/EC/CPlN*PO 

CO*»HON/EF/EHfSS) *GArt(55) »P( SS I • TH (55) » Y (5S ) 

CCWO^/EG/PIN*PP#XL£ 

COMHON/EF/GAHINF*Hl (7) *9 INF 
C^HHOK/fE/Ofl 
w-^KMON/GK/OELX 
COPMOK/WL/ ALPHA, BETA 
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C0HNQM/H|i/ALPN(7,SSI.CPN(7.55I .CPXNISSI »EHN ($51 .GANNtSS I ,HN( 7«5SI , 
ILS.PN<55I.QN(55I,RH0M(55).RN(55>«THN(55I.TN(55I ,MN(55I«XNUN(55I 
C0HNQN/HP/8PTAN(A;* ..^PSEb 
COHHOK/FO/jCHEft.NSP.T(55l 
COMHON/OA/H(7,55l .0(551 «PH0(5St ,XHU(55I 
CCNH0N/CS/RH0P(2) .400T(7«!>5I .ND0TC(7) .UPI2| .XNUP(2I 
C0HN0N/SC/8QN(55I .0ALPN(7.55I .0BQN(55> .0CPXN(55t .DQALPN(7«55) . 
tOTAUN(55I.TAUN(55l 

CONNO^/ST/I13.IRE(I*KS.KFIRST.KKKQ.PSTAR 
C0P8CM/»S/0VlS8,Ctf ISC.IFS.HMM.WIS8.VISC 

COPNON/TU/BQ(5Si.3ALP(7.55l.08Q(55l.OCPX(5SI.OOALP(7.55I.OTAU(55l , 
1TAUI55» 

C0P*(0K/TV/AI.PP(7.3I ,Pf T .BQPIZI . OACHP ( 7.21 ,0 ALPPC 7, 2 • .080P ( Z I . 
10CPXP(2I.OOALPP(7.2I.OTAUP(2) . DTCHP(2 1 .G^HP (2 1 .PP ( 2 1 ,QP(2 1 . 
2TAUPI2I.THPI2».Tf‘(2».VP(2l 
COHNOb/VM/ICONT .tEN0.KT.THBPN.X8PN 
OIHENGICN 00HM(7» 

I=ISIKI 
IT1=! 

8ET=(JETBIk| 

IT11=1 
XXX=l. 

ICCC=C 
I =1 

rF'(K/2l*2.NE.KI L=-l 
NiIS(KI ♦L 

IFCBETA.NE.O. I 00 TO )>210 
TAUN(P)=TAU(H> 

BQNINMSCINI 
OCPXN(PlsCCPX(HI 
0 . AUN(HI=OTAU(NI 
OBQN(HI=OeO(NI 
CPXMNI=CPX(NI 
00 8211 J=1.NSP 
DAlPN(J.KI=OALP(J.NI 
00ALPM(J.NI=0DAEP(J.N> 

A'.PNt J,N» = ALP» ■ 1.I» 

HCOTCIJl^O. 

HOOTNIJ.HI^O. 

HN(NI=MN(II 
8211 CONTIFUe 
8210 CONTINUE 

IF(8ETA,GT.0.»eeT=BETAN(KI 
4 IT=1 
CA;-1. 

5A=0. 

VT=ON(It •COS(BET-THN(in 
UlsQNdl 
U1<A8S(U1I 


•CA*SIN(8ET-THNT1M 
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XMSsRMOKni *U1 
GNsGAKMH 
GP1= IGN«1 .1 
6W1=GK-1. 

RN1=1./RN(I> 

XMlsUl**2*EMINF»*2»t6aR/GAMMI»*RNI/TNtH I 
0XH=1./XM1 

IP<It,£0.nU2=L'l*«(iMl»XMl*2. >/6Pl*0XM 
5 »M2PaXt*S/C2 

P2H=XHS*ll;l-U2»»PN(I) 

V2=VT*»2 

V1=V2*U1*»2 

V2=V2»U2**2 

H6=0. 

00 14C0 JsltNSP 
1I.CC He=HN(j,ll*ALPNtJ.I)*M6 
H2sH6*IW1-V2I/2.*EIN 
im=i 
Tl = TNtn 

lFlini.EC.l»T2=Tl*{2.*GN*XMI-6Pl)*«CHl»XHl*2.) /(2.*GP1)*0XM 
82CC CALL TH€PK0tT2,Hl,CPll 
H 2 PaO. 

00 92C1 Jsl.NSP 
92C1 M2P=H2F*ALPNIJ,m *H1 ( Jl 
ERRl* (H2-H2P>/H6 

IMAeStERRlI .LT.l.E-QSI GO TO 8202 
IITlsIlTlU 

IF«1IT1.GT.15» GC TO 8203 
IF(IIT1.CT.2) GC TO 8204 
ERP 2 aER 4 ‘l 
T22=T2 
T2=T2*1.C1 
GO TO 82tC 
8203 MPITEC6,9iqi» 

9191 FORMATTIHI) 

WRITEiG.820S) 

8205 FORMATf* ERROR IN TEMPERATURE LOOP HSMOCK*» 

STOP 

82C4 0UMaT22-ERR2*«T2-T22l / <ERR1-E«»R2» 

ERR2sERRl 
T22=T2 
T2=0UP 
CO TO 82GC 
8202 CONTINUE 

RH2aP2H*NMM»*6en/T2 
ER=«RF2-Rm2P»/RH0 id 
IFIABStER».LT.l.E-8»GO TO 7 
ITaIT*l 

IF(IT.GT.15.IGO TO 100 
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IFIIT,GT.2K>C TO 6 

U?=U?**99 
GO TO 5 

ICO WPITE«e*9l91l 
NPITFie.cCOI 

2CC FCRPftTC* ERROR IN HUGCNIOT LCOP IN HSHOCK*) 

STrp 

6 0UH=U22-EP?*«U2-U22I/ IEP-EP2) 

FP?=£R 

U2?=U2 

'J2-CUf» 

GO TO 5 

7 continue 

C8=C0?IGET) 

SP=SIM9FTI 

IF (fK/21*: .EQ.K) U2=-U2 

QN2P=-U2»CA 

UV=VT^Cfl-CN2P»S9 

UV-VT«S8tON2P*CP 

PMF2 = ATAN f«V/UWI 

G2=SQPT <UV*UV^MV»MVI 

IFlI'^PBEC.FQ.n GC TO T935 

TNIMI =YIM) ♦ .S*(TaNf0ET9fK) I ♦ TAN I8E T I) •OEL* 

OEt-1. 

IF IBEief II •E0.a..AN0.P£TPl2l .EQ*C ..ANO.BFTBf 3>.EQ*a ..ANO.eETBUJ . 
AEQ.O.IGC TO 777 

IMBETeCl) •GT.O..CR.0ETBC3I .GT.C.I 7T7,77e 
777 CALL LPCINTfH*!.! 

GC TO 775 

776 TMMNH/2*2.EQ.PMMICALL LP0INT|N*G .1 
775 OFL^:. 

S1A=0, 

S2A=C. 

S3AT-C. 

GAHB^GAMPm 

PB=PPC1I 

ca^QPfii 

RMOB^PHCPCll 

THB^Tppfil 

HB^MPdl 

XNUB^XHUPf I) 

Y0=yp(ii 
A1=F1 «PI 

A2^F2(NtSlA,S2A,S3ATI 
IFfJCHEH.EQ.il GO TO 7254 
A3=0. 

GO TO 7257 
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7?54 TPlslTPIll^T2l/2. 

0TCMPai = CTCHPtll/2. 

00 1552 J^l.NSP 
1552 0UPH|J)=CACMP« J,l>/2. 

A3=F3 nPl»OTCHPCll.TP (l),T2,Tf-Plll ,PHE2*DUHP*WPC1I •WNCHI) 
7257 CONTINUE 
0PT=-1. 

IFICK/2I®2.EQ*KI 0PT=l. 

A4=F4fBETA»-0PT*XHUPm •THPfll ,XHUN|NI •TPN(H| I 
A2=CA2«A3t*A4 

PSH-PPUMCOPT»CPhe2-TMPfin*A2*OELXI/Al 
ER3= tPSH*P2H)/P|N> 

IFIA8SCEP3I .LT.1.E-3IG0 TO 19 
IT1=IT1#1 

IFfITl.GT*15>G0 TO It J 
IT11=1T1141 

IFIIT1,EC.2> GO TO 1430 

IFfERl*EP3.LT.3,| GO TO 14 

IFf AflSIEGl-ER3l .LT.5.E-06I GO TC 1492 

IFIAeSIEFll •GT. AeS(ER3ll GO TO 143G 

IMICCC.EC.ll GO TO 103 

xxx=*i. 

ICCC=1 
IT11=1T11-1 
143C ER1=ER3 
PETl=eET 

BET=PET^,Ql*CITll-ll*eET*X<X 
GO TO 15 

1492 0ET2-(BET-0ET1 l♦^0. 

EP1=EP3 
9FTl = eET 
0ET=0ET^OET2 
GO TO 15 

1C3 HRITE(6«9i91) 

WRITE <6*2201 

220 FCRNAT<» ERROR IN SHOCK ANGLE IN hsMOCK*> 

STOP 

14 0UN=8ET1-ER1*<8ET*9£T1M<ER3-ER1I 
ER1=EP3 

PETlsPET 

8ET-0UM 

15 YN<m =Y(NI ♦.5*<TAN<8ETe<K) I ♦TAN IBE T M •DCLX 
I5=IS<K> 

TN .51 =\N<N> 

KS=IS<KI 

LS=KS 

KT»KS 

IF(K*E0.2.0R*K*EQ*4I KT^KT-l 
CALL CPCINT. 
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GO TO <♦ 

19 epTftNTK) =8e- 

YNIMI =.P»(TANt8PTBTKI I ♦ T AN ( 9P T AM Kl n•OELX♦y^Ml 
I5=ISIKI 
YNIISI =YNfNI 
3535 CONTIKUe 

PNIHI =«>3H 
QNINI =C2 
THNTMI =PME2 
PHONTPI =«H2 
tn<ni=T2 

f»N<MI=PC/kNtf*l 
CPXKTPI -3. 

00 lUCl J=l.NSP 

CPN (J,P»=CP1 (Jl 

lAOl CPXN«PI=CPKNJP| ♦ALPNfJ.NI •CPMT J.N» 

GAPNtPI =CFKN|M|/<CPXN«M|-9N»N| /CP INI 
OP»»=l./PN<Nl 

<=^t'NIMI = OA THI •EPINF^SCPT TGAR/CAPMPI»ORP/TN<HI I 
IPTENXTPI .LT.l.OOCH GO TO 14C2 
xmiNTPI -2PlMEPt<PM 
XAC2 CONTINUE 
RETURN 
END 

SueROCIINE SMITCHU.KI 
CONNON/eA/AtP(T,55l .EHINF.HINF 
C0PNON/eC/X**ASS <551 
CnNN0N/CJ/CPT7,55l,CPlC7l .CPXC55I 
CCMMON/FF/EH55 l,GAN<55l,PI55l , TP <55 1 ,Y<55I 
C0PNCN/rC/N(55l .XI55I 
COPNCN/PO/JCMEn.NSP.T 1551 
COHMON/OA/rt<7,55l ,Q<55I .RH0C55I ,X«UI55I 
C0P*!0N/RC/R(55I 

CrPH0N/U/9Q(55l .CALP< 7,551 ,080 <551 , OCPX <55 1 ,00 ALP I 7,55 1 ,0TAU<55I , 
1TAU<55I 


t 

<JI = X 

<KI 

y 

<JI=Y 

(K) 

0 

<J)=0 

iK} 

p 

<JI=P 

(Kl 

T 

<JI=I 

(Kl 

W 

<JI>H 

IK) 

R 

(JMK 

IK) 

EM 

<JI=EH 

IK) 

TH 

<JI=TM 

(K) 

ao 

<ji=eo 

(Kl 

TAU 

<J»=TAU 

(Kl 

oeo 

<ji=Deo 

(Kl 

gam 

<JI=GAN 

(Kl 
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9HC 

CJMPHO IK) 


XHU 

fJIsXHU IK) 


cpx 

CJ)=CPX IK) 


ncpx 

IJI=CCPX IK) 


QTAU 

CJ)=DTAU CK) 


XMftSS (JI = KHftSS <K| 


OC 136 JJ=1.NC“ 


H 


(JJ, K) 

CP 

(JJ.JI =CP 

(JJ.KI 

ALP 


IJJtKI 


<JJ*J|rOAtP CJJ*K| 

COftLP(JJ*J) = OrALP CJJ*KI 

ira continue 

RETURN 

END 

SUBPOtTI^»' P*-C^»L.IPftN*K,OPT*KCPT ) 

CC*«f1CN/aL 

CQf*f^OK/et «LPf7,5E> 

COKNON/CJ ^CP< 7.GE> iCP; (7> *CP^t5Sl 
CCH'lON/CK/HTMOLt C7) 

COPMOK/FC/CPIf *PC 

COPMON/EF / ENf55KGA*^(55l,PC55l*rH(55l»y(5El 
CCrriON/«^r /EIn,p9^ 

CCHNCN/^" -/GANINF*H1 <7 > •PINF 
COP»<0N/ ' :/MC5El *X (35> 

COMflON/PC/JCMEP*NSP*T I5E ) 

C0HH0N/CA/HC7 ,551 ,0 CE 5 » •PMO C 55 I ,X*^Uf55> 

C0KH0K/ffC/9(55l 

IlrL^KCPT 

OP=A ,OG(F CM) /F (I It ) /FLOAT I IF AN- 1) 

IFF= [FAN 

2 J=UN$P 

2 h6=H6^»M J,H) *ALP ( J » 

00 1 LL=?*IFF 
N=tL^l*H 

IF(r'/3)»2,EQ.<) N-M-LL»1 
KK=N-1 

IFC(K/2I*2*EQ.x) K<=N^1 
X CN) sX (KK) 

PCNIsPCKlO/EXPCCP) 

ALNP-rP/CANIKK) 

PH0(N>=PHCtKK)/EXP|ALN9) 

G = 2.*GAMtKO/CGAM(KK) -1.) 

CQ=-G*CPCM/9H0 CNI-PC KK>/oHOCKK) I 

QCN)=CCKK|«QCKK)»QQ 
0CN)=SQPT COIN)) 

H2rH6-0G/2.*ElN 

IITlrl 


F •T' • 

» 
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T1=T(KKI 

IFfIITl.EC.il T2=T1*.99 
8200 CALL THE' 0IT2,H1,CP1I 
H2P* 0 . 

00 4 J=1*NSP 
ALPf J.Nt^ALPt J.KKI 

4 H2PsH2P»aLP(J«N)«Hl(JI 
ERP1= IH2-H2PI/H6 

IFfAeSfERRlI .LT.l.E-081 GO TO 8202 
IIT1=HT1*1 

IFIIIT1.GT.15I GO TO 82C3 
IF<IIT1.GT.2» GO TO 8204 
ERR2=ERR1 
T22=T2 
T2=T2*.99 
GO TO 82C0 
8203 MPITE C6. 91911 
9191 FORHATtlHl) 

MRITt. 82051 

8205 FORMATS ERROR IK TEMPERATURE LCCP IN PM*I 
STOP 

82C4 nUH=T22-ERR?* I T2-T22 » / « ERR1-ERR2 I 
t«R2=FRRl 
T22=T2 
T2=0UH 
GO TO 8209 
8292 CONTINUE 
TIN)=T2 
MINI =0. 

CPX(N|-C. 

00 5 J=1.NSP 
CP(J*M=CP1(J» 

H(J,NI=H1 Ul 

CPX(N)=CPX(NI »ALP(J.NI*CP1(JI 

5 MINI=H(N)»ALP(J*NI/MTNOLE<JI 
MtNIsl./MfNI 
R(NI=RO/MfNI 

GAMfMI=CPXINI/TCPX»NI-RfN»/CPIt 

ORN=l./R(M 

EMfNI*QIM*EMINF*SaRTfGAP 'GAMIN* 'ORN/T«N» I 
XHU(NI=2 mu(EM(NII 

THINI=THfKXI-OPT*ALNR*<COSUMU:K‘' , • SIN ( XMU < KK» I fCOS I XM'J < Nl »*SI 
INIXM><INII )«.5 
M6*H2 

1 CONTINUE 
RETURN 
ENO 

SUBROUTINE CHEHIFASI 
COHHON/eA/ALP(2,55l .fMlNF.MINF 
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COWMON/CA/wnOTNI7^55) tXNf55) 

C0t^H0N/Ce/8ETBI<*» 

C0HH0N/0P/YN(55» 

C0WH0N/£F/EH<55I *GAHI55I fPCSSl • TH f 55> *YC55I 
COMMCN/HI/DftLCHf . ) •nTCMEM 
COHNON/CP/ALPem *PMH55I 
C0MKCN/PC/W<5S> .X 1551 
Cn>f^ON/PQ/jCHrM,NSP.T <55) 

COHMCN/QA/H(7«5 5) • Q 15 5 1 . PHO < 55 > t XMU <55 > 
COMHON/QS/PHOP<e> •MOOT! 7,55 > tWOOTC<7> *WPC2I •XMUPI 
CCMMOK/VT/OACH( 7,551* OTCM< 551 •OYISO.V ISO 
COr'MOK/KV/MPTS^PE ,X5P,XJ 

COPHC^/Y^/8PReSU•CM£MFC♦CPPESU•EHSU0,R^H,XSlrP 

TPCJCHEM.EQ.T) GO TO ^351 

C CHCniSTPY PACKAGE ••••••♦»•♦ 

DC 5355 Lri,NPTS 

00 59 P=l,4 

IFdSIf-I.EQ.CI GO TO 59 
ITEST=ISI^>-1 

IM<P/2>*e.5Q.H) IT5$T=IS<?1) 
I£<L*EC.lT£ST.OR.L*EQ.ITFST^n 00 TO 5398 
89 CONTINUE 

FAT= AflSfPPKLM 

I£<(FflT.LT*.0:KCR*(FAT.GT.lOC.n GO TO 8398 

► =L 

OX=SQPT < <XN<L»-XfK) > »*2^<YN<U-Y(KM**2>*PTM 
DO 8350 J=1.NSP 
835C 4tPBCJI=AtP<J,K) 

CAtt PCCtS <T <K) ,P<K) ♦ Q(K) ♦PHOIK) • ALPB,nX,tl 
OC 5301 J=lfNSP 

83C1 oacm I j,L >= oalcm: J) 

OTCHILI-CTCHEM 
GO TO 8355 
5398 DTCH<L)-0. 

00 8399 J=1,NSP 
WOOT C J,U^O • 

MOCTN<J,U=0. 

8399 OACHCJtLMO. 

5355 CONtINUE 

GC TO 40CC 

8351 OC bice L=ltNPTS 

OTChCLMC. 

00 83C2 J=lfNSP 
HCOTCJ,L>-0. 

WCOTK<J,U>C* 

83C2 OACH<J,U = 0. 

61CC CONTINUE 
•♦oro CONTIMiF 
RETURN 
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ENO 

SUBROUTIKf, SHFA»<I. ASHE**}! 

C0KH0K/CJ/CPt?.55l ,CPH»» .CPX«55> 

CCMH0KXCK/MTM0LE<7» 

C0PK0K/£P/EP(55I.GA>M55»,PI5EI , TH 1 55 » . Y « 55 J 

C0PM0N/&K/3ELX 

CCPNON/PC/HC55) ,X (55» 

CPHHOF/FG/JCMf P.K5P.T 1551 
COPlOF/aA/Ht?,-^ ,a«5 5l .RM0<55I .<*'U<55> 

COPHOt/SS/ALl.ALZ.aOl ,0C2tCX.C2.CMl.CH2,OSl .DP2.031 tO02.0Tl.0Ti ,0W 
Al.0W?,PKl.P)(2.TAl,Ta2,Trtl,TH2,Wl.W2tYl,Y2 
COHMON/ST/113. IKE CI.K ,<FIPST.K<KQ,PSTAP 

C0P«0F/Tt/5OI55» ,0ALP<7,55» ,060 1 55) ♦ 0C«>* C55 » , 00 ALPI 7 ,55 1 ,DTAU»55I . 
1TAU155I 

COHMCN/VT/OaCHI 7,55) .0TCHI55) .DWlSD.tflSO 
CCHHCK/WV/VPTS.PE »X0P, *J 
DIMENSION S3Dt7) 

Vl=V2=tf ISC 

CV1=0V2*CVI5C 

TA1=T«2=TAU(K> 

0T1 = 0T2 = 0IAU «K) 

Yl=r2=Y<kl 

TH1=Tf2=TH(K) 

S10=SHXJ,RE) 

CM?0=C. 

00 1C J=1.NSP 

10 rH20=CM2C*0ALP( J,k»*CP(J,K» 

BGl = eC2=ECU» 

Cl=C2=CPXJKI 

08i=oe2=ceo(K» 

PXlsPX2=CCP» «K» 

CH1=CM2=CH20 

S20=S2<XJ,«E» 

S30T=C. 

00 23 J*1.NSP 
AL1=AL2=CALP(J,x) 

OD1=OC2sCDALPIJ.K» 

S30« Jt=S3«XJ,Pt» 

2C S3CT=£3CT«S3C»J»/WTMDL''IJ) 

PK=l,/PIXl 

SM1=S10/GAM<K) *PK/EM<K» »*2 
OK=l./OCX» 

SH2=- IGAMIKI-1. ••S2D/GA-tiO •PK'OK 
$H3 = “MKJ»S3 DT/RmOUI ’OK 
IE(XJ.EQ,C» SH4=0. 

IFIXJ.EC.O GO TO uO 

IF«K .NE .l.OR. Y«K> ,GT,l.E-6» GO TO SC 

SH4=TH(2>/V<2» 


ORIGlNAr PAGE IS 
OF POOR ^UALTIY 
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CO TO 40 

30 SH4=SlN«TH<KI)/r<KJ 
40 CONTINUE 
00=1. /OEtX 

SH5=-0TCMK»/T«KI *00*005 (TH(K I I 
0UH=0. 

00 5P J=1.NSP 

50 OOM=OUN*CACh<J,KI/MTHOL£<J» 

SH6=-m«J *OUM/OELX*COS(TH«Kn 

SH=SH1*SH2*5M3*SH4*SH5*SH6 

ASHE AP=-GAH<K»*P(K)*E NIKI **?*SH 

RETUPN 

ENC 

SUeROUTINE PRFSS<K.P. TH.TMNI 
COMMON/ AC/ IB 00. PIN 

common/hx/apress.aoresu 

COMMON/ VX/A BOOS .ePRESS.CPRESS 

P=APR£SS»X*(ePRESS*CPP5SS*XI 

P=P*PIN 

THN=TK 

RETURN 

PNO 

FUNCTION XVISIAl 

C0MM0N/BA/ALP<7.55I.EMINF,WINF 
COMMON/Oe/BE T0<4» ,IS<4» 

C0MM0N/FH/XK1.XK3.XP0T 

C0MM0N/QA/H</,55» ,0 «55» . RNK 55 I . XMb« 55 I 

COPMON/MW/NPTS.fiE.XBP.XJ 

data IVI5 :•/ 

I0UM = F' IS 
100=1 

TFI15l4l.tT.OIIOO=IS(4Ml 
IFIISOI.GT.OI IOUM=ISIi»-l 
OUM1=0. 

RU=RHOaCCf*O?IO0) 

IOE=IOO*1 
00 10 I=I0E.I0UM 
OUM=AeS<fiHOni*G«ll-RUI 
IFTDUM.lt. OUMII GO TO 10 
0UH1=0UM 
10 CONTINUE 

IF(ALPI4»1».GT..99 » GO TO 30 

IFIALFT4.H .LT. 1.E-O9.ANO.X0P.LT.XPOTI GO TO 30 
IF! IV IS. EC. 0> XVIS1=XK1*RF*XBP*0UM1*XK3 
IFIIVIS.EC. 01 
1 MPITET 6 . 96961 XBP.xVISl 

9696 FORMATI32M VISCOSITf MODEL SMITCHFO AT x =E1C.3.14H 

tF13.5» 

IVIS=1 


VISCOSITY 
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XVTS-XVISl 

GO TO 

30 CCMTINUE 

XVIS^XKl ♦Rr*xep*OUMl^ XK3 
4C CONTINUE 
i^ETUPN 
ENP 

SU«‘?OOTINE COWL 

CrWMCN/eA/ALP(7,55l *EMlNE,WlNP 
CCM»10N/CJ/CP<7,55) tCPK 7) ♦OPXiS^) 

C0HM0N/0e/B£TR(4) » IG< 4 ) 

CC^HC^/OF/»^M 

CCMM0N/EF/E'^(55 >tGA»1( 551 ,P(5SI , THC5SI ,r(55 ) 

CCHHC^/HJ/KOUNT *LZtMPT 
COPMON/HL/ALPHA,0ETA 

C0PHCN/Hf/ALPM'',55) *CPN(7,55l * JPXNC55I »EHN < 55 ) «G ANN ( 5 5 ) tHNC? *55) * 
ILS*PN(55),QNC55) *PH0Nf >5) ,RN(55 ) • TmN ( 55 • • TN ( 55 ) * WN ( 55 ) ♦ XMUN ( 5 5) 

cohmon/hp/petan ( 4) • iEwe:o 

rOH«ON/PC/W<55) *x<5t>) 

CCPMON/PC/JChFm,nSP*T (55) 

COP»«ON/GA/HC7,55) *.)(55) *RH0(55) *XHU(55) 

COPPON/PC/R(55) 

C0HM0N/TL/8Q(55) tCALP (7*5 5) *000 ( 55 ) . OCpx (55 ) *00 ALP I 7, 55 ) *DTAU(55) * 
1TAU(55) 

C0PWCN/WV/NPTS*PE*X1P .XJ 
ALPH A=l. 

PETA=C. 

NPTSSS-NFT 

NPT=NPT-3 

IE(P(^FT♦HM)-P(NPT)) 3CC1*ZCCZ*20j3 
20C2 WRITE (6.9191) 

9191 FCRMAT(lHl) 

WRITE (6.2004) 

20CA FCRflAT(lHl,67H ERPO^ IN INPUT DATA - NO PRESSURE DIFFERENCE ACROSS 
1 SPLITTER PLATE/23H SET INPUT - INTACT = 0) 

STOP 

20C3 0PT=l. 

K = 4 

IS (K ) -NPTSSS-1 
L=IS(K) 

M=L^HM 
GO TO 2005 
20C1 CPT=-1. 

K=3 

IS(K) =NPT4MH 
M=NPT 

L=Mf WK 

2005 IFAN=NM*2 

K0PT=CPT»1,5 
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N13 = NPTSSS-2»K'>‘>T 
IF<IK/2)*2.EQ.K» M^=NPTSSS 
00 3PC M2 = N1J,MPTS 
K5 = NP1S*M J-N12 
J5=«<5*1 

CALL SWITCH CJ'^ .KSI 
300 CONTINUE 

NPTS=NPTS*1 
DO 331 111 = 1.*. 

IFUll.EC.KI r.O TC »C1 
IFTISTIlll.&T.NJTI r..<Ill»=IS(IlH* 
3*:i CONTINUE 
LL = L«-KCPT 
TftL'«tL»=C. 
pn<LI=PTLI 

IMNIU |r Ih (LI 
Q NTLMC <LI 
T N(LI=T (LI 
P N(LI=R (LI 
PMOMLI=PFO(LI 
CAHN (LI =GflM(LI 
PQ(LLI=-0. 

DCPX (LLI=C. 

0TAU(LH = a. 

08Q(LLI=0. 

CPX( LLI = C . 

WN(U=W(LI 
00 2106 J=1,NSP 
HN<J,L»=H(J,LI 
0ALP( J.LLI=0. 

OCALP(J,LLI=0. 

20C6 ALPN(J,i.*-ALP(J.LI 
ITT = 1 

PtTBIf ,-(TH(LI-OPT*XmuU: I*!.!! 
BET=Beie(Ki 

2c:r iEMaeo=i 

CALL HSHCC<(KI 
IEM3EC=C 
KK=LL*KCFT 
P TLLI=P NILLI 
0 (LL)=C N(LLI 
T (LLI=T N(LLI 
W «LLl=N M(LLI 
P (LL}=fi NILLI 
TH ILLIsTH NILLI 
EH ILLI=EH NILLI 
XMUILll=XPUNlLLI 

cahilli-gahnilli 

»HOILLI=RHONILLI 
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CPX<LL1=CPXN(LLI 
P (KKI=P N(LLI 
Q (KK|=C NtLLI 
T (KKIsT N(U) 

H (KK|sH N(LU 
9 (KKMR NILL> 

TH IXK»=TH NILLI 
6H <KKI=EP N(LLI 
XMUIKKI = XhUK«(Ll.l 
G«P(KK)=G/>MN(I.LI 
RHOfKK)=RHON(LLI 
CPX<KK»=CPKN<l,UI 
00 20C8 J=1.NSP 
H NU.LLI 

CP <J.LL»=CP N(J,Lu» 

ACPI J.LU=ALPN( J.LLI 
H (J«KK)=H N(J.IL> 

CP IJ,KKI=CP N(J*CLI 
2606 ACPI J«.^K|=ALPN(J«CCI 
X(LL)=X<L) 

X(KKI=XILt 

tms=thkki 

CALL PH(P.L,IPAN,X,CPT,KOPTI 
NNNsIFAK-1»M 

IFHK/2I*2.EQ.ki NNN=H-IFflN*l 

THPM= tMINNNI 

FRR=THS-THPM 

xFCAeSIEKRi .LT.1.E>0<.| GO ; , 15 
ITT=ITT»1 

1FI1TT.CT.15I GO TO 1C2 

1M1TT.GT.2I GC TO 1<4 

EP1*EPR 

BET1=8ET 

OET=1.01*BET 

eCTtIXMBET 

CO TO 290? 

1Q9 MR1TEI6.203> 

2C3 FORHATI* E9P0P l^ BETA SHOCK IN COWL*l 
CALL EXIT 

14 OUHl«BETl-eRl»l8ET-0ETlI/(ER‘' -Efill 
EI»l*eRR 

8ET1»EET 

aET«'JUMl 

OETBIKIi^eET 
GO TO 2007 

15 CONTINUE 
NPT*NPTSSS 
RETURN 
END 


Oi^IGlNAL PAGE IS 
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SUBROUTINE OPOINTtK-* ! 

COHH<N/A /GAR. GEN 

cocrok/e... alptt.sgi .eninf.hinf 

C0PN0K/CJ/Cr(7.SSI*CPll7I.CPX(55t 
COMHON/Ci' 'HTNOLE 171 
COHNON/DP/VNISSI 

CCMN0N/EF/E*M5St«GA!t(9SI,P(S5l*THI5S| .VIS5I 

C0PHCN/EF/6AHINF.H1(7I,RINF 

COPNQK/GK/OELX 

COP'OK/PC/MtS^I.XtSS) 

COPPOK/PG/JCHFP.NSP.T(SSI 
C0PH0K/CA/HI7«S9> .0(591 •PHOISSt .XNUfSE) 

COPHOK/TU/BO<S5I.OALP(7.95I.OPQ(S5I.OCPX(95I.OOALP(7,55I«OTAU(S5I. 

iTAUfSSI 

C0PKCN/TV/ALPP(7.2>*nET.RaP(2l.0ACHP(7.2l,0ALPP(7.2l ,0BQP(2). 

1' 'PXPI2I .C0ALPP(7.2).0TAUP(2I*CTCHP(2)*GAPP(2I.PP(2I.QP(2I. 

2. AUP(2I.T»-P(2».TP(2I.YP(2I 
COPHON/VT/OACH(7.9SI.DTCh(95I .OVISO.VISO 
IT=1 

»o=(yp(ii*rp( 2 ii/ 2 . 

16 RAIs|T 0 -»F(in/(TP( 2 l-TP(lll 

AlAM 0 =TAK(THP(ill«RAr*(TAN(THP( 2 n-TAN(THP(in I 

yAT=y6(tl-ALAM0*0eLX 

ERR=ABSi(yAT-Y0i/iyp(2i-yp(ini 

1FIERP.lt. l.E-051 GO TO 19 

yo=yAT 

IT=1T*1 

IF(IT.LE.lG) GO TO 16 
HRITE(e*9191» 

9191 FORMATIIHII 
NRITE(6.2C2) 

202 F0RBAT(« ERROR IK 0 POIMT ITERATION*! 

STOP 


VIKI 

• YO 




P 

iKIsP 

P<l>*RAf*(P 

PI2I-P 

Pllll 

Q 

IKIaC 

P(1I*RAT*(0 

Pt2l*Q 

Pin) 

T 

CKI«T 

Pm*»AT* IT 

PI2I-I 

Pin) 

TM 

IKIsTM 

P(1)»RAT*(TH 

PI2>-TM 

Pill) 

BQ 

CK)^ec 

P(1I»RAT*|B0 

Pt2i*eo 

pini 

TAU 

IKI^TAU 

PI1I»RAT*(TAU 

PC2I-TAU 

Pill) 

OBO 

CKI«0eQ 

P(1I*RAT*(09Q 

Pf2l-O0O 

pini 


OCPX (KlaOCPXPd) »RAT* OCPXP(2l -OCPXPI 1 1 1 
0TAU(K»*CTAUP(H*RAT*(0TAUP(2»-9TAUP(in 
OTCH(KI*OTCHP(ll»RAT* (0TCHP(2l>0TCHP(ltl 
CPX(KI*e. 

H(K)sC. 

CALL THERNOITIKI.Hl.CPlI 
DO 1 J«ltKSP 
H (J*Kt«H KJI 


C ^ 2- 
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CP(J«Kl=CPltJI 

ALP tJ.KIsALP P(J«t)*PAT»(ALP P(J«2)-ALP PU.lll 
□Al tJ«KI-OALP P*J,H*»AT*<0ALP P*J»2I-0ALP PfJ*in 

CACHIJ,KI=OACHP(J,ll»PAT«(OACHPIJ«2l-OACHPf J,ll) 
OOALP<J«K)«OOALPPIJ.ll*RAT*ICOALPP(J*2)-COALPP(J*tn 
CPXItCMCPX(KI«ALP(J*KI«CP|J.KI 
1 MIK>=WfK)»ALPU*Kt/UTHOLe(J) 

M<KI =l./MfKI 

«?MC<lO*P«KI*mKI*G?M/TCK» 

RETURN 

END 

SLPROLTIAE STEPIVISI 
COHNOK/AC/IPOO.FIN 
C0KN0N/AX/JSUBL«JSU3U 
COMNOK/CA/HOOTN (7.951 *XN(55i 
C0PNCN/C8/BETe(<*l .15(1.1 
C0PHCN/CP/YN(55l 

CCPN0N/£F/E*«(55 I.GAN(55t .P(55l . TH(5SI .V(5SI 

CCPHON/EG/EIN.PP.KLE 

CCNNON/GE/RAO.RCO.UIN.YISINF 

OOPHON/GX/OELX 

COHHON/HJ/KOUNT.LL.NPT 

COtiriON/HI'/ALPM(7.55l.CPN(7.55l .CPXN(55I .EHN(55I .CANNISSI .HN(7 .551 . 
ILS.PN(55).QN(55I .RH0N(55 I .RN(55I .THN(55I .TN(55I .riN(55l .XNUN(55I 
COHHOK/CR/ThBP. Y 5P. YBPN 
COPMON/FO/JCMEM.NSP.T (551 
COPNON/GA/H( 7.551 .0(551 .PU0(55I .XNU(55I 
C0PP0N/ST/U3.IRECI.KS.KFIRST.KKKQ.PSTAR 
CCPNON/UV/IIll. lEPR. IP9ESS. IPPESU. ISUfl 
COPPON/V ICONT.ICNO.KT.THBPN.X-JPN 
COPHON/M\i'NPTS.P€.XnP.XJ 
C0PM0N/MX/APRFSS.APRE5U 

CC««ON/Y7/3PRESU.CH5MFC.CP9ESO.eMSue.^TH.XSTEP 
nIP5NSIC^ OELLK (55I.KM(2».YH(2I .THH(2l 
DATA I11/C7 
CATA IPEGl/l/ 

ISPOaC 

ISPA=0 

NSAVEsZ 

KM(1I«C. 

XH(2|slOOOO. 

YMdtslOSGO. 

YM(2lxlQGC0. 

THHdUl. 

THM(2)sQ. 

ISUB»0 

JSUBL*NPTS»1 
JSUPU*NPTS*1 
OC 910 Isl.NPTS 
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IFfEHtll.GT.EMSUai GO TO OtC 
JSUOLsl 

isue«i 

GO TO aoe 

910 CONTIItUE 
GO TO 601 
600 CONTI^UE 

00 602 Isl.NPTS 
II*NPTS-m 

IFtEHIIll.GT.EHSUei GO TO 6C2 

jsueu*ii«i 

GO TO 631 
602 CONTINUE 
aCl CONTINUE 

IFdSUe.EQ.QI GO TO 10 
113=1 

KFIRSTsKCUNT 
10 CONTINUE 
jouN=jsuec-i 
NP2=NPTS-1 
00 499 K=1«NP2 
0ET=Y«K*1I-YIKI 
IFOEY.lt. 1.E-06T GO TO 496 
IFfK.GE.JSUOL.ANO.K.LE.JOUMI GO TO 49S 
EHlsKNl<1..0..TH|KI,XNUfKI.O..O.) 

EH2=XH2fl..0..TH(K*ll .KNU(K»1I « J.«0.t 
OELLXtKI=(Y(K»ll-YIKI I/(CM1-EN2> 

IFI8ETefll.Eq.O..ANO.eET6(2l .EQ .0 ..ANO.OETBTSI .EQ.3 .. ANO.SETOf 41 . 
AEQ.O.IGO TO 499 

IF(eETeai.GT.O..OR.OETP(3i .CT.I.ICO TO 499 
IF<K.EQ.IS(lt.OR.K.EQ.IS(3MOELLXIK;=DELLX(K-ll 
IF«K.EQ.IS(2I.0P.K.EQ.IS(4II0ELLX(K-1I=0ELLXIK>2I 
GO TO 499 

496 OELLXfKI=l.E*06 
499 CONTINUE 

OELXNsCElLXm 
00 501 K=2.NP2 

IFOEtLXIKI .LT.OELXNI OELXH=DELLX(KI 
501 CONTINUE 

IF<aETem.EQ.O..ANO.BeTe<2t.EQ.O..ANC.eETB(3I.EQ.O..AN0.8ET8(4). 
AEQ.O.IGO TO 777 

IF(8ETB<1T.GT.0 ..0R.BETB(3) .CT.C.I 777.776 
777 OCHAR=OELXN 
GO TO 775 

776 DCHAR«2.«0ELXN 

775 IFfXLE.EQ.O..OR.VIS.EQ.O.I60 TO SO 
VI=1./VIS 

0EtW=.5»PR*RE*WI/XLE 
60 TO 51 


IOKIGINAS page is 
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50 OELll*3. 

51 00 5C2 K»1.MPTS 
OCLVl*l.C»i: 

IFtK.NE.ll OELV^yIKI-VIK-II 
IFtK.FE.NPTS) DELVYsV«*l)-YtK) 

IFIK.EC.I) OELY-DELYV 
IFK.EQ.NPTSI CELYYsOELV 
IFIOELYY.lt.OELYI OELVsOELYY 
IFIOELY.lt. l.E'QOl GO TO 502 
OELV1 sOFLV»PHO(KI*OIK 1«OELY*«2*COS(TH(KII 
502 OELLXIKI«CELVt 
OELYHzOELLXUI 
00 50<I Ks2*NPTS 

5C4 IFfOELLX(K>.LT.OELYHIO£LXHsOELL<(K> 

OSHEAFeCELXH 

OELXsl./(l./OCHAP»l./OSHEAP) 

OELXsOELX/XSTFP 
IFTISPA.EC.tl GO TO i» 

ISPAal 

call SPACEIISPP> 

IFIIIll.EQ.il PETU<tN 
IFdSPP.EC.ll GO TO 10 
4 COKTIKUF 

IFUCHEN.EQ.OI GO TO 4275 
00 505 lal.NPTS 
0TEST*CIII*UIN*4.E-7/RTM 
0TE5T*CMEPFC*0TEST 
505 IFTOELX.GT.OTESTI OELXaOTEST 
4275 CONTINUE 

IFII13.KE.1.0R.K0LNT.NE.KFIRSTI GO TO 4545 
KKKCaKCU4T»2C 
4545 continue 

IF lENTJSUeui.LT. 1.051 KKKQsKOUNrtl 

xHT*xePtc'Lx 

RAsl./RAC 

IFIXHT.LT.XHI1) I GO TO 741 
IFIXHT.LE.XH(NSAVE)) GO TO 5209 
THOPNaO. 

XaPNaXH(NSAYE) 

VBPNaYH(NSAVE) 

II11>1 
GO TO 5210 

52C9 CALL TBLIXHT.THBPN.XH.THH.NSAVE) 

GO TO 5204 

741 TMBPNaTFeP*OELX/COSITHBPl*RA 
5204 OELXaCELX/COStTNBP) •(COS(THBP) »COS(THePNI I /2. 
XPFNaXMT 

YOPNaYEP* CSINITFBPI »SIN(THBPN) I •.5*DELX/COS«TH8P) 
5210 CONTINUE 
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00 5211 I*1»NPTS 
XNCIIsXePK 

YN<II=VIII«TAN(rH(III«OELX 
5211 THNm>Th«II 

lF(inC0*EC.ll CALL ROOVtXNdl.VNIll.THNfil.O) 
tFIIPRESS.EQ.il call PQESSfXNIll .PNdt.THlll.THNf II I 
IF (IPRESU.EQ.il PN(RPISI=PIK*(APRESU*XN(NPTSI*<BP«ESU*CPRcSU* 
IXNINPTSIII 

IF (IPRESt.EQ.il HETURK 

IF(Y(NPTSI .EQ.YBPI GO TO 6211 

CALL eOOV(XMr.VN(NPTSI.THN(NPTSI«ll 

XN(NPTSI-KHT 

RE TURK 

6211 XN(NPTSMX8PN 
VN(MPTSI=YflPN 
THI«(NPTSI=THePN 
RETURN 
ENO 

SUBROUTINE SSONICdOGI 
COPMON/AC/lROD.FIN 
CCPHON/AX/JSUBL.JSURU 
C0HN0K/eC/XHASS(55l 

COPPON/CC/AUP.BUP«CUP.OTSPRI(55l,OUP,EUP,JCONWtTHPRI(55l«YPRI(55l 

COPHON/OP/YN(55I 

CCRN0K/EF/€R(55I.CAMI55I,P(55I«TH(55I.Y(S5I 

COHROR/GK/OELX 

COPHONi/PJ/KOUNT.LL.RPT 

COnnON/HL/ALPHA.SFTA 

C0Nn0K/HP/ALPIi(7«55l«CPN(T«55l*0PXN(55l«ENN(55l*GANN(55l*HN(7.55l * 
ILtPN(55l«QN(55I.RH0N(55l,RN(55l *THN(55I *TN(55I .MNIGGI .XNUN(55 I 
CCHN0N/RC/AP9 • API ,AP2 

CCPN0N/ST/I13.IRFCI.K .KFIRST.KKKO.PSTAP 

COPRON/UV/IIll. lERR* IPRESS. IPRE3U. ISU8 

COPNON/VM/ICONT.IENOtKT.THRPN.xePN 

COMNOK/MV/NPTS«RE*KRP*XJ 

OIPENSICN E2(2I.£RET(2I,FRTHL(2I 

DATA LKIF/1/.LNE/1/.LICQ/0/.LICK/0/.L1JUHP/9/ 

DATA OTHC/0./ 

IF(I13.NE.1.0R.K0UNT.NE.KFIRSTI GC to 1777 
XKFcXeP 

call 0PCTH(THS80T*JSunLI 

ATH8sV(jSU0LI 

BTHBs T AN ( TH USU8L 1 1 

1F(IRE6I.EQ.OICTHB«THS5CT/COSITH(JSU8LII**3 
IF(IREGI.NE.0ICTH8>(CTHRn 0TH8*X0ELI 
OTHBoOTHO 
XMKFsXMASS(JSUeui 
1777 CONTINUE 
Ksjsueu 
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OQsXePUHkKF 

0Ssl./6. 

00*1 

YNIKIs<IUP48OP«C04CUI>»0Q*»2»#5^0UP«0Q»*3«OS^CUP*DQ*»ii*0Q 

THIllKI*ATANCeUP4CUP*OQ«OUP*OQ««?*«S^£UP«OQ»«3»OSI 

0S*2*«0€LK/fC0S<ThllO l♦COSITHNIKII I 

THGN*THKCIO 

VGH*VNCKI 

ieooi*iecc 

ieoD*i 

IPRES1*IPR€SS 

IPP€SS*0 

IPRE Ui*IPRESU 

IPRE$U*0 

AtSV* ALPHA 

8ESY*e£TA 

CALL LPOINTCJSUeUtO*) 

K*JSU8U 

THNIKI*THGH 

YN(KI*YGH 

ALPHAs«S 

PETA«.5 

CALL LPCILTiJSUeUtO.I 

alpha^alsv 

eETAsRESV 

K*jsueu 

THHCKI’ThGH 

YNIKI*YGH 

1800*18001 

IPRESS*1PRES1 

IPRESU*IPREU1 

DS*a • *OEL X/ CCOS C TH ( K M ♦C OS f THN ( K I M 
PSA* (PNIKI*PCKI I/OS 

THSA*fCUP^OOP*CQ^€UP*OQ»»2»*5l*COSfTMN<Kn*»3 

PNA«-CAHMKI*PNCKI*EHNIKI**2*TMSA 

PYU*CCSnHNIKII«PNA«SlN(THNfKI »*PSA 

K*JSU8L 

00*X8PN*XXf 

YNlKI»ATH0^BTH8*COACTH8»OO**2*.SfOTMB*OQ*^?*OS 

fMNIKI*ATAN<BTMe^CTHB»OO^OTM8*OQ»»?^#5l 

TMS8* <CTH8tOtHe*OOI ♦COStTMHU) I ••S 

OS*2.«OELX/(COS(THIKI I ♦COSITHNIKM I 

PY8*AP1 

ITElal 

I£T*0 

1790 AP2* I PYt'-PYBI 5/ ( YN C JSUBUI -VN ( K II 

AP0«PNCjSUBUI-*9«CPVe^PYUI*CYH<JSUBUI*YNIKII 

AP1*PY8 

PNtKI*APO 
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ICONTsl 

ICNO^O 

KPsK 

L*K 

KT»K 

PSTARsPNIKI 

VGHsVPIKI 

TH6M=THNIKI 

CALL GPC1^T 

YN(KI«r6H 

TNNCKIsTHGH 

PNB>>CAHMKI •PN(KI«€NN(KI ••2*THS0 
PVeZsCCStTHNlKl l«PNB 
eTsPVE-PY82 

IF(A0S<ETI.LT.1.E-CSI GO TO 17A9 
IET=1ET*J 

IEfIEf.LT.20t GO TO 6532 
MRITEI6*€533t 

6533 ECRHATI* £T LOOP IN SSONIC»» 

STOP 

6532 C0KTI6UE 

ERETIITEIMET 

IFfITEl.GT.il GO TO 356 

1TE1*2 

PYBl =PYB 

PYP3PY82 

GO TO 1750 

356 PYeo=PYai-E»ETfll*CPY0-PYOll/ieRET(2l-ERETIin 
PY01«PY8 
ERETfllsERETfZI 
PYBsPYBO 
GO TO 1790 
1769 C0NTI6UE 
ICCNTsl 
lENOsC 

JSU8L1«JSU0L»1 

JSL'8UlsJSUSU-l 

DO 1736 KK<JSU6L1*JSU0U1 

K=JSU8U-KK*1 

THNfKIsTMKI 

YNIt<l*TA6(TMfKII*0ELK*Y«l 

PNfK ' =AP0*APl«fYNtKI-YN(JSU8LII *AP2 YNf Kl -YNf JSU8LM *«2 

KIPal 

NE*1 

KP«K 

L«K 

KT*K 

PSTARsPNCKI 

DS>2.*OFLX/.ICOSfTHfK> MCOSITHNI Kill 
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VN(Kl3V(KI».6«(TANITHIKn»TAN(THN(KII t«OELX 

VGHsVKtKI 

THGH= THNfKI 

CALL CPCINT 

VNlKl3VGt< 

THN(KI«THGH 

VN(K|sV(KI*.5*fTAN(TH(KI)*rAN(THN(Kll»«0ELX 

KP»K*1 

TERKs (RHCN(K) •ONI KI »C0S<THNUI I »RHON<KPI •QNlKPI*COS (THNiKRt II /2. 

XKOUMs)i|'ASS(KPI»TEPM«IVNtKI««(t.»XJI-VN(KPI«*Ct.«XJII/<l.*XJI 

QlJsl.^XJ 

VN<K|3|V|«|KPI*«C1J»Q1J*(XNASS(K>-XHASSCKPI l/TERMI ••U./QIJI 
173<» CONTIKUE 

IF(KOUNT.RE.KKKQ>ll RETURN 
XMOIEFsXPKF-XIIASStJSURUI 
JK«JSUeL*l 
OC 3V7 t=jK«jsueu 
3<>7 KHASSai = KNASSm«XNniFF 
jsuisjsueu'i 
MRITEf6«lAiai 

1A18 FORKAmCX«*CCRfiECTEO INTERNEOIATE STREAPLINES»/2X*«STREAHLINE NO. 
1*,7X,*X*«12X,*Y*,IIX,*TM*I 
DC 3«€ KR«1*JSU1 
KsuSURU^KK 
KIP«1 
HE«t 

OTERMsO. 

OVOX«TAN(fNPftI(Kll 

02VCK2>C1SPRI(KI/C0S(THPRI(KII««3«.9 
IFIK.E0.1I02VOX2» CTH0*.9 
X0ELsX8PN-KKf 

6030 ySTAR«YPRHKI*OYOX*XOEL*02YOX2*XOEL**2*Of£RH*XOEL**3*OS 
TH8rAR>ATAN(OYOX»2.*n2VOX2*XCEL*OTERH*XOEL**2*.9l 
KP»KR1 

TERN* (RHONfKPI*QN(KO| •COStTHNIKPI l♦RHONfKI«aN(KI •COS< THSTAP.I I /2. 
XtJ*l.«XJ 

xhoun*xmass<kpi*tern*iystar**xij-ynikpi**xiji/xij 
E2 INE l«XI*OUN-XNASSI KI 

lF(AeS<E2(NCII.LT.l*E«06 I 60 TO 603« 

KIP*KIP*1 

60 TO <60Al«6C4*2lfHE 
60 Rl NE«2 

OTERNl«OTERN 
OTERN«-.01/XOEL**2 
CO TO 603C 

CORE 0TERN0«0TCRNl-EZtll*{nTERN*0TERNllxfE7(2l-E2<lll 
OTERNl*OTERN 
OTERN*OTERND 
EZ(il*EZt2l, 
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IFIKIP.LE.2QI GO TO 603C 
MRTTEt6*60ai> 

60«1 FODHATt* TOO MANY ITERATIONS FO% ONE POINT IN SSONIC*) 

STOP 

6034 IFlK.EQ.t) GO TO 6036 
YNIKUVSTAR 
THNIKIsTHSTAR 
X0ELl>XCEL/4. 

DO 1417 1*1.4 

XOE*XOELl*FLOATIII 

XPRNT»XKF»XOE 

YPRNT«YPR1IKI»OYCX*XOE»02YOX2*XOE**2*CTERH»XOE**3*OS 
THPRNT*ATANI0Y0X*2.*02Y0X2*X0E4DTERH*X0E**2*.5I 
NRITE 16.1419) K.XPPNT.YPRNT.THPRNT 
1419 FOPNATISX.I9.SX.3E13.S) 

1417 CONTINUE 
GO TO 306 

6036 ERTHsTHSTAR-TMNIJSURU 
JCONV*0 
113*2 

ERTHL(LHE)*ERTH 
IFIlIJUNP.EQ.l) 60 TO 29C1 
lFIAf)S(EFTHL(LNE)I.LT..OCl) CO TO 2S01 
LKIPalKlP*! 

GO TO 12902. 26031. LHE 

2502 LHE«2 
0 TH 82 * 0 THe 
OTHQ*CTHe*.05 
GO TO 2604 

2503 IFaiCC.EQ.l) GO TO 2506 
IFfERTHL«l)*EPTHH2).ur.0.) GO TO 2505 
IF(LICK.EQ.l) GO TO 2606 

LICK*1 

RTHL*-.5 

IFIA8S(ERTHL(2) I.GT.AeStERTHLtini RTHL«-PTHL 
2506 IFUeSIERTHLten.GT.AeSIERTHLdll .AN0.LKIP.6E.4) GO TO 26J7 
DTHai*OTHF2 
OTHOEsOTHB 
OTMQ»OTHO*RTHL 
2509 ERTHX*ERTHLf 1) 

ERTNLfll*ERTHLI2l 
IF(LKIP.LE.IO) GO TO 2504 
WRITE 16.2508) 

2500 FORNATI* TOO HANV ITERATIONS IN LONER HALL LOOP IN SSONIC*) 

STOP 

25C7 PNB«ERTHL(l)*DTM8**2“ERTHLf2) *0TH82**2*ERTHL (2) •OTH0l**2-ERTMX* 
10TH8**2«ERTHX*0THB2**2-ERTML«1)*DTH81**2 
PHC*E9TMH2)*0TH82-ERTMLIl)*0THe*ERTHX*0TH8-ERTHL<2)*0TN81 
1*ERTMH1)*OTH81-ERTHK*OTMB2 
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0THQ*-PHe/«2.*PHCl 
LIJUNPal 
GO TO 26G4 
2S0S LICQ^l 

OTHO«OTHe2>ERTHLtll«(DTHR-OTHa2l/(ERTHL(2l*ERTHUlll 
OTHais|3THe2 
0THa2>0THP 
OTHOaCTHC 
GO TO 2SC9 
29C1 JCCNV-1 
U?aO 
LKIPal 
LMEsl 
LICQ-C 
L1CK=C 
LIJUHP*0 
IREGI»1 
OTHQ«C* 

25C(> IFUCCNV.EQ.CIOTHOsOTHO 
366 CONTINUE 
PETURN 
ENO 

SUOROUTINE POOV(Xt«Y«TH,IO) 

CONNON/ AC/IP OOt PIN 
COPNON/XO/KOO 

COHNON/2Y/APOO*aPOO*CaOO«E900«FPOC«GeOO.IAVEtIPUNCH,JBOOtKKKKK 

X»X1 

IFdO.EQ.ll CO TO 4 
IFlieOO.EC.OI GO TO 1 
XxXl-XOO 

Y>Aaoo»x« taaoo* »«caooi 

THa4TAN<eeOO*2.*CBOO*XI 
GO TO 2 

1 VxD. 

TM»0» 

GO TO 2 

4 VsEBOO«X«<Feoo»X«GBOO) 

THxATAN<FeOO«2.*GROO*XI 

2 RETURN 
ENC 

SUeROUTINE OPOTHIOTOS.II 
COPNON/AX/JSUBL * JSU3U 

C0NHnN/EF/ENI99t,GANfSSI*P<56)tTH(56) *V(SS) 

CONHON/RC/APO*AP1,A02 

CALL SNEARfl.ASPEARI 

ASHsASHEAR 

PV«APl*2.*AP2*<YIIi-Y USUBLI) 

02«VCI*1»-Y<I» 

ni«YCII-YII.-l» 
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SUMs01»02 

«ATl»Cl/02 

«»AT2aC2/Cl 

IFd.GT.JSUBUI Pt*IPII*l»*'?ATl-P<II*IRATl-l»AT2»-PII-H*«AT2»/SUH 
TMVs(THa«ll«RArt>rH(II*fPATl>RAT2l*TH(I>ll»RAT2l/SUN 
OPOS «IASH«C0S<THIIM**2>GAHm*PtI>«EP(II«*2«C0S(THIIt l*rHY 
1-SIN <TH« 111 *PYI/ ic’d n**2*C0S<f HI ft I ••2-1.1 
OPON«PV/CCS<THtIII-TANITHtI>I^DPOS 
CS«l./6APtII 
SMal,/ENHI*^2 
OTOS« -CPCK^6S/P I 1 I •SH 
RETURN 
END 

SUBROUTINE CPCINT 

CCHNON/AL/GAR.GEH 

COMPON/eA/ALPI7«5SI .EHINF.UINF 

COMMON/eC/CAHP.PBt4T*RHQO*THEtHH«XHue«Y8 

COMHON/CA/MOOTNI7.SSI .XNC55I 

nOHNQN/CJ/CPI7,95l *CP1I7I .CPXIS6I 

CCMNON/CK/NTHOLE (71 

CONNON/Oe/RETei<*l,IS(i:t 

CONNON/OP/VN(SS) 

CONHON/EC/CPIN.OO 

C0NN0N/EE/E**(S9».GAH(55»«P(59I«TH(55I.V(5SI 

COPHON/EG/EIN,PR,XU 

COPNON/EP/6ANlNF«Hl(7t*RINF 

COHHON/GF/OELV.OVISA*KOUNT3.VISA 

COFNON/CK/OELX 

CONMON/HL/AtPHA.BETA 

C0MN0N/HF/ALPN(7,9SltCPN(7,SSI *CPXN(S9I tENN(95l «GANN(99 l« HN (7 *951 « 
IL*PN(56l*CN(59t •RH0N(95I«RN(99I «THNI55t *TN(99I »UN(59I *XNUN(55I 
CONMON/HN/CHC(2l *CPB(7| ,CPXP(2I tOALDIF 171 *0ALPB(7 t •OOALP8(7I *OELS* 
AEMPI2l*Hei7»«HC(7l«RP(2t«S3A(7),S3Bm«S304 7l,M00T8(7t*XP(2l 
C0PN0N/PC/H(55I.X(99) 

COPNON/PO/JCHEN*NSP.T (991 
C0NN0»i/CA/H(7,99|«(|(95ltRNQ(95)«XNUI95t 
C0PNCN/CS/Rtl0P(2)«M00T(7*55l*M00TC(7l«NP(2l*XHUP(2l 
CONNON/SC/BQN(9S)*OALPN(7«99)*CeQN(55l«CCPXN(95l*OOALPN(7,59t « 
10TAUN(S9ltTAUN(99l 

CONNON/SS/ALltAL2tBQl«B02*Cl*C2*CHl*CH2*0O1.0B2*OOl*0O?tOTl*0T2*OV 
Al«DV2«PXl,PX2,TAl,TA2,THl,TH2ttfl*V2*YlfV2 
C0PN0N/ST/I13«IREGI*K tKFIRST.KXKOtPSTAR 
C0PNCN/TS/0VXSB.0VISCtIFS«NMP*VIS8.VISC 

COMNON/TU/eO(99l«DALP (7,951 ,OBQ (991 , 0CPX(59I *00ALP(7,55I *OT AU (551 . 
tTAUI59l 

C0NNQN/TV/ALPP<7,2l,BET*aaP(2l*0ACHP(7,2l«0ALPP(7*2l*0aQP(2l* 
lOCPXP (2 1 «00AtPP (7*21 *0TAUP(2l*0rCHPI2l,GANP (2 1 *PP(2 1 *qP (2 (* 
2TAUP(2l*TPP(2)tTP(2l«VP(2l 
COMNON/VT/QACH(7*99t«OTCH(99l,OVXSO*VISO 
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C0NH0N/«H/IC0NT*IEtl0tKT«TH8PN*XBP»l 
CONH ON/MV /NPT S*RE*XaP»XJ 

OIMENSIOP ALPSSI7l,OALPSS(7l,OOALPS(ri*HSS<n.CPSSm 

OIMCNSIQK 0UNCHP(7|,0ACHSS(7l 

IFCIENO. EG.lt GO TO 601 

EH3SXP3 (ALPHA, BETA, THTKt tTHNU 1 1 

XN(LI»XnPH 

IFdFS.EG.l.ANC.L.EQ. ISfHHMM £(13sTAN(8£T8(NHHn 

IFdFS.EQ.Z.ANO.L.CO. ISIMMNII £M3s.S« dANtBETBIHHHM »TAN(8ETt I 

YNat«V(Kt»OELX»EN3 

IFdCCNT.EO.il GO TO 601 

KP«1 

EHEPsXHl (ALPHA .BETA, TH(KT»1 I ,KHU(KT»1 I ,THN(LI .XHUNILI t 
EH2L3XN1(AL»HA.BETA,TH(KTI.XHU(KTI,THN(LI .XNUN(L) I 
351 EN2K»0.5*teM2L*EM2Pl 
XP(KP|sXBP 

VP(KPI=yNIH-0ELK*EN2K 

KIP2=0 

IF(YP{XPI.LT.T(KT*ll*l.E-25.ANO.YPIKPI.GT,V(KTl-l,E-05l GO TO 2C1 
4150 FCRHAT(3I5.6E13.5I 
MRITE(6, 91911 

WRITF(6«415CI KP .L.XT . ALPHA,YP(KPI ,YN(LI ,THN(LI ,XHUN(L) , PN (LI 

HRITE(6«11H ' Y(<TI,MKI,OELX,eN2K,XBPN,TH(KI .XHU(KTI ,XHU(KT»1I. 
IDELS, THBPN.OELY.Y (KT*1I 
WRITE (6.20001 

20(i0 FORNA' 97H Y LOCATION OF CHARACTERISTIC ON ORIGIONAL OATA LINE IS 
lOUTSIut CF BOUNDING STREAMLINES IN CPOINT I 

STOP 

201 RATB*(YP(KPI-V(KTII/(Y(XT*1I-Y(KTII 
EN2»EM2L*RATB*(EM2P-FN2LI 
YBT«YP(KPI 

YP(KP|»YK(LI-0ELX*EH2 

IF((AeS(YP(KPI>VBTI/A6S(Y(KT*ll-Y(KTIII.LT.O.bllGO TO 2.2 
KIP2«KIP2*1 

IF(KIP2.LE.20I GO TO 201 
WRITE (6,91911 
9191 FORHATdHlI 

WRITE (£,20011 

20C1 FORNAT(56H UNABLE TO LOCATE Y LOCATION CF CHARACTERISTIC IN CPOINT 
II 

STOP 

2C2 RAT9*(YP(XPI -YIKTI 1/ (Y(XT»1 l-Y (XTI I 
QP(KPI--5CIKTI*RAT8*(<J(KTFH-C(KTII 
PP(KPI»P(KTI*RATB*(P(KT*ll-P(KTl I 
TP(KP| »T(KTI*RATe*(T(XT*ll-T<KTII 
THP(KPI»TMXTI*RAT8*(TH(KT*1I-TH(KTI I 
TAUP(KPI«TAU(KT|4RATB»(TAU(KT»1I-TAU(KT|I 
BQP(<'pl«8C(KTt«RATe*(aQ(KT»ll*8C(KTII 
0CPXP(KPI>qCPX(KTI*RAT8*(0CPX(KT»ll-0CPX(KTII 
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OTAUP(KPI>OrAU(KTI«<tAie»IOTAU(KT»ll-OT«UfKT)l 

08QP(KPtaCBQIKT)»RATe*(08Q(KT*ll>0BQ(KTn 

OTCMPIKPI*DTCMtKTU9AT0*<OTCMIKT*ll-OICH<KTII 

CPXPIKPlsO.O 

MPtKPI»a.O 

CHCfKPtaO. 

CALL THeRPOITPIKPl.HP.CPBI 
no 4020 4=1, KSP 

'LPP(4«KPI=ALP(4, <T| »RATn*(ALP(J, KT«tl-ALP<4* KTII 
'IALPP(J.KP)=0ALP<4* KTMRATe*(OALP(J. KT«H •OALPC 4* KTII 
00A(.PF(4,KPI=0CALP(4, KTt »RAT8* <00ALP ( 4 , KT»1 l*00ALP(4« KTII 
CHC<KPI=CHC<KPI*CiLPP«4,KPI*CP8<4l 
CPXPI KPI«CPXPIKPi ♦AtPP(4»KPI •CPBI^I 
MPiKPIaPPIKPI* AtPP«4,KPI/MTM0LEC4l 

0ACHP<4,KPI=0ACH(4,KTI»QATa*(0ACHI4,KT«l|-0ACH(4fKTII 
IF(K».NE.2I0C TO 4020 

MOOT8(4lal>OOT(4,KTI »RATfl* <H00T( 4. KT^l I -HOOT (4,KT| I 
4020 CONTINUE 

HPCKPI=1./HP(KPI 

PPIKP|=9C/HP«KPI 

CAMP<KPI=CPXPIKPl/<CPXP»KPI-RP<KPI/CPINl 

RK=l./RPIKPI 

RM0P«KPI=PP(KPI*HPIKPI*CEW/TP<KPI 

FMPIKPI=CPIKPI*EHIN?»SQRTIGAR/G4MP<KPI*PK/TP<KPII 
XNUP(KPI=ZNU(FNP(KPI I 
IFIKP.FC.2I GO TO 6n 
KP»2 

IFdFS.tC.QI KT*L 

CM2P=XP2<ALPHA,BFTA.TH(KT*1I ,XHUIKT*ll,THN(Ll,Xt1UN(U I 
CN2LaXP2<ALPHA,8ETA,TH(KTI,XPU<KTI,THNai,XHUN(LI I 
GO TO 351 

C GET ALL the PPOPERT.es AT THE C POINT 
6C1 CONTINUE 

IMIFS.EQ.OI GO TO "6C0 

y ss«y CKi 

P SS»P <KI 
Q SS«Q <KI 
T SS*T <KI 
TH SS=TP IKI 
80 ?S«fiC IKI 
W SSsM (Kl 
TAU SS=TAU IKI 
RHO SS=RHC IKI 
CP' S«=CFX IKI 
080 SSsCOG IKI 
OTAUSSaOTAUIKI 
OCPXSS«CCPXIKI 
CTCHSS«0TCHIKI 
00 1555 4«1,NSP 
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CPSSt J»*Cf>(J*KI 
HSStJIcHiJ.iO 

OALPSSUMOALPlJtKI 
OACH^S(JI«OACH|J*Kt 
15SS OCALPS(JI>OOALPO*KI 
CALL CPClMIK.Lt 
KT*K 

S6;C CONTINUE 
CH20*C.O 

00 A330 jsl. NSP 
CH20sCP2C»0ALP(J.KI*CPIJ.KI 
40 3C CONTINUE 

IFieETA.NE.C.riGC TO 4036 

TAUNfLMTAUtKI 

OCNlLlsOCtKI 

OCPKN(LI»OCPXtKI 

CTAUNtLI*CTAU(KI 

OOOKIU>CPQIK) 

THNILIsTNCKI 

CPKKILMCPXtKI 

TN<LI=T<tU*OTC*<tit> 

HN(L>>‘M(K) 

Ch2C*CH2C 

00 4?»5 J=l. NSP 

OAtP»Mj.L>*OALS>IJ.Kt 

OOALPNIJ.LI^OCALPU.KI 

HCI JtsHIJ.KI 

HDOTC(J)=N0CTNIJ*L> 

4C35 CONTINUE 
4CJ6 CONTINUE 

IFIOETA.EC.C.Ot GC TO 3C? 
CH2CSC. 

DO 331 Jsl.NSP 
HCIJtsHN(J.K) 
NOOTCUIsKOOTNU.LI 
3?1 Cm 2C=CH2C*OALPN(J.U ’CPNU.LI 
3C2 CONTINUE 
Wl*tf ItA 
V2»*ISC 
QVlsOVISA 
3tf?«0VISC 
TAlsTAUPIlt 
TA2sTAUN(LI 
0T1«0TAUP<1> 

0T2*0TAUML» 

0Ql>8QPIlt 

BQ2sSGN(LI 

Tl*YPm 
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t2atN*H 

TM2«TI-N(LI 

Cl*CPXPItl 

C?*CPXN«LI 

oei«oeQPii) 

oes«oeQKiL» 

PXt*OCPXPItl 

PX2«OCPXML» 

CHt*CHC(i> 

CH2«CP2C 

IFItCONT.EQ*!! GO TO 43C9 

S1«*S1IXJ,PE> 

S2««S2(XJ*9e» 

%3Q9 Vl«VtSO 

ovis aviso 

TA1«TXUIK) 

OTisOUUINI 

8Ql«aC(Kt 

TNl«rHCK} 

CisCPXIKt 

OBlsOeOtKI 

PX1«0CPX(KI 

CN1SCN20 

S10>SUXJ«RFt 

S20«S2IXJ*Rei 

XPlICCNT.lQ.l) GO TO 6927 

If tL.EG.NPTSI GO TO 6927 

vt«vise 

ovisDVise 

TA1«TAUP(2) 

0Ti«0TAUP;2i 

aQi-eopcE) 

Vl*YP«2i 

TH1«TNPI2) 

Cl«CPXPf2l 

oei«oeQPi 2 i 

PXt«0CPXPI2l 
CNt*CHC(2l 
S18>81UJfRFI 
S2B«S2(XJ«REI 
6927 CONTINUE 
SJAT>0.0 
SSST-e.O 
SSOT*C.O 

00 9090 U>1« N8P 

AL2«0ALPNfJ*LI 

002>00Al.PNIJ*U 
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IFCICGNf •EQ.i) GO TO 4311 

Vl-WISii 

OVl^OVISA 

ALl=OALPPtJ,l) 

OOlsOCAtPFCJ.l) 

901«BQP(1I 

THlsTHPfll 

Y1*YPI1I 

SiA(J|sS3<Xa»REl 
S3AT=S3A1«S3A C Jl /MT^OLE ( Jl 
4311 Atl^OAtPU«KTI 
001*OGALP(J»Ktl 
V1=VISC 
□Y1=0VISC 
QQlsBCCKTI 
THlaTMKtl 
Y1=YIKTI 

S30CJI-S3CXJ«PEI 
S3CTsS3CT«S30CJI/»iTN0L£UI 
IPCICCNf.EQ^llGO TO 494C 
ira.EQ.NPTStGO TC 4040 

vi^visa 

ovi^ovise 

AL1^0ALPPfJ«2l 

0Cl«0CAtPFfJ«2l 

SQ1«BQP<2I 

THl*ThPI2l 

Y1»YPI2I 

S3PIJ|sS3CXJtPEl 
S3et*S3eT#S3BUI/WTHOtEUI 
4t4P CONTINUE 

IMICCNT«FQ.l) GO TO 6429 
GAHB«GAPF(1» 

PBsPPin 

OesQPCll 

PHOBsPHCFllI 

THe*THPIll 

we^MPfll 

XHUBsXHUPTl) 

YB*YPfll 

Al«FlfLI 

A2»F2a,SlAtS2A«S3ATI 
IFUCHEH.EQ.il GO TO 7252 
A3«0. 

GO TO 72S5 

7252 00 1712 J^ltNSP 

1712 OUHCHPUUfOACHPUtll ♦0ACHU*KII/2. 
OTCNP ( II » t OTC HP 1 1 1 ♦ OTCH I K 1 1 /2 . 
TPl«mLMTpm^CTCM(tll/2. 
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A3«F3<TF1.0rCHP<l)*TP(ll.TI«a)»TNPai,THN<LI*0UHCH(».MP(ll.HN(LII 
fZSS /U>«f4(eeTA«l.,XHuPm •THP<tl*XHUNfU«THK(Ul 
A2s<A2*A3l»A(» 

IFAL.CQ.FFTSI GO TO 6429 
GAHesGAPPI?) 

PSsPPtZ) 

OesQPCZI 

RH0S«RH0P(2I 

THQsTFPtZI 

HesMPIZI 

xHue«xpuP(2i 

YB«TPI2I 

9l=Fl»U 

B2=F2<L.S19»S2P.S38TI 
IFIJC^FF.FQ.ll GO TO T253 
93*C. 

GO TO 7256 

7253 00 1713 J*1.NSP 

1713 OUKCHFtJI=(OACHPf J,2I »OACH( J,Ktl/2. 
OTCHPt2lsfOTCMP(2l*OTCHlKlt/2. 

TP2=ITaMTP<2l*0TCH«LI»/2. 

033F3ITP2*0TGHP(2»tTP(2l*TNILI tTHP(21«TPN(LI«0UHCHP,WP(2l,HN(LI I 
7256 9<i3F4fe6TA»>l.«XHUPI2t«THP(2>*XNUN(U (THNILII 
B2s(e2»83l«n«i 
6629 CONTINUE 

IFdEN0.EC.5l GO TO 630 
THNILIsTeePN 

PN aisPPm*<THPm-THNINPTSI-A2*CXN«NPTSI-XPIim/Al 
GO TO 631 

630 IFdCCNT.EQ.OI 

1PNIU*IA1*PP<1I*01*PP<2I*TMPI1I-THPC2I“ 

P(A2»e2l«XMU»A2«XP(ll»e2*XP(2ll/(Al»eil 

IFdCCNT.EQ.il 

1PNH.MPSTA9 

IFdCCNT.EQ.OI 

lTHNILIaTNPm-Al*«PNdl-PPIlll-A2*<XNU.I-XPflll 

631 CONTINUE 

OCLSa 2. • (XNIL I «X ( XI I / tCOS d H IK 1 1 »COS ( THN (L 1 1 1 
TC«N2=RI<CIKI*CIKI 

IFiaETA.GT.O.ITEPP2a<TEPN2»RHONai«UNILII«.5 

0T*1./TEPF2 

QNIL I «IS1C*0ELS-PNILI *01X1 I *OT*CIKI 
IFISETA.eC.O.CICPXNILIaCPXIKI 

OTCHEPaOTCHfl.l»IPN(Ll-PIKII*(QNILdQIKII/ICPX<KI«CPXNILII*CIN«OT 

DT0IFF*S20*0EtS*£INa2./ICPXIKI «CPXNILI I *0T 

TNILI«T(KI»OTCHEH*OTOIFF 

CPXNILMO.O 

HNILIaO.O 

CALL THERPOfTNILI .HltCPlI 


ORIGINAI! PAOQ ]g 
OF POOR (^UAUTV 
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00 <i«Si J-1, NSP 
OAtOIF(Jt>S30IJI«OELS«OT 

AtPNIJ* LI«ALP<J«KI »OAL0IF(JI»O«CHf J.I.I 

CPNIJvLI XPtUI 

MhtLIaNK (tl»ALPN(J* O/HTN0LE<JI 
CPXNILM ;FXNU.I»ALP*HJ* LI*CPNIJ, L) 

<>090 CONTIFUL 

<tN<Ll*OC/MI|tU 

GA >1K ( i I sCPXN (L I / (CPXN a I -ItN (L I /CPINI 
0RF'«1./RMLI 

PHOI<U|afN<L>*MtltLI*GEM/TN(LI 

ENN(UsOi<(LI*EnikF«SQRTtGAR/GAMN(LI«ORN/TN(Ut 
IFtEHFILl.LT.l.QOail GO TO OOu 
mo XHUKfU«ZPUIEnNILII 
OOO CONTIKGE 

1F(IFS*EQ.0I GO TO 1361 

V K»=t ss 

P <K»*P SS 

Q «KlaQ SS 

T rxMT SS 

H 'KMM SS 

TH IKMTH SS 

BQ < (l«80 SS 
TAU 4Kt«TAiJ SS 
080 IKMOOQ SS 
CPX tKI«CFX SS 
RHO (KlaRMO SS 
OCPXfKMOCPXSS 
OTAUCKI«OTAUSS 
OTCHIKMCTCHSS 
00 1596 J«1*NSP 
AEPtJvKMALPSS Ul 
OALP IJ » K I sO ALPSS U » 

DACH(J«Kl«01CHSStJI 
DOALPtJ.KIaOOALPStJ 
CP<J«KlsCPSSIJI 
1956 H(J«KI«HSSIJ1 
1361 CONTtFOF 
111 FORHATIlOX«9Ei...9) 

RE TURK 
ENO 

suorouttfe space (ISPPI 

COPMOM« mC/IBOO. PIP 

COPR I/AL/5AR.GEH 

C0> jK/AX/jSUBLtJSUau 

rCHN0N/BA/ALP<7«55».SHlNF.MlNF 

C0MM0N/B0/XMASS(9*:| 
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CONHON/CH/MOOTNtr.SSt ,XN<SSI 
COHHON/C J/CP C 7 « I • CP 1 ( 7 I *CPX < SS I 
COHNON/Ck/HTnOLE (71 
COHMON/Oe/BETe ( 41 • IS ( 41 
C0HN0N/0P/VN(S91 
COHMOK/EO/CPIk.PC 

C0PH0N/EF/EH(55I«GAH(9SI«P(5SI.TH(S5I *V(55I 
CCKHOk/EP/CAHINF.Hl (71 «PINF 
COHHON/GK/OcLX 

COPNON/HP/«LPM(7,55I,CPN(7,5SI .CPXN(S9I .EHNISSI •6AHN(SSI«HN(7*S91 « 
IL.PN(S5I,CN(55I«<(H0N(SSI«RN(S3I,THN(9$I,TN(99I«MN(55I«XHUN(5SI 
COPNON/PG/M (SSI , X (S5 I 
CSHHON/PC/JCHEN*NSP.r(SSI 
C0Pk0K/CA/H(7,SSI .<1(SSI «RH0(S5I .XHU(SSI 
C0kH0k/9C/R(SSI 
CCHMCk/RC/APa.APl*AP2 

C0PH0N/SG/3QN(SSI «0ALPN(7*SSI .QQQK(SS1 •CCPXN(5S1*00ALPN(7*5SI * 
10TAUN(SSl.rAUN(5Sl 

CCHNCA/S1/I13.IREGI.K •KFIRST.KKKG.PSTAR 

C0HH0N/T(j/8Q(SSI.CALP(7.SSI .08Q(SS1 *OCPX(SSI »00ALP(7«SSI .OTAUCSSI . 
1TAU(SSI 

C0ni10N/UV/tUl«IERR*IPRESS«IPRESU«ISUe 

COPBON/MV/NPTS«RE*XaP*XJ 

COHHOk/kX/ APRESS* APRE SU 

CCMNOK/XY/APRS*APUS*OELTAV,EeOOS«iaODS. INTACT. IPRStIPUS.ITVP* 
ljeOOS«MHAX«RHEAT.X<2*XK4.YeOT.VTP 
CCPN0N/VX/A80n«*ePRESS*CPRESS 

COHNON/7V/AeOO«eeOI}.CBOO*ESOO.FBOn,GeOO*lAVE.IPUNCH.J800.KKKKK 

OINENSICN ISN(4>*XHASSN(SSI 

XJ2=1.-XJ 

XJ1«1.*XJ 

OTaOELTAY 

OC 1 Ial.4 

IFdSdI.EQ.Ol GO TO 1 
Na-1 

IF((I/2I«2.EQ.II N«1 

istaism 

0YSaAeS(TdSl*NI-Y(ISI*2*MH 
IF(0VS.ET.2.*0VI GO TO 2 

C********* ACC PT ON COHNSTREAH SIDE OF SHOCK ••••• 

TSN*TS(I1-1 

IF((I/2I«2.EQ.II IS'1aIS(II»2 
00 3 KK«ISN.NPTS 
KaNPTS*ISP-KK 
JsK*l 

CALL SMITCH(J*K» 

3 CONTINUE 
NPTSaNPTS^l 
ISPPal 
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IMJSU8L.CT.ISm t JSUBLsJSum.»l 
iF(jsceu*GT.ismi jsueu>jsueu«i 

DC •* 

IfflSIJI.GT.ISdit ISIJI>ISIJ»»1 
<» CONTINUE 

IFI(I/2I*2.NE.II 

L-ISM-1 

MaISH»l 

<sISH 

RATs.f 


p 

CKI»P 

CU^RAT*IP 

INI-P 

clm 

TH 

IKIsTH 

ILMRftT*CTM 

INI-TH 

an 

X 

(KlaX 

ILMRAT’IX 

IMI-X 

CLI> 

V 

«K»ay 

a>^RAT*IY 

IMI-Y 

(til 

Q 

CK)aC 

(LI*RAT»(Q 

INI-Q 

ctii 

T 

iKIaT 

tLMRAT*CT 

IHI-T 

(Lll 

PG 

(K)aeg 

Ctl^RAT^CPQ 

INIoBQ 

(til 

TAU 

(K)aTAU 

(LMRAT^ITAU 

fM»-TAU 

(til 

oeo 

iKtaCBQ 

fL)^^AT*(ORQ 

IHI-OBQ 

(Lll 

OCPX 

(KIrOCPX 

a»*RAT»IOCPX 

CNI-OCPX 

(til 

CTAU 

CKIaCTAU 

ili»rat*iotau 

(MI-OTAU 

(til 


XHflSSIKI»XMASSlL»»R»T*»KMdSS(M»-XHASS«U> 

CPX(KI*0. 

WIKIaC. 

CALL THERPO(TlK)tHl,CPll 
00 $ J>1»NSP 

ALP fJ*Kl>ALP (J»LI*PAT»fALP (JtMI-ALP UtLII 
DALP (J*K)sDALP tJ.LI *RAT«(DALP IJ*n»-OALP IJtLU 
OOALP(JtK>>DOALPIJtLI «RAT«(OOALPUtN)»OOALPf JtLI) 
HtJ*KtaHilJI 
CPIJtKIsCPlUl 

M(K)«MtKMALP(JtK)/HTPOLE(JI 
CPX«K»*CFXJK»*ALPCJ.K>*CPJ JtK» 

5 CONTINUE 

N(K|al./t>(Kt 

RtKlaRO/H(K> 

GAP<KI*CPX«K»/tCPXK»-RltO/CP IN» 

0R*1./R(K> 

RHO(K)aPCKt •NIX>«GEH/T(K) 

ENCK»»CCK>*EMINr*SQRT tGAR/CAM|lO*OR/TIK»» 
XMU(K>*2 KUIcN(K) ) 

2 CONTINUE 

YSN*VtISt>«TAN(BETB(in«OELX 

IFIII/2»«2.EO.n GO TO 6 

J>IStII 

K*J»1 

EHPsXH2U.*C.*TH(KI «XPU(IO«0.tO>> 
ENL>XK2(l.*0.«THtJ>.XHUU) t0.*0.> 
EHPl*XHlllktO.»TH(K>t XHUIK> tO.tO.I 
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GO TO 7 

6 JsISfl) 

K*J-1 

EHPsXni(l.*C.*lnMI*XnU(JltO.«G.I 
Et*L3X«11fl.,C.*TH(KI *XrHMKI«Q.*Q.I 
cHPlsKH2<l.»0.*TH(J),XMUtJI,0.,S.I 
enLl>Xn2ll.*0.«THKI.XNU(KI.0..0.l 

7 EH3=XH3ft.«0.*TH(Kl,G.I 
VCNsV(Kt«EN3*0ELX 
YCTsVCN-CELX»«EMPl*EMLll*.5 
VST=YSN-0ELX*IEHP»E1LI*.5 
OYTs-N*IYCT-YSTI 

If (8CTe(ll.EQ.0..««iO.eCTe(2l .EQ.O..ANO«eETP(3».EQ.a..ANO.BETei(»t. 
AEQ.O.IGC TO 777 

IFieETe(ll.GT.0..09.eETem.GT.?.l777.77e 
777 IF<OVT/AeS(V(JI*Y(Kl) .GT*. 111*775 
776 IF<N*IYSN-YCNI.GT.0Y*.5IG0 TO 1 
ISIPN»ISm-N 

IFIISIPN.EQ.NPTS.OR.ISIPN.EQ.IIGO TO 1 
C********* SUFTRACT PT FROM FREE STREAM SIDE Of SHOCK •••• 

775 L=K»1 

00 S K*L*KPTS 
J=X-1 

CALL SMITCH<J*KI 

8 CONTIMCE 
MPTSaNPTS-1 
tSPP=l 

IFTJSUOL.GT.ISIIII JSUOLsJSUBL-1 
IFUSU8U.GT.Ismi JSUBUsJSUeU'l 
SO 9 Jsi.ii 

IFlISUl.CT.ISdII ISUI»ISUI-1 

9 continue 

IFUI/2»*2.FQ,H ISm«IS<II-l 

1 continue 

IFdTVP.NE.ll GO TO 850 
IFTNPTS.LT.MHAXI go to 2121 
II11«1 
IPUNCHsl 
MRITE16.9191I 
9191 FORNATdPlI 
WRITE (6.8511 

851 F0RHAT(7iiH REQUESTEO NAXINUM NUMBER OF FLOW FIELD PTS. EXCEEDED. P 
lUNCH FILE OBTAINEO/97H RESUBMIT RUN WITH REDUCED NUMBER OF FLOW FI 
lELD PTS. OR INCREASE INPUT FOR MAXIMUM NUMBER OF PTS. I 
RETURN 
850 CONTINUE 

IF(Y(1). EC.YBOT.OR.ITYP.EQ.il) GO TO 2190 
IPRESS«l 
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ftPP5SS»P<l»/PIN 

2100 IFCYfNPTSI.EQ.YTP.QP.ITYP.eQ.SI GO TO 2101 
IPPESU*! 

APRESU«Pi*>PTS)/FlN 
21 Cl CONTINUE 

ifinpts.it.mnaxi go TC lOCO 
ISC»0 

00 roc 1=1. <* 

rcc iFiism.KE.o) isc=isc*i 

IFtITVP,fE.2.0R.tSC.NE.F.0R.ISU0.NE.3.)G0 TO 701 
YQ»YTP 

IFCYIFPTSKLF .YiP-2.*0EtTAY»YQsY«NPTSI*0EtTAV 

IK=1 

IXal 

NP»NPTS 

IFHNPTS/2»*2.NE.NPTS» GO TO 702 
8CC J»NPTS*l 
K=NPTS 

CALL SWITCHU.K) 

Y«J»=YG 

NPTS*NPTS*1 

YFUN»IYIJ»*<XJ2*YIJ»*«J»-YIK»*IXJ2*YtKl=XJH/XJl 

PQAV=PHCIK»*QIKI*COSITH«<n 

XMASSUI = XNASStK) ♦ROAW*YFUN 

THtNPTS>»C. 

APRESU*PIKPTS>/FIN 
GO TO (702.903I.IX 
7C2 J*1 

00 7C2 K=3.NP,2 
J»J*1 

CALL SHITCHIJ.KI 
7C3 CONTINUE 

0ELTAY»CELTAY*2. 

IFIIK.EQ.2l 60 TO 704 
NPTS»NPTS72»1 
GO TO 1000 

m IF(ITYP.NE.3.0R.ISm .EQ.C. OR. ISUe.NE.O. OR.ISC.NE.il GO TO 303 
IK=2 

NP=IS(3I*1 

TF((NP72I*2.NE.NP| CO TO 702 
10*1 

IF|Y(1).6E.YBOT*2.*OELTAY*1.E-03» 60 TO 706 

YTsYOOT 

GO TO 707 

7C6 YT«Y(1I-CELTAY 
707 00 700 KK*1.NPTS 
K«NPTS*1-KK 
J=K*1 

CALL SWITCK(J.KI 


OBlGlNAi; PAGE IS 
OF POOR QUAIOT 
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7CB CONTIKUC 
NPTS*NPTS^1 
00 7C9 Isl«4 

7C9 IFaSfII*NE*OI ISCIMlSm^l 
IFflSUa.EC.OI GO TO 710 
JSUBLsJSUEL^l 
jsuBu^jsueut^i 

710 mi*VT 

YFUN«IYC2l«CXJ2«V(?l»XJt*Yf ll»(XJ3»Ytll»XJM/XJl 

PQAYsRHOI2l«QC2l«COS(THf2n 

XHASS CllsXMASSf2l*RQAV*YFUN 

THIIMO. 

APRPS?«Pm/PIN 
GO TO (711t907l »IQ 

711 NP=NP^1 
GO TO 702 

7C4 ISM(3l«NF/2^2 
IO-ISC3l-I$Nf3l 
ISS«ISI3) 

ISC3I-1SM3I 
00 705 K«ISSfNPTS 

CALL SMITCH<J«KI 
7C5 CONTIKUE 

^PTS«NPTS*IO 
GO TO 1006 
303 ICT»ISC 
ITCP*I«PTS 

IFfIS(3KN£«0l ITCP*ISm*l 
IB0T«1 

IFIIS(4I»NE.0I I80TsISt4)^l 
OTY»YIITCPI-Y(ieCTI 

0EtTAY*CTY/Ft0ATUHMAX-(NPTS-IT0PI-ISl4H/2-ICTI 

IB*IBOT 

iSNfii«ism 

ISN(2I>1SI2I 

ISM3)sXSC3l 

ISNl4lxISt4l 

JSUBLNsJSUBL 

JSUBUKsJStBU 

iBe«ie 

IPCG«i 

IFfIS(2).EQ«0) GC TO 501 
I*-IS<2I 
GO TO 502 
SOI IPEG*2 

IFCISUB.EO.O) GO TO 5C4 
IFUSOet^EQ*!! GO TO 503 
ITsJSUBL 
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GO TO 502 
603 CONTIMUE 
iBsjsueu 
iae«ie 

6CG IREGsG 

IFdSdl.EQ.OI GO TO 505 

1T»IS«1I-1 

GO TO 502 

505 WEG*5 

IFIIS(3).EQ.0I GO TO 506 

IT»IS«3I-1 

GO TO 502 

506 ITsNPTS 

5C2 MP=«y llTI-r«iei l/OELTAY 
Lsie 
J 2 »l 

DELS (V( IT I'Y (let » /FLOAT (NP» 
5932 CONTIMUF 

j»iee 

K»I« 

X NOIsX IKI 
Y NIJIsY IKI 
Q N<JMC (Kl 
P N(J|sp (Kl 
T NIJMT (Kl 
M N(JI«M (Kl 
R N(J|30 (Kl 
EM N(J|aEN (Kl 
TH NUIsTH (Kl 
80 N(J|s80 (Kl 

TAU N(JI«TAU (Kl 
080 N(J|sC8Q (Kl 
GAN NUIsCAN (Kl 
RHO NfJIsRHO (Kl 
XNU N(J|bXNU (Kl 
CPX N(JI«CPX (Kl 
OCPXNIJIsOCPX (Kl 
OTAUN(JI«OTAU (Kl 
XHASSMJ|sXNASS(KI 
003108 JJsl,NSP 
H N(J,|«JI«H (JJ,KI 
CP N(JJ«J|sCP (JJ«KI 
ALP N(JJ*JI«ALP (JJ«KI 
OALPN(JJ«JI«OALP (JJ*KI 
DOALPF(JJ«JI*OOALP(JJ.KI 
3108 CONTINUE 

GO TO (2201«290AI,JZ 
2201 00 600 KKsltNP 
IsKK«iee>l 
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YNII»llsVN(II»DEL 

6C2 IPIYN(I«ll*GC.V(Lt.ANO.YNlI»ll.LT.V«L«in GO TO GOl 
L*L*1 
GO TO 602 

661 PftTs<VNtI«ll>Yai l/(Y<L«ll*Y(LI > 

IPdT.EC.JSUBUI GO TO 1200 
p«<a«ii-f ai«RAT*cpa*ii>paM 
THNII«llsTHIt.l»RAT»CTH<L*ll-TH(l.l) 

GO TO 1201 
1200 CONTINUE 
YY»YN<I*1I 

PN^I♦1»*AP0♦YY•*AP1♦YY•AP^I 

THNtI«l»<C. 

12C1 CCNTIKUE 
M=t»l 
<=I»1 


X 

NiKisX 

fll^RAT^IX 

<HI-X 

an 

Y 

N(KI«Y 

<LI^RAT»ir 

CHI-Y 

ftn 

Q 

NCKUG 

lU^RAT^n 

<M»-Q 

CUI 

T 

NCKMT 

lU iRAT^IT 

(Ml-T 

tun 

eo 

NiKi^eo 

CU^RAT»(eO 

iHi-eo 

tun 

TAU 

NfKisTAU 

CU^RAT^ITAU 

(HI -TAU 

tun 

Deo 

NfKIsOBQ 

ai^RAT^tDPQ 

(HI -080 

tun 

OCPXRtKtsCCPK 

aMPAT^inCPX 

(HI-OCPX 

tun 

OTAUN(KlsCTAU 

CU^RAT^fOTAU 

(HI-OTAU 

cun 


XMASSMKI>XHASSai*RAT*(XHASS(RI>XMASSail 

CPXN(K»sO. 

WNfK)«0. 

CALL TPERPO(TNlKt.Hl«CPtt 
0055 Jsl.NSP 

ALP N(J«K»sALP tJ«L)«RAT*(ALP (J.HI'ALP U«LM 
OALPNU«X)>OALP (J.LI«RAT«(DALP U.m-OALP UtD) 
OOALPMJ«KI«OOALP(J«L»»RAT»(COALP<J*NI>OOALPU«UI 
HNUf K|3Hl(Jt 
CPNU«KMCP1CJI 

MNIKI«MMK|«ALPRtJ»KI/MTHOLE(J) 
CPXN(KI«CPXN(KI*ALPN( J»KI«CPN(J*KI 
55 CORTINUE 

MNIK>3l./MN(K) 

RN(KI»RC/hN(Kt 

GANN(KI«CPXNIKI/(CPXN(KI>RN(Kt/C.>lN> 

RKel./RNIKI 

RHCN (Kl «PMKI*MN(Kt •gem/TNI Kl 
EMNIO QR(X)*ERINF*SQRT(GAR/GAMNIKt«RK/TN(KII 
IF (ENN«KI.GT.l.t 
1XM0N(KI«2MU(ERN(K)> 

6C0 CLRTIAUE 

GO TO (22C0«603«604(&C5t6<'6t*IREG 
2209 ISM2t«I«l 
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iesISI2Ml 
leo^isNceui 
GO TO 501 

6C3 

ie«jsueu 

jsueuNsJSu<^tNfjsueu-JSURt 
IBBsJ<(ieUN 
6C4 CONrI^UE 
GO ro 504 
6C5 ISNClMl^e 
iQ«ism 
IBe>ISMll 
GO TO 505 

6C6 IFIISf31.N?*0» IS^!(3I-I♦^ 
IFC1SC3KFQ.0I NPTS«Ifl 
NP*NPTS 

IFfISCSI.EO.OI GO TO 2eC3 

IO«IS(3l-ISNC3l 

ISS»1SI3> 

00 23G4 K«ISS«NPTS 
J»K*IO 

CALL SNltCHfJ^KI 
2314 CONTINUE 

NPTS»NPTS-IO 

NP^ISNm-1 

iFcisue«eQ*0) go to 2903 

JZ*2 

J2IC<JSueL«i 
jzi«jsueLN-i 
29C4 JZI>JZI^1 
JZIOSJZIC*. 

IFf JZI«EG.JSUOUNI GO TO 2903 

leesjzi 

ie*jzio 

GO to 5932 
2903 CONTINUE 
22C1 00 22C4 l«ldOT#NP 
J«I 
K»I 


X 

CJI«XN 

IKI 

y 

CJ)sVN 

CK) 

0 

(JI«0N 

(Kl 

p 

(JI»P 

NCK) 

T 

<J>*T 

NCKI 

N 

(JI*N 

NCK) 

9 

IJI«9 

NCK) 

£P 

CJI*EH 

NCK) 

TM 

IJI*TH 

NCKI 

OQ 


NCK) 
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TAU UI»TAU NIK) 

080 IJIsOSQ }|KI 
GAM CJUGAH NIK) 

RHO IJIsRHO NIK) 

IJ)«XHU NIK) 

CPX IJ)3CPX NIK) 

OCPX IJ)sQCPXNIK) 

OTAU IJ)sOTAUN|K) 

XMASSIJ)sXHASS^*IK) 

DOlilOa JJslvNS) 

H IJJtJ)-H NIJJ*K) 

CP IJJ«J)sCP NIJJ«K) 

ALP IJJ«J)»ALP NIJJ*K) 

OALP IJJtJ)sOALPNlJJf K) 

DOALPIJJfJ)sOOALPN(JJtK) 

4108 CONTINUE 
2204 CONTINUE 

00 607 I«i«4 
6C7 ISlIlsISNII) 

JSUeL«JSIi8LN 
jsuBu^jsteuN 
1000 CONTINUE 

IFIITVP*EQ«3) GO TO 903 
If IYINPTS)«EQ.YTP) GO TO 903 
YQ«YTP 

IFlYINPTS)*tE#YTP*l.*OELTAY)YQ«VlNPTS)^OELTAV 

IX«2 

L*NPTS-1 

Mat«l 

IMAeSIPIH)*Pa))/PlL)*.Q01) 9C0f900f 630 
9C0 IPIABSIQIN)-Qa))/QIL)-.C01) 90lf9Cl«600 

901 IPIABSITIP)*Ta))/Ta)«.QCl) 902f902f800 

902 lFIABSIAtPI6tH)*AlPI9«L))-«001«ALP(3tU) 9e3f903t60C 

903 iriYI|).EG#YBOT) GO TO 9C7 
IPIITYP.EC.4) GO TO 907 
L»2 

H«3 

• 0*2 

If %M«)*6E.YBOT«DELTAY^1*E*03) GO TO 910 

YT»YBLT 

GO TO 912 

910 YT8Yli)*OELTAY 

912 IP|ABSIPIN)*PIL)|/Plt)-.031) 904«904«7C7 
9G4 iriABSIQIH)-QIL)l/OIL)-«001) 906»905«707 
909 IMABSITIN)*TlL))/TfL)*.001) 906f906«707 

906 IPIABSIALPl4fH)*ALPl4a))*«001»ALPl4fL*' 9C. *9079707 

907 CONTINUE 

IFIITYP.E0.4) 60 TO 2102 
IFIYID.NE.YBOT) GO TO 2102 
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iBcosieccs 

A80Qs«eCCS 

APPCSSsAPfiS 

21C2 IFIITVP.EC.3I GO TO 2103 

IFtVtMPTSl.NC.VTPl 50 TO 2103 
J800-J8CC« 

CBOOsCeCCS 

ipp£su*:pts 

APPESU^APUS 
21C3 CCNTIAUF 

lFIV(ll.EG.VeOT.A»iO.ITYP.EQ.3l ITVPcl 
IF(T(ll.EG.VBOT.ANO.TTVP.EQ.2l ITYP=A 
IF(Y(»iPTSI.EQ.VTP.ANO.tTVP.EC.AI ITYP=l 
IF(YIKPTSl.EQ.YTP.ANn.ITYP.EC.2l 1TYP*3 
2121 CCNTIKUE 

IFdSUB.EC.il RETURN 

JSU8L=NPTS»1 

JSU8U=l«PTS»i 

RETURN 

End 

SUBROUTINE RSET 

CONN OK/Al /GAR. GEM 

CCNN0N/BA/ALPd.99l .ENINF.HINF 

C0NN0N/8G/XHASST9SI 

COPNON/C J/**® «7,S5I.CPldl.CPX<55> 

CONNON/OP/VN(SEI 

CONNOK/EC/CPIN.RC 

C0NN0N/EF/ENISSI,GAN(55I.P(5SI«TH<SSI .YTSBI 
CCNHON/EP/GANINF.HI 17 1 .PINF 

CONNON/Hr/ALPN17,S9l«CPN17.SSI .CPKNlSSt .ENN 1591 .GANN! 991 ,HN 17.5SI « 
IL.PN199I.QN199I.PI-ON195I.RN199I .THN195I.TN199I.HN159I.XNUN15SI 
CONNON/rR/TNBP.YeP.VNPN 
CONNON'/PC/M199l .X199I 
C0«'N0N/FC/JCHEN*N$P,T159I 
C0NN0N/QA/H17.99I .Q195I .RH0I99I .XNU199I 
C0NN0N/PC/R159I 

CONHON/SC/SQN199I .0ALPN17.S9I .0BQN159I .OCPXNlSSI.OOALPNd.SSI # 
1DTAUN199I.TAUN199I 

CONNON/TU/BQ 1951. 0ALP17.95I *080 1591 tOCPX 1991 «00ALP<7«55I *07 AU:55I« 
1TAU19SI 

C0NN0N/VM/IC0NT.IEN0.KT.THBPN.X9PN 

CONNON/UV/NPTS.RE.KBP.XJ 

00 Site lal.NPTS 

THd)»THMII 

XlIMXBPN 

Y1II=YN1II 

QIIMQMII 

PlII^PNllI 
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rtiisTMii 

XKUIII=XrUN(II 
TAUlIlsTAUNf II 
‘'QfllsPCMIl 
OCPXf II>CCPXN(II 
OTAUCIMCTAUNfll 
onodtsoecNm 
oo<*a9Cj-i« NSP 

ALPtJ* IlzALPrUJ. II 
0ALP(J« II-OAtPNU, II 
OOALPU. IlsDOALPMJ, II 
CPU. IMCPNCJ, II 
HU. II^HMJ. II 
4a9C CONTIIcue 

HIIIsMMIl 
PdMRMIl 
SANf IlsCAPNdl 
CPXdlsCPXMIl 

stio CONI H ue 

XJlcl.»XJ 

IMYdl.EC.P.I XHASSIllsO. 

00 10 I*2.N«»TS 

vFUN=lvdi*d.-xj*xd i*xji- vi:-ii«d*-xj»rd-ii«xjii/xji 
RQAV=UtHCdl*Qdl«COSdHdlM RMO d-1 1 •Cd-ll *005 dHd-1 1 1 1 /2. 
XHASS dl = XriASS d>l I »RQAV*rEUN 
10 CONTIKUE 

00 A*»C9 I=1.NPTS 
CPXlIlsC. 

CALL THERNOddI .Hl.CPlI 
00 «(>10 Ki=l.NSP 
CPII^.IMCPldAI 
Hd4.IlsHldAI 

0410 CPX(IlsCPX(II»ALPd4.II«CPld4l 
RH0dl=6Eh*Mdl*Pdl/Tdl 
GAMd|sCPXdl/fCPXdl-Rdl/CPlNI 
RI«l./Rdl 

EHdlBCm’EHINF^SOPTICAR/CAHddRI/Tdll 
IFIEHfll.LT.l.OCOH CO TO 8409 
XMUdl=2H0lEMdll 
84C9 CONIIFUE 
XBPaXOPN 
Y8P*y INPTSI 
THeP^THINFISI 

IFfVaPN.EC.VNtNPTSII RETURN 

yOP*Y0PN 

TMeP>THeFA 

RETURN 
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€N0 

SUeftOUTIKC SHCAAKCFFtVISOI 

COHMOK/AC/ieOO.PIN 

C0MN0N/ea/ALP(7.S5l .ENINP.HINF 

COPN Ott/C J/CP ( 7 • SS I • CP U 7 1 *CPX ( S9 1 

C0PN0II/EF/EN(S5I«CAH(S5I.P(56I.TH<S5I«V(55I 

COPMO»l/PC/JCHEH,NSP*Tl$fi| 

C0NN0I«/QA/HI7,S5I .QfSSI •PH0(59) *XMUtSSI 

C0mt0N/SC/aQN(95l*0ALPN(7.5SI *0?QN(SS I *0CPXN(S9I «00ALPN<7«S5I • 
lOTAUNISSI •TAU(i(S5l 

CCPHOK/TU/8Q($SI.OALP(7.SFI.Oeo(6Sl,OCPX(65l*OOALP(7.5SI.OTAUt55l* 

1TAUC5SI 

COPMO»l/NV/NPTS«RE*X8P.XJ 
OIPENSIOK LOCStSI 
KKl-0 

00 IQP K:l«8 
100 LCCS(KI-0 
LASTsKFTS 
LASri«NPTS*l 
LAST2=NP?S-1 

r «LASni*2.*V IIASTI-Y (LAST2I 
Q tLASntsQ CLAST2I 
Y ILAST1I=T aA^f2) 

CPX ILASTltsCPX ILAST2I 
PaASTl»»P«tAST2» 

THaASTlt:TH(LAST2l 
TAU (lASTIsO. 

80 fLASTlsQ. 

OCPX (LASTMO. 

00 6290 J»1«NSP 
OALP IJ.iASTI^O. 

6290 ALP (J.LASTII-ALP fJ,LAST2l 
00 62V2 K«2«LAST 
0ELY2=Y <K*1I-Y (Kl 
0ELY1*Y *KI-Y tK-l» 

IF40ELY2.lt. t.E-06. OR. 0FLY1 .lt. l.e-C6t CO TO 1301 

SUH«0FLY1*0ELY2 

RATt01sCELYl/0ELY2 

RATIO2»0ELY2/OELYl 

SUsl./SUN 

RHRaRATICl-RATIC2 

TAU4K|sCQfK*ll*RATX01>QtKI*RMR>04K«l) •RATI02I«SU 
00«1./0ELY2 

OTAU4K»»2.*404K*1I*OELY1*SU-04KI*04K-1»*OELV2*SUI/OELY1*00 
eQ(Kls<T(K^ll*RATI01>T4KI*RPR-T(K>ll*RATI02)«SU 
O804K»*2.*«i 4K*1»*0£LY1*SU-TIK» *T IK-l»*0ELY2*SUI/0ELYl*00 
OCPX4K|a|CPX4K*ll*RATIOl-CPX«l*RPR-CPKCK-ll*RATI02l*SU 
00 6291 J*1.NSP 

0ALF4 J.KI«f4LPU*K»ll*RATI01-ALP( J.Kt •RPR>ALPU«K>lt*RATI92l*SU 



118 


OO«tPtJ•Kl>:2.•fALP(J*K*n•0CLrl•SU-ALP(J•Kt»ALPCJ.K-ll•OELV^•SUI 

l/OEtVl*CO 

6291 CCNTI6UE 
GO ro 6292 

13C1 KKIsKKI*! 

LOCS(KKI)=K 

6292 CONTIAUF 
TAU (11=0.0 
cy=»i 2 i-v(i» 

IF(IBC0.CC.1ITAUI1I=CFF*PE*RM0(1I*QI1I**2*.5/WIS0 
90 (1I=u.C 
OCPX (11=0. 

DTAU (1I = (Q ( 21 -C (lM=?./(y ( 21 -V (1II*=2 

IFdROO.EC.ll OTAUdI =4.*(Q(2l>U(lll/DV««2-2.«(TAU(lt*rAU(2M/0y« 
lCTAu(2t 

nPQ (1I=(T I2I-T (1M*2./CY (21-V (lll**2 

IF (I9C0.EC.il 090(11 -4.*(T(2l-T(in/OV**2-2.*80(2l/OV*08Q(2l 
006293J=1. NSP 
CALP (J. 11=0. C 

CCALP (J, ll=2.*(ALP IJ. 21-AtP (J. 1II/(V (21- 
IV (11 •••2 

IF (loco. EC. II OOALP(J.ll-4.=(AtP(J,2l-AUP(J,ll|/oy**2-2.*OALP(J.2 
1I/CV*C0ALF(J.2I 

6293 COMI6UE 

00 101 P=1.8 

IF(LOCS(PI.EQ.OI GO TO 102 

K=tCCS(PI 

t=l 

IF((M/2I*2.NE.HI L=-l 
YNK»y (Kl-y («f*LI 

BC (KI=2.*(T (Kl-T (<»UI/yNK-e2 (K*LI 
TAU (KI=2.*(Q (Kl-Q (K»LII/VNK-rAU (K«LI 
OCPX (KI=2.*(CPX (XI-CPK (K»LII/VAK-OCPX (K*LI 
OTAU (KI=2.«(TAU (KI-TAtl (K*LI I /yNK-OTAU (K»LI 
OnO (KI=2.«(8Q (KI-9Q (K«UI/YNK-ORQ (K«LI 
00 103 J=1.NSP 

OALP (J*KI=2.*(AIP (J.KI-ALP ( J .K»LI I /VNK-OALP (J.K»LI 
1C3 OOALP (J.KI=2.*(0ALP (J.KI-OALP ( J«K*LI I /VNK-OOALP (J»K»LI 
ICl C0NTI6UE 
102 CONTINUE 

OC 70C0 I=l«LASr 

rAUN(II=TAU(II 

eQN(II=*>G(II 

OCPXNdMCCPXdl 

OTAUNdlsCTAUdl 

OPON(II«CPQ(II 

00 70C1 J=1.NSP 

OALPN(J.II«OALP(J*II 

reel ooalpn(j,ii=odalp(j.ii 
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TOCO CONTINUE 
NPTS«t*ST 
OETURH 
END 

SUOROCTIKE SNE«R2(CFF.VIS0I 
CONN ON/ «C/ leOO » PIN 
CONNON/OP/VNtSSI 

C0NN0N/NT/«*wPNI7,SSl«CPN(7«5SI*CP«NiS5l<ENN(5SI.6AHN(35l*HN»Tt?n« , 
ILS«PN(S»I.QNI5SI«RH0N<55I.RN(3SI.THN(SSI«TN<56I,HN(55I.XNUNI5SI 
CONNQK/ PC / JCNEH .NSP • T I I 

COHNON/SQ/eQNI5St«OaLPNI7,SSI.03aN(S5I.OCPXN(SSt,OOALPNf7,5SI • 
lOTAUNtSSl.TAUNtSSI 
C0HN0N/NV/NPTS*RE«XNP»XJ 
OINENSICN LOCSISI 
KKIaO 

00 IOC Ksi,0 
ICO L0CStK**6 
LASTskptS 
tASTlaNPTS*! 

LASTZsNPTS*! 

VN CL ASTI » 32 .*TNCLASTt -TN CLAST2 1 

ON CL ASTI IsQNC LAST2t 

TNCLASTtl3TN(LAST2l 

CPXNCLASTll3CPXNCLAST2t 

PNCLAST1I*PNCLAST2I 

THNCLAST1|sTHNILAST2I 

TAUNCLASTIsO. 

eONCLASTlsO. 

0CPXNCLASTI>0. 

00 3001 J«t«NSP 
OALPNiJ.LASTIsO, 

3001 ALPNCJ*LAST1|3ALPN(J,LAST2I 
DO 6002 K32,lAST 
OELY 2s YN (Xni -VN C K) 

OELTlsYNlXl-YNCX-ll 

IFCOCLV2.LT.l.E*Oe.OR.OCLYl.LT.l.E-e6l GO TO 1301 

SUP«0ELV1*0ELV2 

RATI0isCElVl/0ELV2 

RATI02S0ELV2/0ELV1 

SUsl./SUN 

00»1./CEIT2 

RHRSRATIC1-RATI02 

TAUNfX|sC0N(K4ll«RATI01>QN<KI*RMR«QN(K>ll«RATI02l«SU 

OTAUN (K Is2,s CQN CK«1 1 •0ELV1*SU-QN(KI 4QN CK-ll •0ELV2 «SUI /OELYISOO 

00Nf<l«(TMK»ll«RATI01oTN(KI«RHR>TN(K-ll*RATI02l«SU 

0eaN(K|32.«ITN(K*l|S0ELVl*SU-TN(K»4TNCK>l)*0ELY2*SUI/0ELVi»0O 

OCPX N f K 1 3 CCPKN f K«1 1 »R ATIOl'CPXN (K I «RNP>CPXN f K*1 1 *RATI02 ) *SU 

DO 4001 Jsi.NSP 

0ALPNUtX|sfALPNtJ«K*ll«RATI01*aLPN(J*KI«RHR>ALPNf J«K*1I*RATI02I 


rgsKiINAS PAGE IS 

OT POOR GDAiire 
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1*SU 

OOALPKU.K|x2.*(aLPN(J.K»ll*Oei.Tl»SU-ALPN(JtKI«ALPN(J«K-ll •0Et»2 
l*SUI/0CtYl»00 
AO 81 CONTIKUE 

GO TO 6002 
13tl KKT=KKI*l 

LOCSCKKI)=K 
60 C 2 CONTIKUC 

rSUNdlsC.O 

OV*lfM2l“VNfll 

IFIieCO.€C.HTAoNU>*CFF»Pt-*RHONlU*aNlll*“2*,S/VISO 

BQNdlsO.O 

OCPXNUIxQ. 

0 TAUNfll=(QNI 2 l-QN(ll l• 2 ./(VN(^l-VNIln••^ 

IFdeCO.EC.lt 01A'JN(ltsA.*fQN(2l-CN(lll/0Y««2-2.«(TftUt(lt»TAUN<2l 
1I/CY4CTACM2I 

0B0Nllla|TNdl-TNIlH*2./»YN<2l-YMlM**2 

IFfIBOb.EC.lt OeOM ll=A.*(TNf 2 l-TK(lt t/CY*« 2 - 2 .* 33 lf( 2 l/ 0 Y» 0 eQN< 2 l 
OOA 082 JS 1 * NSP 
OALPNfJ. llsQ.e 

OOALPMJ. ltx2.«fAlPN(J, 21-ALPNfJ* Itt/fYNfZt- 
lYNflt t**2 

IFIiaCO.EC.lt COAlPN(J«ll=A.*IALPNIJ» 2 l-AtPNIJ,lll/OY** 2 - 2 .*OACPN 
If J. 2 t/CY*C 0 ALPNIJ« 2 t 
AO 82 CONTIKUE 

00 181 Psl .8 

IFfLCCSfHt.EQ.Ol GO TO 1G2 

K^LOCSfet 

L*l 

IFf fM/2t«2.NE.t1t L = -l 
YNK=YAfKt-YNfK*Ct 

BQNfKI«2.*fTNfKI>TNtK»LII/VNK-8Q NfK»Lt 
TAUNfKts2.«tQN(Kl-QN(K»Lt l/YNK-TAUNtK»Lt 
0CPXNfK|s2.*fCPXNfKI<'CPXNtK*U) /VNK-0CPXNfK»LI 
DTAUNfKt32.«fTAUNtKt-TAUNtK«’LI)/Y»iK>0TAUNtK»Lt 
0BQNfKts2.*f8QN(KI>BQNfK»Lt t/VNK>OBQNfK*LI 
00 103 J=1»NSP 

0 ALPetJ«Kt 32 .«tAlPNf J,Kt«ALPNfJ*K»U I /YNK -OALPNfJ, K»U 
103 ODALPAf J«KI> 2 .«( 0 ALPNtJ.KI- 0 ALPM J.K^LI t/VNK-OOALPNfJ*K»LI 
101 CONTINUE 
1 C 2 CONTINUE 
NPTS»LAST 
RETURN 
ENO 

SUBROUTINE LPOINT tl .OPTPt 
CONHON/AB/EPP.EPQ.EPT 
COPNON/AC/ieOO.PIN 
CONNON/AL/GAR* GEM 
C 0 MH 0 N/eA/ALPtr, 55 ) .ENINF.MINF 
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ccfM0K/pe/sie,s2P.S3eT 

CCPHQKi/BC/G»NB.Pe.Q8«9HCq«THe«Ma.Xt1Ua.Ye 
rOHNOI«/CA/HQOTN{2.SSt .XNIBSI 
CCPN0K/CJ/CP(2.5Sl»CPl(7t ,CPX(56| 

QOHHOK/CK/MTttOUP <71 

CO«?10K/OP/YN<55| 

CCPNON/EC/CPIN*PQ 

C0Pt«QN/CF/£H<5SI«GAqt59l .PISSI , T h<G5I ,Y<5SI 

CCHHON/EC/EIN«PR«XtF 

C0P*t0K/€F/6At1INF,Hl(7),PINF 

COMMON/FE/OFU 

CCPHQK/CK/OELX 

CCPNQN/Hl/ALPHA.eETA 

C0PM0K/HP/ALPNI7,S9)«CPN<7,SSI •CPXNfSSI »EHN<5SI *GAHN<55 ltHN<7*55) * 
ILS«PN<S5I «QNtS5> «RHON<9qi *RN<S6I (THNCSSl »TN(55I «MN (SSI «XNUN <5SI 
CCPqOK/hK/CMC(2l .CP9(7|,CPXP(2I .OA(.OIF<7t,OALPa<7l,OOAI.PB(7ttOELS« 
AENP(2l,HP(7t.HC(7t,RP(2l •S3A(7) ,S3B( 7 1 ,S30(7I •MOOTS(7l ,XP(2I 
C0PH0K/CF/ALPB<7I tPHI (551 
CCPMON/PC/M ( 55 1 < X <55 1 
COPHON/PC/ JCMEM ♦NSP . T (55 > 

COPNQ^/QA/H(7,S5l «<)<95I .RH0(55I .XKU(55I 
COPNOK/CS/PHOP(2I «MQ0T(7.55I «HC0TC(7t «MP(2I •XHUP(21 
C0NN0N/SG/BQN(55I •OALPN(7,55t >r8QN(55l .0CPXN(5SI»00ALPN(7.55I • 
10TAUN(55l«TAUN(55l 

CCBM0N/SS/ALl«AL2.B'U.qC2.Cl.C2.CHl,CH2,081.0B2«a0l.002«0Tl.0T2,0V 

Al,OV2«Pxi,PX2.TAl,Ta2«THl,TH2,Vl,V2*Yl«Y2 

C0PHQK/YS/0VIS8.0VISC,IFS,MHP,VIS6.VISC 

COPHON/TU/SQ(55l«QA(.P(7,55l .09Q(5SI«0CPX(S5t*00ALP(7,55l*3TAU(95l • 
1TAU(551 

C0PM0K/TV/ALPP(7,2I •BFT,8QP(2I •CACHP(7«2I .0ALPP(7*2I *OaQP(2l « 
13CPXP(2l,C0ALPP(7»2i*0TAUP(2).0TCHP(2ltGAHP (2 1 .PP(2I *QP(2 1 « 
2TAUP(2I tTHP(2I.TP(2l,YO(2» 
ccpHON/uv/nii«iEPP.iPRess«iPP€su*isua 
COPMON/VT/OACH ( 7*551 *DTCH (551 tOtf ISO* vise 
COPMON/MW/MPTS * RE *X5P *X J 
OIMENSIOA 0ACMe<7| 

<PPFS£»0 

L*I 

K*I»1 

IF(OPTP.EC.C.I GO TO 2000 
K»I-1 

IF(IFS.NE.2I GO TO 5500 
IF((PPB/2I*2.NE*PPBI GO TO S5CQ 

EPtl(«XN2(ALPHA*8ETA*TMI l*XMU(I l*THN(I l«XNU»l(I )l 

eNlL«XP2(ALPHA*eETA,TH(K l*XHt(K l*THN4l l*xqUN(t II 

CO TO 0501 
8500 CONTINUE 

ENlRsXPl(ALPHAtBETA*TH(I)*XNU(II*THN(II *XNUN(III 
INl(.«XNl(AVPHA*8ETA*TH(tO *XNU(KI*THN(II *XNUN(II I 
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•SOI CONTIKUE 
2000 CONTIKUE 

ve>(VfKi*vmi/2. 

KIPbsQ 

8372 CONTINUE 

PATGs <ve-y ( I II / (V IKI -V < I I I 

TH8=TMfII *08TG*ITHIKI -TMIII I 
XNUft«XNU(II»PArG«(XMUIKI-XHU(ln 
EN2sXN2fALPMA.eETA*THS,XNUa.TriNftl.XNUNtIII 
ir lOPTP.KE.'" . I EP2=EP1L*PATG»IE'*1P-EH1LI 
»BT=te 

XR*XOF 

T«»rN«II-EN2*0£LX 
TESTT*IY6-YPTI /ITMXI -¥l II I 
IFtAnSITESrri.LT.O.Jll GO ' «37l 

<IP4sKlPG*l 

IF<KIP4.L£.2CI GO TO 837? 

WRITE (6*91911 
9191 FORNATIIHII 

WRITE(6.2C2?I 

2020 FCRHATISEP UNA9LE TO LOCATE Y LOCATION CF CHARACTERISTIC IN LPOINT 
II 

STOP 

8371 RATGs(Ve«Y(Itl/(V«0*Y(III 
51 TH9«Th(H*RATG*(TH|KI-TMmi 
qe*Q ( II »RAT6*(a(K 1-0(111 
PecP(II*RATG*(P(KI-P(III 
IT»T(H*R4TG»(T(KI-T(III 
TAUa«TAU(II*RATG«(TAU(KI-TAU(III 
8Qe«00(IT»RATG*(eO(KI-80IIII 
OCRX8>OCPK(I I »RATG*(OCPX (Kl -OCPX (III 
DTAU8«0T8U(II»RATG«(0TAU(KI-CTAU(III 
OBOOsOeO ( 1 1 YRAf G* ( 09Q (K I -09Q( 1 1 1 
OTCHP(ll«OTCH(IMRATG*(OTCN(KI-nTCH(III 
OTCHe>(OTCHP(ll«OTCH(III*.S 
CPXB«C.O 
HB«0.C 
CH2a«G. 

ru9R« f ^ I* 

CALL THERNO(TT«He«CPRI 
008060J«1« NSP 

ALPBUI«ALP(J* IMRATG*(ALP(J.KI-ALP(Jt III 

OALPB ( Jl «0ALP ( J * 1 1 YRA TC» (0 ALP ( J . K I -0 ALP ( J« 1 1 1 

OOALPBU) *OOALP ( J« 1 1 »RATG* (OOALP ( J*KI -OOALP ( J« 1 1 1 

CN20«CH2C*0ALPUt tl«CP(J. II 

CH2B« CN2 e«0 ALPO ( J I POPS IJ I 

CPXB* CPXB4ALP8 ( J I •CPB I J I 

HP»MB4ALP0(JI /MTHOLE I Jl 

NOOTB (JI«HQOT (J 1 1 lYRATC* (HOOT ( J* KI-HOOT (J« II I 


ORIGINAU PAQSjg 
OF POOR GUAIOT 
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DACNH«O.CH(J*I)»R«T6«fO*CHtJtKI-OACMU*Il» 
OACNB ( J) * (OaCHN^DACH U« I ») *.S 
OACHP«J«1I«OACHH 
V060 CONTIMIE 
Neat ./Me 
IMI*RO/Me 

6*Me«GP»e/ (CPxe*Re/cp in* 

ORalc/RB 

RH0eBPe«lie*6EN/TT 

EHB*Oe*EPlMF*SQ«1(C»P/6AHB»OR/TT» 

XHUBbZNUCEHBI 

IFfBEL.EB.O.I GO TO 0392 

Y PIllaY 8 

X PUlaX e 

Q PU)>0 8 

P PI1)«P 8 

T PflMTT 

N Pai«U 8 

P PUMP 8 

TH P(1I«TH 8 

EN PfllsEN 8 

8B P(i}«8Q 8 

RHO Pm«RHO 8 

XHU P(1I»XNU 8 

CPX PIllsCPX 8 

GAN PI1)«GAN 8 

TAB PIllaTAU 8 

080 P<1)«0BQ 8 

OTAUPliMOTAUB 

0CPXPUMCCPX8 

00 3939 J«1*NSP 

ALP P<J*1I*ALP BUI 

OAIP PtJ»tl*0ALP 8(J» 

3939 OOALPPtJ«i)*OOALP8(JI 
0S9B eONTtOUE 

XFtBETA.AE.OtOlGOTOLOrO 

TAUNIX)«TAUm 

eoNix)* 80 (Xi 

ocpxNfxi*ocpxai 

OTAUN<XI*BTAU(X) 

TN(XI*TIXIY0TCN(XI 

MN(XI*NIXI 

OIQN(X)«CBO(X) 

CPXNIXI«CPX(Xt 
CHBC>CH20 
OOAOriJal* N$P 
OAlPNIJf U»OALPIJ«Xt 
OBALPPIJ* LMOOALPUtll 
HCIJ»*HfJtXI 
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HOOTC(JI-MOOTNU,II 
<»e71 CONTINUE 
4070 CONTINUE 

IFfBETO.EC.O.Oi GO TO 4C72 
CHZC^O.O 
00 4073 J^ltNSP 
HC(JI=NN(J*II 
NOOTCUtsKOOTN(J«II 
4CT3 CHZC=CH2C«0ALPN(J,T)«CPN(J«II 
4072 CONTINUF 

visvise 

tfZsVISC 

ovi=ovise 

Ctf2»0WISC 

TAl^TAUe 

TA2=TAUN<II 

0T1*0TAUE 

0T2S0TAUMII 

OQl^eCD 

nQ23acMi) 

Ti*ye 

Y2sYN(II 

TH1*THB 

TH2sTHNm 

$lR«Sl(XJ.aEI 

CtsCPKO 

C2«CPXNUI 

oBisoeoe 

ce2soeQN(ii 

pxi«ocpxe 

PXZsOCPXMIl 

CH1sCh28 

CH2«CH2C 

IF(0EL.EC.a.Il«2 

1 S2BsS2(XJ«REI 

vi«viso 

0V1«0VISC 

TAlaTAUm 

OTisDTAUlII 

eQl«0Q<It 

Y1«YII» 

TH1«THII» 

S10«S1<XJ,RE» 

60 TO SO 

2 S2B«S2(XJ«9Et 

SO IF(OEL.EC.O.t3«4 

3 Vl-VISO 
OVl-OVISO 
ci>cpxm 
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eQl«8C(II 

oeiaoeom 

TAlsT»U(II 

THlaTHfll 

n*T«ii 

PXl«OCPHCI) 

CH1>CH20 
S20«S2CXJ«PEI 
It sserso. 

SSOTsC.O 

00A97SJ«1. NSP 

ALl»0«LPeUI 

ALZsOALPNtJ.Lt 

001«00ALPeiJI 

OOZoOOALPNIJ.LI 

vtotfise 

0V1»0VIS8 

eoisaoB 

TMlaT»«B 

V1*VB 

S38(Jt«S3fXJ«<tEI 
IP(OEL.eO.O.I5«8 
5 Vt*«ZSO 
OV1>OVISC 
ALlsOALPtJtll 
001*00ALPU*II 

8Ql«8Gm 

THi«rp(ii 

VlaVltl 

S30UI«S3(XJtRCl 
S30T>S30T«S30(JI/MTHOI.EUI 
8 SS8T«S3ei«S38O»/Nr,<(0lEU> 

U079 CONTINUE 

ZPIOEl.NE.O.I RETUPN 

Bi>pim 

OUN-XJ 

iiS«o 

VAX«Ve*VNCXI 

ZMVAX«LT.i.E-0€.ANO.XJ«NE.0.t 115*1 
XX«XJ 

ZPIX15. €0.1160*61 

60 XJ*0. 

61 82>P2a«Sie«SZB*S38TI 
XJ*XX 

XFUCHCH.Ea.il GO TO 7256 
83*0. 

60 TO 7257 

7256 TPl>(T(XI«OTCHm«TTI/2. 

e3*F3(TPlt0TCHe*TT*TN<XI*THa*rHN<Xt*OACHB«N8*NN(XI> 
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7257 0PTT«1. 

IFtOPTP.NC.O.I OPfT«-l. 

B«*F<i(eETAvOPTT«XHUe«THetXHUN(I)«THNIIII 

B2s(S2»e31*B<» 

IF{OPTP.NE.O..ANO*IPRES U.EO.O) GO TO 7<i4i4 
IFfOPTP.I>E*O..ANO.IPffES U.EQ.l) 60 TO 7462 
IFdBCO.CC.tl 60 TO 7444 
XFdPBESS.EQ.lt 60 TO 7482 
AXsl. 

lFtI15.EG.0t 60 TO 100 
AK«XJ»SINtXHUB» /SINtTHB-XHUBI 
lFIBETA.GT.0.0IAXstAX4XJ*StNtXPUNfll> 
l/SINtTMNtII«XNUMim*.5 
AXsl.*AK 
ICO CONTINUE 

PNtII«Pe*tTHB«AX»B2«tXNtII>XBI)/Bl 
GO TO 7445 
7402 CONTINUE 

KPRESS«KPBESS*1 
XFtKPPESS.LT.ei 60 TO 3232 
XERR«7482 

NRITEt6t3131l lERft.I.TNNIXI .PN( II •VNtll .TNB.PB.VB 
3131 F0RNATf2f5«6F13.5l 
STOP 

3232 THOUN«TMNtXI 
RIP4-C 

tHNtII*THB40PTT*Bl«tPNfXI*P6l»0PTT*B2*tXNtXI-XBI 
tFtA8StTHNtXI-TH0UNI.6T.l.E*04l 60 TO 0372 
VNai«Vm«.5«ITANITNtIM»TANITt<NtIIM*0ELX 
60 TO 7445 

7444 PNtI I •PB«OPTT« < THN < XI -THBI /Bl>62/ei* t XN t X I^XBI 

7445 CONTINUE 

XFIABSlPliaioPail.LE.EPPI PNtU«PtU 

0ELS*2*«(XNCXI«XtXII/CC0StTHfIII«C0SITNNtXIII 

TeRN2«BK0ltl«QtXt 

tFtBeTA.6T.O.OITCRN2><TERN2*RHON(XI 

1«QN(III«*S 

0T«i./TERN2 

QN«II«CSlC«0ELS*PNIXI4PtXII«0T*Q(X) 

XFtABSI0N1LI-QUII.tE.EPQI QNU>*QUI 
XFtBETA.SC.O.OICPXNtXUCPXtXI 

0TCNCN«0TCHIXI«tPNtXl-PtX}|*(QNIXl4QtXII/tCPXtXl4CPXNtXII*EXN*0T 

0T0XFP«S2C*0ELS«E XN*2 ./ tCPX I X I ♦CPXNI X II «0T 

TN t X I «T t X 1 40TCHEN40T0 XFF 

XFtABSITNttl-TtLII.LE.EPTI TNtU*TtLI 

CPXNtXI«0.0 

MNtXIaO.O 

CALL THSRNOtTNClIfHlfCPlI 
004000J«1« JISP 
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DALOIFI JMS30( Jl •0ELS/TERH2 

ALPN( J« IMALPI J« It ♦OALOIFtJI ♦OACHf J* It 

HNtJ.ItsHlUl 

CPN(JtII«CPl(JI 

MNat«NNIIt»ALPP«J* It/MTHOLE<JI 
CPXNCIIsCPXN(It»ALPNIJ« II»CPN(Jt It 
4080 CONTINUE 

MN(II«l./4Nm 

RNm«ftC/MN(II 

6ANN<It»CPXNm/(CPXN<II-PN(II/CPINI 

0R>1./PNIII 

RHON ( 1 1 sPM 1 1 •HN 1 It •6EN/TN ( It 

ENNf IlsOMII«EHINF«SQRriGAR/GANN(: )•OR/TNIItl 

XNUNdlsZPUIENNIItt 

PETURN 

ENO 

FUNCTION FERVfXl«X2*X3t 

C0NN0N/QR/0ELl«0EL2*RATltRAT2*SUN 

OERV»lkl*RATl-X2*fRATl-RAT2|-K3*RAT2l/SUN 

RETURN 

ENO 

SUOROUTIAE THSSSfTHSSI 
COPNON/AX/JSUeL*JSUnU 

CONN ON/EF /EN <55 1 «G AH < S5 1 , P ( 95 1 « TH (55 1 • Y < 55 1 

C0PH0N/CR/0£L1«0EL2*RAT1*RAT2«SUN 

CCNHON/RS/GPN8*PStTHS«THSL. THSU 

JSU8P«JSUEU*1 

JSU8H«JStPU*l 

CALL SMEAR(JSU8P*ASN1I 

CALL SNEAR(JSU8U«ASH2I 

CALL SHEAR(J8UBN«ASH3t 

EMS«0. 

0EL2>Y(JSL8PI-VUSUeUI 
OELl* V USUeU I • V ( JSU8N I 

SUN*0EL2«0EL1 

RAT1«0EL1/0EL2 

RAT2«CEL2/0EL1 

AY«OERY IASHl«ASH2*ASH3t 

C0STH«C0S(TH(JSU8UI t 

TERPl»-At*C08TN 

ENY>0ERV(EH(JSU8PltEH(JSUeUltEH(JSU8Ntl 

TANTH«TAN(TH(JSU8UII 

ENNN«ENV/C0STH«EHS*TANTH 

GPHbGANI JSU8UI Ppl jsueui •ENUSUPLI **z 

TERH2>2.«C08TH«EH(JSUeui*ENNN*PS 

CPMY«0ERV(GAN(JSUePI*P(JSUBPI*EH(JSU8PI«*2«6PN«GAH(JSUBHI*P (JSU8HI 
l»EN(JSU8HI»«2t 
6PHN«GPNV/C0STH>6PNS« TANTH 
THV>OERV (TH ( JSU8PI • TH ( JSU8UI • TH ( JSU8H 1 1 
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thhn»tmy/costh«tms*tanth 

TeRH3«6PHN*C0STH*THNN 

TE«HA*-6PMS*C0STH*C0STH*jemjSU9UI**2-l.»*THS 

TMSYsOE9Y(T“SU,THStTHSU 

TEftNSaGPM*SIN ITH ( JSUaU> t •tHSY 

THYPa DE9Y I TH ( JSueP»l I * TH ( JSU0P ) * TH ( JSUBU) > 

THYC> CEPY ( TH ( JSUBU) « ) H ( JSUBN) t TH C iSUBM* 1 1 1 

THNNPaTHYP/COSTH-THSU*TANTH 

THNNtaTHYt/COSTH-THSL*rAHTH 

THNYsCEBY(THNNP«THNN*THNHL) 

TEPH6»-GPH*C0STH*THMY 

0*CPM* lEP I JSU9U) •*2*C0STH*C0STH-1 . ) 

XNUH«TERPl»TE9M2*TE9M3^TERH4*T£9PS*TEfcM6 

THSSaXNUP/0 

PETURH 

END 

PUNCTIOH 2NU(EH) 

2HUaATAMl*0/S0RT(EH*EH-1.0ll 

return 

END 

SueROUTtNE THERHC(TI*H.Cei 
CCHMON/EC/CPIN«RO 
C0HNON/HK/RCO2* RHZOtMPUEL 
COHHON/TN/TIN 
OIHENSION MTH0L€(9) 

OIHPNSICK HtT),CPIT) 

OIPENSICN 0(9)*API9) 

MTHCLE(1)«1.098 

HTHOLE(2)al6.0 

HTH0CE(3)«1«*01€ 

HTH0LEI(»)«2.01G 

MTHOLEIS)*32.0 

MTHOLEI6)«12.008 

MTHCLE<T)«28*C14 

MTH0LEI8l«(><»«Ctl 

MTMOLEISIaMfUEL 

T«TI*T1N 

C1-RO/CP1N 

C2«C1/TIK 

no 10 j>i*R 

Hl«C2/NTNOLPIJ) 

H2*Cl/NTPOLE(J) 

CALL COEPP(J«T.A«a«C*0*F*FfG) 

Q1J»«T*U*T* C8* ,5*T*(C/3.*T*(0*.26*E*.2*T) )))♦£ 
Of JI>C(J)«H1 

APU»«A*T*«e*T*fC*T*<0*E*T))» 

APIJt>AP IJ)*H2 
10 CONTINUE 
H(1)*Q(1) 





129 


H«2I»CI2I 

H(3I«M«2Q JSI««C02*QI8) 

NC«I*Q(«I 

4ftSI*Q«SI 

(II6I«QC6I 

H<ne0(7l 

€eui««Mii 

CBC2I*«PC2I 

C8 ( J I *RH20MP < S I ♦RC02 • AP 1 9 1 

CeiAIMPIAI 

CS(5I««P«S» 

C8l6l«APiei 

cs(7i«iwm 

PETUBN 

END 

SUmOUTIAE feLlt::,TEHPY.X«V,MI 
OIMENStOK XUl.Vtll 
00 10 J5«1*N 
XF(TX»XIJSII «,4.10 

0 

TEPPY«V(^»*IV(JS>-VIJ6ll*irx«X IJ6tl/<X(jSI-XCJ91l 
GO TO 11 
9 TEIWY>YtJSI 
60 TO 11 
It CONTXKUE 
11 RETUON 
END 

FUNCTION XNl (M.PN««OETA.TA«XAtTC. XC> 

XNlaAlPHA«TAN«TA«XAI 

XF(0ETA«6T.Q.IiNl«XNl«0ETA*TAN(rC*XCI 

NETUNN 

ENO 

FUNCTION XN2(AL«e«TA«XA,TC*XCI 
KNO«Al*TANtTA*XA» 

XF(0.6T«0«I XHE>XH2*8*TAN«TC«XCI 

RCTUNN 

ENO 

FUNCTION XNKAtO.TOtTCI 
XN9«A«TANfTOI 

IF l••6T••••OI XNS«XN3^e*TANtTCl 

RETURN 

ENO 

SUMOUTXNE HOCUSITX.PltUltRNOlfALPHAfOXtLI 
CONNON/AC/ieOOfPIN 
CONNON/OC/IOCNEN 
CCNNON/CA/NOOTNfTtSSI tXM9SI 
CONNON/GE/RAO. R00« UIN tUXSINF 
60NN0N/Hl/0ALCN(7I tOTCMEN 
CONNON/Pr/I^LPNN (7 } .IF’jt L* PRES 
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CQHMON/QS/RMOPtSI •NOOf I7.SSI .110070(71 .MP<2I •XNUPCZt 
COnHON/TK/TIN 

OINENSIOit ASAVE(7I.MTN0LE(7I,M.PHA(7I 

MTKCLEdlel.OCA 

MTN 0 LE( 2 |s 16 . 

MTMOLE(31«lA.Oie 

MTNOL£(AI-2.0i6 

NTM0LE(S»-32.C 

MTPO(.E(€M17.00A 

NTNOL£(7M2«.01A 

TXXsTI 

PXK^Pl 

UXX«U1 

TERHxRN01*Ul 

TI=TI*TIII*.OOl 

P1»P1/FIA*P«»ES/2116. 

UlsUt«UIA 
0£LTATa%.E-7 
0ELTAX«U1«0ELTAT 
JERs I KT ( CX/OE L T A X I 
IF (JER.EQ.O) ^R«l 
DEL X«OX/FLOA T ( JER ) 

TSAWE*TI 
00 201 Jsi,r 
201 ASAVE(J)>AIPHA(JI 
OTsOELX/Cl 
PsPl 

0PS2116./09917. 

RM*P*CP/TI*..'.l 

00 10 JERRYsl.JER 

P*P1 

OUM»O.C 

00 96 J«l*7 

96 OUM«OOP«ASAVEUI/NTHOLEUI 
RHOIbRH/CUH 

IfdCCPEN.EQ.OI 

INRITE (6,2501 Tl*P.RHOI.ASAVE,OT,rN,ALPHN 
250 FORPAT(« POCUS FROM HnCUS*,10E11.3/17X,10E11.3/l 
call FOCUS(TI,P,RHOI,ASAVE,DT,TN) 

IPdCCPEN.EQ.OI 

IMRIfE (6,2501 TI,P,RH0I,ASAVE,0T.TN*ALPHN 
IFdCCHEN.EQ.OI 
IMRITE (6,2321 
232 FORMAT (//I 

IF(JERRY.AE.1I GO TO 100 
DO 110 J«l,7 

110 MOOT(J,LI«TERM*(ALPHN(JI-ASAVE(JII/OELX 
100 CONTIMUE 

IF(JERRY.EQ.JERIGO TO 10 
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TI*TN 

00 20 J«1.7 
20 aS«VE(JMALPHN(JI 
10 CONTIKUE 

0TCHEP8(Ta-TS«VEI *1000. /TIN 

00 %c 

OALCH ( Jl aaiPNNt J1 oALPHA ( Jl 
AO tfCOTNU.L»«TEPH*fALPHNfJI-ASAVEIJII/OELK 
TIaTXX 
PlsPXX 
U1*UXX 
PE TURK 
ENO 

SUmtOUTIAE COEPPtI.r.A .8 .C «0 «E *F ,G I 
IF«T-1000I1C.10.2C 

10 GO TO (15«1G«13«11«12«17*1A.1A.19).I 

11 A = 2.0A600A9E Od 

9 = A.1932116E-C3 

C --9.6119332E>G6 

Q = 9.Sl22662E>09 

E »-3.JC93A2lE-12 

F =-9.6725372E 02 

G S-1.A11785BE 30 

GO TO 40 

12 A 3 3.71S9946E 00 

a =*2.91672002-03 

C = O.SO37393E-06 

0 =-0.29987l6E-C9 

E » 2.7Q02100E-12 

F =-1.0S767C6E 03 

G = 3.9C0C7C4E OO 

GO TO 46 

13 A = 4.196S016E 00 

8 =-1.7244334E-03 

C = 9.6902316E-06 

0 =-4.5930044E-09 

E = 1.4233694E-12 

F =-3.0200770£ C4 

G =-6.0616246E-Ot 

GO TO 40 

14 A = 3.6916140E 00 

B =-1.3132552E-03 

C = 2.E5C31P0E-06 

0 *-9.7600341E-10 

E =-9.9772234E-14 

F =-1.062S336E 03 

G * 2.2074903E 00 

CO TO 4C 

15 A « 2.SC000Q0E 00 


0 _ _ 

,(^]nal paqbis 

POOR QUAIiro 
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8 c Q.O 

C = 0.0 

0 s o.Q 

E s 0.0 

r « 2.SO70%97C 04 

6 «*4.6Q01090E*01 

GO TO 40 

16 A = 3.0210094E 00 

8 s*2.175r249E*03 

C = 3.7S422Q3E-Q6 
0 S-2.9947200E-C9 

E » 9.0777547E*13 

F = 2.9137190E 04 

G s 2.E460076E CO 

GO TO 40 

17 A s 3.A234700E 00 

8 «*1.1187229E*03 

C = 1.246E819E-06 
0 **2.1C3S096E*10 

E «*5.2546591E*14 
F s 3.SaS2707E 03 

G « 5.62S3029E*ei 

GO TO 40 

18 AS2.1701 
8»1.037eE*02 
Cs*l*C7339£«05 
0*6.34592E«€9 
E«-1«E2B07£*12 
Fs«4.83$26£^04 
GslO *6644 

GO TO 40 

19 A»2*4912S 
8s7.64362£*03 
C«7.97754E*06 
0«*1.29S70E*08 
E«9.03078E*12 
F&-S421.ee 
G»C. 

GO TO 40 

20 GO TO C2Sf26t23«21«22f27«24»28«29)«I 

21 A s 3.0436897E 00 

8 8 6.1187110E-04 

C s*7«3993591E-09 

0 8-2.0331907E«11 

£ a 2.4593791E-15 

F a*8.8491002E 02 

G a«1.6481339£ 00 

GO TO 40 

22 A a 3.6976129E 00 
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e > 7*8U5603£«04 

C =• e.a3a6670E*07 

0 « 4»e49GlS9E-ll 

E s*3»346e2Q4£«lS 

P =-l*19279iee 83 

G « 3.74926S9C 00 

GO TO 40 

23 A = 2.6707932E CO 

0 = 3.0n7115€-03 

C =*a.5351570e-07 

0 • l.l790a53E^10 

€ «-€*197356aE-l5 

F **2.9aaa994€ 04 

6 - 6.aa3a39iE ca 

GO TO 40 

24 A « 2.aS4S7eiE 00 

8 s i.5976316E*03 

C **a.256G254E*C7 

D « i.i31Sa49E*iO 

E **7.e097070E*lS 

f »-a.9Gi7445E«02 

G a 6.3902879E OC 
GO TO 40 

25 A s 2.5CQ0000E CO 
8 > 0.0 

C ^ 0.0 
0 >0*0 

E « O.O 

r « 2.5470497E 04 
G «*4.e00l096E*01 

GO TO 4C 

26 A * ^.5372567£ GO 

8 «-l.a422190E-05 

C =*a.801792lE*C9 

0 « S.964362iE*12 

6 *-5.574360aE^16 

F s 2.9230007E 04 
G * 4«9467942E 00 
GO TO 40 

27 A s 2.aa95544E 00 
8 s 9*9a35061E«04 

C s*2.18799C4E-07 

0 « l«9a027a5E-il 

E «*3«a452940€-16 

F > 3*a0il792E 03 
G = 5.S597016E 00 
GO TO 40 

26 A&4.41293 
B«3.19229E*03 
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C«-1.2978C*06 

0«2.«14i7E-ta 

E«ol.67«3E*l«» 

F«> 4 . 99 «<iE« 0 « 

G»«. 7767C 
GO TO «0 
Z9 *s3.1E9<il 
B«1.0227GE'>02 
Ce-3.aS032E«06 
0»6.77198E>i0 
E«-G.5013?E-1% 

F*-S«%5.93 

GsC. 

9ETURN 

END 

FUNCTION 52 tXJ,FEI 
CONNON/EG/EIN*PR*Xt£ 

CONNOK/SS/ALl«AL 2 .SQl«BQ 2 ,Cl«C 2 *CHl,CH 2 ,Oei*OB 2 «Oai« 002 *OTl.OT 2 ,OV 

ai,OV>«PXl.PX2,TAl,TA2.THl.rH2*Vl*V2«Tl.Y2 

RPRsl./PR 

TERNl=Xl«Cl*OBl»PPR»V 2 *C 2 *OB 2 *RPR 
TEPN 2 =C 1 * 0 W 1 * 801 •RPR* C 2 * 0 W 2 *802 *RPR 
TERH 3 = < W 1 * 801 1 *V 2 * 8 Q 2 »CM 2 » *X l£ ♦PPR 
TFRN 4 =IW 1 * 8 Q 1 *PX 1 *W 2 * 9 Q 2 *PX 2 I*RPR 
TERN 5 =fWl*TAl»* 2 *V 2 *TA 2 ** 2 t *tlN 
IFtXJ.NE.C.I GO TO 10 
TERN 6 * 0 . 

CO TO 2 
10 YT=Vl*V 2 

lF<YT.LE.l.e- 10 l GO TO 20 

TERH 6 *Vl*Cl* 8 Ql*C 0 S<TMll/Yl*RPR*W 2 »C 2 *a 02 *C 0 S»TH 2 »/Y 2 »RPR 
GO TO 2 
20 CONTINUE 

TERN 63 V 1 * 01 * 081 •RPR*Y 2 *C 2 * 082 * RPR 
2 S 2 *«TERNl*TERN 2 *TERN 3 *TERH 4 *TERM 5 *TER« 6 l/RE*. 5 /eiN 
RETURN 
ENO 

FUNCTION S 3 (XJ«RFt 
C 0 NN 0 N/EG/EIN«ORtXL£ 

CONH 0 N/SS/ALl«AL 2 .BQl*BQ 2 «Cl«C 2 .CHl»CH 2 «DBl« 082 «O 0 ltO 02 tOTl« 0 T 2 ,OV 

Al*DV 2 .PXl.PX 2 ,TAl,TA 2 tTHltTH 2 ,Vl«V 2 «Yl.Y 2 

RPPal./PR 

TERNlaVl*C 01 *V 2 *CC 2 
TERN 2 «OV 1 «AL 1 *OV 2 «AL 2 
IFtXJ.NE.e.l GO TO to 
TERP 3 S 0 . 

GO TO 2 
10 YT*Y 1 *Y 2 

IFtYT.LE.l«E> 10 ) GO TO 20 
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TERns*COS (THl I vVi**!! /VI *COS t TH 2 I •V2«*L2/V2 
GO TO 2 

20 continue 

TEffNS«TERNt 

2 S3« (TERNl«TENH2*TERH3 » •XLE«RPR/RE • .5 
RETURN 
ENO 

FUNCTIOK Ft (HI 

CONNON/eO/GANa«P8«qB«RHC8«THe*Me«XNUa«ve 
GONNON/Mi/*tPNA .OETA 

C0NH0N/MH/«|.PN(r,SSI«CPN(7.SSI,CPXN(5SI.ENN(5SI,G«HN(551«HN(7.55l* 
1L*PN(55I«ON($5I,RHON(SOI«RN(5SI .THNfSSi •TN(SSt«MN(5Sl (XHUNISSI 
RP«1./P8 

F1«S INI XPV8I *C0S ( XHU8 I /6AM8*RP 
RPN^I./PNIHI 

IFI8ETA.6T.0.1Fl3|Fl*SIN(XHUN(NII*COS(XNUN(HII/GANN(NI*RPNI«.5 

RETURN 

ENO 

FUNCTION F2«N*Sil.S2t.SSll 
CONNON/80/6ANB.Pe«qB«RHOB,THO«Ma*XNUa,V8 
COHNON/OP/VN(S5l 
CONNON/NL/ALPN A .BETA 

C0NN0N/NP/ALPN(7.5SI .CPNITtBSI «CPXN(55I .ENNISS) .GAHNIS5 l*HN(7«S5l . 
IL.PNI55I .GNISSI .RHONISSI *RN(SSI t THNI5SI«TN(55l,MNI55l .XHUNISSI 
CONNOK/CA/N(7.SSI.Q(5SI«RHO(S5» tXKUCSSI 
CONNON/MV/NPTS* RE • XBP« X J 
IFIXJ.EQ.0.0ITERN1«0.0 
IFIXJ.N€.0.ITERNlsSIN(TH8l/V8 

IF(XJ*NE.O».AN0.8ETA.GT.O.ITERM1«.$«|TERN1*SIN(THN(H| I/YN(MII 

0S*l./O8**2 

TERH2«S11/RH0B*QS 

SQ»t./GN(NI**2 

IF(8ETA.GT.O.)TERH2s.5«(TERN2«S11/RMON(NI«SQI 

Pi«t«/P0 

TERNS*S21« (GAHB'l • I /G ANa«Pl/QB 
P2«t</PN(NI 

IFI8ETA.GT.0.ITERN3«.S*(TFRH3»S21«(6AHNINI>1«I/6ANN(HI*P2/QN(HI| 

RQ*t*/Q8 

TERNA«S31«NB/RHOB«RQ 

Q0«t./QN(N| 

IF (BETA.CT.e • I TERRA*. 5* ( TERNA^S 31 •HNI HI /RHON(N) «Q0I 

F2«(TERM1«TERN2 TERN7-TERHAI 

RETURN 

ENO 

FUNCTION FA(e.0PT.XHUl*lHl.XNU2«TH2l 
FA«SINIXNU1I/C0S(TH1*0PT»XNU1I 

IF (B.GT.O . IFA* (FA»SIN (XNU2I /COS (TH2«0PT*XNU2I t *.5 

RETURN 

ENO 
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SU8R0UT ll€ HERMAN < VHt OT • A • Y «C1 1 BR tCC t SC ALE I 

OINENSION PC10«10l»SNALe(lQI«QC10l«AC10«10l»VI7l«YNC7l»CI(4l«fINKI 
lAI 

TIN1»OT/2«0 

TINE«OT 

TOsTIMl^^E 

T2«IOT»«3^T1H1»TOI/3#0 

Y3»T0*.5 

f4»TINlMO/3.0 

K*l 

00 19 

DO 10 J«l«4 
PII(«J|s«A«l«JI»T3 

10 P<K^l%J|s-All*JI*Tl 
19 K»K^2 

K«1 
00 20 

DO 11 JsltA 
P(K»J^A|s*Aa»JI«fT4l 

11 PIK41^Jt4l**AIItJI*<T2l 
?Q K«K«2 

J«1 

00 12 Isl«B«2 
S«1»/SCALE 

PfI«JI«PatJI^TlMl»S 
PCIf J«4»8PI1*J«I||4T0»S 
K«I»1 

P<KtJ|8PIK«JI#OIP2*riHlMS 

PfK»J^4>«PfK»J^4>42.*Tl*S 

12 CONTINUE 

00 13 X«1»0 

13 Q(IM0*0 
fINKIll8YCl) 
flNK<2l«YC2) 
fINKI3|8VC6l 
FXNKI4|8Y<3I 
K«1 

00 IS UUk 
DO 14J»lf4 

14 QHO«C(KI^A(X«J)»fXNKUI«CTlH2-TXMll 
QlKtll^QCK) 

15 K*Kf2 

00 IS I«l»4 
J«2»I 

OU^lMOU-ll^CXai^lTXHZ-TXMH 

16 QIJI«QIJI»CXm^CTXM2-TIHll 
00 202 I«l*e 
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QCII«0CII/i.0E*5 

00 

2€2 PCI»JI«PCI«JI/i.O€-S 
CALL CLEHId«SHAL6«<»tQ» 

CALL SCL?(SV1ALe«0T«CC«BetV«YNI 

RETURN 

END 

SUBROUTINE CLENCH«X«8tDI 
OIHENSION ATIiOtliUXIiOl 
OINENSION BClOtlOltOClOl 

00 12 X»i«H 

12 xm«c*o 
00 aoc i«i«N 

EOO ATCItPlMCm 
00 201 I«1»N 
CO 201 J»1#H 
201 ATIItJ)«eiI«JI 
00 32 N»1«H 
OsATIN«NI 
IT=0 

00 9 IsNtH 

IPlABSUTft«NII*AeS<OII 9«9t0 
a 0«ATCI«M 
IT«I 

9 CONTINUE 

IFaT-N>7t/f7C 
ro 00 71 JsNfHl 
TENP«AT(N»JI 
ATCN»J»sAT(IT«JI 
71 ATflT«JI«T£HP 
7 00 10 X«ltHl 
10 ATCNtIlsATCNfll/G 
IMN*NI90f9CflO 
IS NlsN^l 

00 so X*N1«N 

0»ATCI«NI 

00 so J»N«R1 

30 ATCI«JI«ATIX«J)*ATIN«JI»0 
32 CONTXNUE 
90 XCN)«ATCHtNfl) 

00 69 N*2tN 
NRsH^l«N 
0»AT(NR«N«1I 
00 60 X«NR«N 
60 0«0-ATCNRtI)«xm 
69 XCNRI«0/ATCNRtNR) 

RETURN 

END 
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SUmOUTIItC SOLT (SHALB «OT *CC •BB. V , VN) 

OINCNSIOM SHAL6(10)«V(7I,VN<7I 

TIHEaOT 

TNK»TI«€**2 

YN I 1 I eV I 1 I aSNALB ( 1 I *T XNE» Sn AL B I 5 I X 
VN(2l*VI2»«SNAL8l2»*TINr*SNAL8(6l«TNX 
VN(6l»V(e)«SHALB(3|*TINE«SHALB(7l«TNX 
VN (3 1 » V (3 1 »SHAI.B ( Al « T IHEaSMALB < b I «TNX 
VNCAIaCC-lYNtlMYNIBI l•.S•VN(3l 
VNtStsBB- (YNI2I*VNI61«VN(3II*.5 
RETURN 
END 

FUNCTION SKXJ.REt 

C0HH0N/SS/AU*AL2*BQl«BQ2.Cl*C2*CHl*CH2*0ei,0e2«0ait002.0Tl.DT2,0V 
Al,0V2*PXl«PX2«TAt«rA2.THi,TH2*Vl*«2*Vl«V2 
TERNlaVl*0Tl»V2*0T2 
TCRN2«0V1*TA1«0V2*TA2 
IFtXJ.NE.O.I GO TO to 
TERN3B8. 

GO TO 2 
10 YT»Y1*Y2 

IFIVT.LE.l.E-lOl GO TO 20 

TCRH3«C0SITN1I •Vl*TAl/VlAC0S(TH2t •¥2*TA2/V2 
GO TO 2 
20 CONTINUE 
TERN3«TERNi 

2 Sl*(TERNlbTCRN2«TERN3l/RE*.S 
RETURN 
END 

SUBROUTINE PUNCH 
CONNON/AC/ IBOO* F IN 
CONNON/eA/ALP(7,SSI .ENINFfWINP 
C0NN0N/0e/BCTB(AI *IS(AI 

C0NN0N/EP/EN(SSI,GANf9SI*PIS9l«TH(5$l«V(55t 

C0NN0N/E6/E1N*PR*XLE 

C0NH0N/FH/XK1«XK3*XP0T 

CONNON/HU/KOUNT«LL«NPT 

COMNON/PC/ALPHN 171* IFUEL t PRES 

CONNON/PO/JCNEN*NSP«T(99t 

CONN ON/TH/TIN 

CONNON/MV/NPTS«RE*XBP.XJ 

CONNON/XV/APRS*APUS*OECTAV«EBOOStI80DS*INTACT«IPRS*lPUS*ITVP, 
lJ800StNNAXtRNEAT*XK2*XK4tVB0T«VTP 
CONNON/VX/ABOOS*BPRESS»CPRESS 
CONNON/YI/BpRESU*CHENFC«CPRESUtENSU8tRTH,XSTEP 
CONNON/ZY/A800»6BOO«CSOO«EBODfFflOO*GBOO*IAVE«IPUNCH.JBOO*KK<KK 
RCNIND 7 
too FORNATIIEIS) 

101 FORNATIOE10.3I 
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200 fOOHAmStSX.7eiO.3t 
1C2 fORHAT(7ei0.3tflO.St 

103 FCRHATISflO.St 

104 F0RHAT(7E11.4I 

MRiretr.iooi kkkkk ,ll 

MRITEt7.200t IPUNCH.XSreP 

INTACToO 

ISHOCMO 

DO 1111 l>lt4 

1111 IFdStll.NE.Ot 1SH0CK«1 

MftITE (7. loot NPTS.NPTt ITVPt ISHOCK.NHAX tKOUNT 
MRITEI7.100t JCHEHtlAVE. INTACT 

MRlTE(7.1C2IXJt EHSUA.RTH. OELTAVtVeOTtTTP,CHENfCtXBP 

RO>RE/RTH 

HRITE(7tl01IRQtPRtXLEtEHINFtTIN.NINFtPRES 

NRlTEI7tl01IXPOTtXKltXK2tXK3tXK4 

MPITE (7,200 1 t8ooStA80os,aeoo.ceoo 
MRlTE(7t200IJBOOStEBOOS,F800.6800 
MRITE(7t 2001 tPRS. APRS .8PRESS.CPRESS 

MRtTE(7,200ttPUStAPUS ,8PRESU.CPRESU 

IFItSHOCK.EQ.OI 60 TO 5 
HRITE(7tie0MIS( It .1*1.41 
MRITEIT.lOlt (BET8ai.I*l,4l 
S CONTINUE 

to 13 Isl.NPTS 

A *p(ii/p:n 

ALP7*ALP(7.II-ALP(4.1t«(l.-RNEAT| 

ALP4*ALP(4,It/RHEAT 

MRITE(7.103tY(II.A ,TH(II ,EH(II .Till 

WRITE (7.1041 ALP (1. It .ALP(2. 1) . ALP (3. It . ALP4. ALP(S. It .ALP (6.1 1 . ALP 
17 

10 CONTINUE 
RENINC 7 
END 

SU8R0UTINE INDATA 

C0NH0N/AC/I80C.PIN 

CONNON/AL/GAR.GEN 

C0NH0N/8A/ALP(7,S5I,EHINF,MINF 

COHHON/eO/XHASS (SSt 

CCNN0N/CJ/CP(7,55I.CPl(7t,CPX(5SI 

C0NH0N/CK7MTH0LE (7) 

C0HH0N/0B/BET8(4I « IS( 41 
COPNON/EO/CPIN.RO 

C0NH0N/EF/EN(55t.GAN(SSI.P(5SI.TH(SSt .Y(5SI 

COHHON/EG/EIN.PR. XLE 

COHNON/EP/6AHINF.H1(7).RINF 

C0HN0N/FH/XXl.XK3,KP0r 

CONNON/GE/RAD.ROO.UIN.VISINF 

COHHON/GF/OELY.CVISA.KOUNTO.VISA 
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CONHON/H J/KOUNT f LL t NPT 

COMNON/HK /RC02 1 RH20 1 MFUEL 

C0MN0N/0R/THBP*VeP«V8PN 

C0HN0N/PC/MISSI*XI5SI 

C0PN0MPC/ALPPNI7ttIfUCL*PRES 

C0NN0N/PC/JCHEN»NSP*T(55} 

C0PN0N/CA/H(7«$$» »ll(95l «RH0I55) .X»< JI55) 

C0MM0|C/RC/R(99I 

CORNON/TM/TIN 

C0NN0MUV/mi*IERR.tPRESStIPRESU*ISU6 

C0MN0N/M«/NPTStREfX8P*XJ 

CCNHON/RX/APRESStAPRESU 

C0RH0W/XY/APRStAPUSt0£UTAy,EiOOS*l800S» INTACT, IPRS. IPOS. ITYP, 

1 jeoo S , HP A X , RHE A T , XK2 , XK A , VBOT , V T P 
CONNOAyYX/AnOOS.ePRESS.CPRESS 
C0NN0N/Y2/BPRESU,CHEHFC,CPRESU,EH$Ue*RTH,XSTEP 
COPHON/ZY/AOOO,eeOO,CeOO,EBOO,F90C*6ROO,IAVE,IPUNCH,JBOO*<KKKK 
IIN*S 
ISU8«C 
XPPsO, 

Yepaioooo. 

TH8P<iC. 

RAO>0. 

IFUEI.«1 
MFUEL«2.016 
100 FORNAmOISt 
ICl FORNATI8E10.0I 
4CA FORNATtmi.AI 
2C0 FORNATfIS,SX,7ElO.O) 

REAOlIXN.lOCt KKKKK.LL 

REAO(II»>,2001 IPUNCHtXSTEP 

READ (XIN, 100) NPTSfNPT,XTVP,XSNOCK,HNAX,KOUNT 

REAO(XIN,100) JCHEN.XAVE. INTACT 

XFfKOUNT.LT.l) KOUNTsO 

KOUNTC»KOINT 

MRXTE (0,111) KKKKK.LL 

111 FORMATIOHIKKKKK «XS,9X,4HU >13/) 

NRXTE (6,112) XPUNCH.XSTEP 

112 FCRNAT(9H XPUNCH sX2,SX,7HX$TEP *E10,3/) 

NRXTE (6,113) NPTS,NPT,XTYP,XSH0CK,HNAX 

113 F0RHAT(7P NPTS bX3,9X,5HNPT >X2,9X,6HXTVP «X2,9X, OHXSHOCK >I2,9x, 6 
IHMNAX aX3/) 

NRXTE(6,11A) JCHEH.XAVE, XNTACT 

llA »URHAT(OH JCHEN «X2 .9X,6HXAVE «I2,9X, 
lOHXNTACT «X2/) 

XF(XTYP.KE.2.AN0.IT»P.NE.A) CO TO 12 
XF(XNTACT.EQ.0.ANC,XSHOCK,EQ,O) CO TO 12 
NRXTE(6,9191) 

9191 FORNAT(lHl) 
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NRXTG(6«102I 

i02 F0RNATt91H TYPE 2 OR TYPE h PLOWS HAY NOT START WITH SHOCKS OR HAV 
6E SHOCKS CONXNG OPP SPLXTTER PLATES/A3H RECHECK INPUTS AND SU8NIT 
INXTH PROPER TYPEI 
STOP 

12 CONTXNUE 

10% REAOCXXNtlOll XJ, ENSUe«RTH*OELTAY,YeOT»YTP*CHEHFC,xeP 

REAOf IINtlOl) RE.PRtXLE.ENINP.TINtWINP.PRES 
REAOaXN.lQll XP0T«XKl*XK2*XK3«XK<i 
REAO<nA*200l ieOCtA800.aeOO.C800 
READ t UK. 200 1 JBOO.C800.P800.C800 
READ! IIN. 2001 IPRESS. APRESS.8PRESS.CPRESS 
READ < nw,200l IPRESU. APRESU.8PRESU.CPRESU 
iPtxep.LT.o.i xapso. 

J«XU».S 

XJ«J 

WRITE <6.1191 XJ.ENSUe.RTH.OELTAV.VeOT.VTP.CHENFC 

115 FCRHATI9H XJ *E10.3.2X.7HENSU8 -E10.3.2X.5HRTH SE1G.3.2X.SH0CLTAV 
1«E10.3.2X.6HYBOT >£10 .3.2X.5HVTP «E10 . 3.2X.SHCHEHFC >ElC.3/t 

WRITE <6.1161 RE.PR.XLE.EHINF.TIN.WINF.PRES 

116 FCCHATlSli RE >ElO .3.2X.4HPR «E10.3.2X.5HXLE «E10.3.2X.7HEHINF >E1C 
1.3.2X.SHTIN SE10.3.2X.6HWINF >El0.3.2X.6HPRES >E10.3/I 

WRITE<6.1171 XP0T.X<1.XK2.XK3.XK4 

117 F0RNAT<7H XPOT «E10.3.2X.5HXK1 ^E10.3.2X«SHXK2 >E10.3.2X. 5HXK3 «E1 
10.3.2X.5HKK% «E10.3/1 

WRITE<6«110I 1800. A800. 8800 .C800 

110 F0RNAT<7H 1800 >I2«2X,6HA800 >ElQ .3.2X.6H8800 >ElO .3.2X.6HC800 «E1 
10.3/1 

WRITE <6.1191 je00.E80O.F80(l.6eOO 

119 F0RHAT<7(4 JBOO «I2.2X.6HE80D «E10.3.2X*6NF800 «E10.3.2K.6HG800 >E1 
10.3/1 

WRITE <6.1201 'PRESS.APRESS.epRESS.CPRESS 

120 F0RHAT<9P IPhuSS >I2.2X.0HAPRESS >£10 .3.2X.0H8PRESS «E10. 3.2X .OHCP 
IRESS >E10.3/1 

WRITE <6.1211 IPRESU.APRESU.8PRESU.CPRESU 

121 F0RNAT<9W IPRESU »I2. 2X.0HAPRESU >E10. 3.2X. OHBPRESU «E10.3.2X,SHCP 
IRESU «E10.3/I 

All ie00S>I800 
AB00S«A800 
XPRSb ipress 
APRSbAPRESS 
jeoosBueoo 
EBOOSbEBOD 
IPUSb IPRESU 
APUSbAPRESU 

IF<ISHOCK.EQ.01 GO TO 5 
REAO<IIW.100I <IS<I1*I«1.%1 
REAO<IIK.101I <BETe<IltI«l*41 
NRITE<6.1201 <IS<II«I*-i.%l 
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128 FORHAtfSt* ISIll «I3*2X*7HtS(2» aI3.2X«7HlS(3l •Z3*2X«7HISUI > 13/1 
MRirr<6*i29t tecTeai«i«i* 4 i 

129 FCR.rATUOH eCTeci) >E10.3*2X«9H8ETBI2I s£10.3«2X«9HeETe<SI «E10.3* 
12X.9HaETeUil sEl0.3/t 

S CONTINUE 

00 10 IsltNPTS 

REAoiiiNtioi) Yfii»pm.THii)«ENm.rni 
REAO(IIN.404I <ALPfJ*It«Jel,NSP» 

10 CONTINUE 

1P<1TVP.EQ.1IG0 TC 4201 

IBOOsO 

JBCO«C 

IPIITVP.E0.3I jeOO^JBOOS 
IF(ITVP.EC.4I IBOOsIBOOS 
4201 PPfATal. 

9«20=1. 

RCC2aO* 

4204 HTH0LEC4)aNFUEL 
RE«RE*RTH 
lERftaC 

CALL C0EFF|9,TIN,AZ,BZ.CZ«0Z.NZ«FZ*GZI 

CPIN»<AZ*eZ*TIN*CZ*TIN**2*OZ*TIN**3*HZ*TlN**4l*RO/NTNOLE»5» 

CALL C0EFF(7,TlN«»Z*PZtCZ*0Z*HZ*FZ«GZI 

CPIIafAZ*BZ*TIN*CZ*TIN**2*0Z*TIN**3*HZ»nN**4l*R0/MTH0LFC7l 

CPIKa .232*CPIN* . Z68*CPI I 

RINF*RO/HINF 

GAMINF>1./(1.-RINF/CPINI 

RINFFaRCC/MINF 

UINaEFINF*SQRTI6ANINF*RINFF*TiNl 

RF«1./RINFF 

RNOINF«PRES*RF/TIN 

WISIMF«R»<CINF*UIN*RTH/RF 

GAR«GANINF*RINF 

CEW*GAHIKF*EMINF**2/NINF 

ElNa (GAF.NF-l.l • ENINF**2 

FHS«1./ENINF*»2 

PlN»l,/6ANIF?®aEHS 

WRITE (6«6896l 

6B98 F0RNAT<////48X«31HP ROGRAN VIS'CHAR //GCX.rHN I T N 
1//42X.43HE HBEOOEO SUBSONIC FLO W//S3X,21HS K 0 
1 C K NAVE S//33X.63HA NO FINITE RATE H2>AI 
IR CHEMISTRY! 

IFfXJ.EQ.O.I WRITE(6tS610l 
IFtXJ.NE.O.I WRITE(6t 96111 
IFIJCHEN.EO.OI HRITEI6fS612! 

IFfJCHEH.EQ.il WRITEf6*6613l 
5610 FORNAT(///10X.31HTYPE OF FLOW IS TWO OIHENSIONALI 

9611 FORNATI///10X.20NTYPE OF FLOW IS AXISYHNETRICI 

9612 FORNAT(10X«19HCHEHISTRY IS FROZEN! 


OF POOE <»JAUi:w 
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5613 F0R16T«1«X*24HCHEHISTRV IS FINITE RATEI 
MRITET6t560C» RTH 
56Dt FORNATI 

tlOX,?8HJET OR NOZZLE RADIUS (RTHI « E13.S.6H FT. I 
MRITEI6.5601I ENINF.UIN.TIN.PRES*RH0INF.6ANINF,MINF*RE.PR«XLE 

5601 FORNATI///20X.20NREFERENCE CONOITIONS/20X.20M— — — — — — / 

110X.16HNACH NO. (ENINFI = E13.5/10X»16MVELOCITV (UINI s E13.5* 

1 7H FT/SEC/10X.19HTENPERATURE (TINI »ri3.5*10H DEGREES K/IOX.IFHPR 
lESSURE IPRESI «El3.5*9 H L8/FT«*2/10X«l6MOENSirT (RHOINFI >E13.5* 
112M SLUGS/FT««3/10X,37 

IHFROZEM SPECIFIC HEAT RATIO IGAHINFI «E13.5/1GX,25HN0L£CULAR NEIGH 
IT TNINFI «E13.5/10X*22HREYNOLOS NUN«ER CFEI «E13.5/13X.21HPRAN0TL 
INUHDEF fPR» sE13.S/l0X.2OHLENlS NUN6ER CXLEI sEl3.5l 
MRITE(6*56«2I 

56i2 FORNATC///20X«15HOUTPUT HEADINGS/20X«15N~ 

IlOX* 5HX - X/RTH/lOX* 5HV - V/RTH/10X.16HQ - VELOCITY /UIN/ldX. 
119HT > TEHPERATURE/TIN/10X.17HP - PRESSURE/PRES/10X«3CHTN > FLON D 
1EFLECTI06 IRAOIAKSI /10X,16NEH - NACH NUNDER/lOX.EdHRHO - DENSITY/ 
1RHOINF/10X,19HGAH • SPECIFIC HEAT 
1 /10X.33HXNASS • NON*OINENSIONAL NASS FLON 
1 /10X.23HPHI - equivalence RATIO/IOX. 

120HN • HCLECULAR NEIGHT//10K16HHASS FRACTIONS/15X10HALPHI - H/15X 
110HALPT2I - 0/15X*12HALR«3I - H20/15X»11HALP(4I - H2/15X.11HALPC5) 

I - 02/15X.11NALPI6I - 0H/15X«llHALPf 71 - N2I 
613 DO 1776 Isl.NPTS 

XfllsXDP 

pm«pai*PiN 

ALPf6»II»RNEAT*ALPI6*II 

ALPI7.Il«fl.-R*<EATI*ALPC6.II*ALP(7*II 

DO 7S0 J»1.*(SP 

766 IFIALP(J*1»«LT.1.1Eo10I ALPIU«II>1.1E-10 

ALP(7.1l«l.>(ALP(l*IMALP<2»I»*ALP(3»II*ALP(6«I»«ALP(5«II»ALPf6«II 

II 

1776 CONTIMiE 

00 0083 I«1«NPTS 

CALL THERNOaai.Hl.CPlI 

CPXIIIsO.O 

NfllsO.Q 

00177E J«l« NSP 

CPU.IMCPKJI 

HIJ«II«H1 (Jl 

CPX(II«CPX(I)»ALP(J*II*CPU«1I 
177e N(IIsMIII«ALP(J*II/NTNOLE(JI 
Nm«l./Nfll 

RH0IIlB6EM*NfII*Pa>/TIII 

RIII«RG/H(It 

6ANm«CPX(ll/ICPX<tl-Rm/CPINI 

0N«1./EHINF 

ORbI./RCII 
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QCII«eMII«0H/SQRTI6«lt/6«NCtl*0f>/T<IM 

IFtEHdl.CT.i.l 

IXNUa >b«T«NU ./SORT If HCl I ••2-1. 1 1 
86S3 CONTtnUE 

tFltRTACT.EQ.il CALL CONL 

XJl«i.»XJ 

KNASSUIsO. 

00 178S I*1.NPTS 
tFtt.EQ.ll 60 TO iros 
XJl*t.«XJ 

TFUN««TCtl*lt.-XJ«Yttl«XJI> V«t>ll*(l.*<LI»Ya-il«Xjn/XJl 
RQAV*IRHOftl«Q(tt«COS(THCtll« RN0lt*ll*0f I*l>*C0StTH(t*llll/2. 
XNASS<I|sXHASSIt-ll»RQAV*VFUM 
1705 COHTtNUE 

OELT« CY8R>Y 1 1 1 1 /FLOAT tHPTS-1 1 

RETURN 

ENO 

FUNCrtON F3CTPl*0TC«ri*TC.THl*THC.0A«Rl.HC» 

CONNON/CK/MTROLE f 71 

C0NN0N/6X/0ELX 

COHNON/HL/ALPHA.OETA 

OlNENStCN 0AI7I 

NSP«7 

A« ALPHA 

ObOETA 

TERNl«0TC/CtA-0l«TPl»0«(Tl»rCI I 
TERH2«0. 

00 10 J«1*NSP 

TERN2*TERN7« OAUI/miHOLE(JI 
10 COHTtNUE 

TERNS sA*Hl*e*NC 

TERN2«TERN2*TERN5 

TERN3«A*GCS<TH1I«0*C0SCTHCI 

F3**<TERPl»TERHtl •TERH3/0ELX 

RETURN 

ENO 

SUOROUTtNE POCUS < T ttPRESS t« RHOt * ALPH I • Of « TNI 
C0NN0N/P0/ALPHM7 1 « tFUEL* PRES 

OtNENStON ALPHII7l*AOUP«10l*Ct(10l«V<7l*YNf}|«ALPHAf7l 

OtNENSm T0t7I.TlC7l ,0171 ,CI7I .0(71 •E<7l«6t7l«Z(7l 

T0I1I«6.0 

T0I2I«6.0 

T0I3I«0.5 

T0«ila0.5 

T0(5I«0.5 

TOfOI«O.S 

T0C7|s0.S 

Tl<ll>6.0 

T1I2I«6.0 
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Tt(3l«3.02$9 
Ti 14 1 Ml. 0960 
T1(SI»2.9Z02 
U(6l«3.e392 
Utn«2.400e 
9tll«39*709S 
B«2I»2.56?4 
e(3l«3.«961 
9C4I *27.4023 
8(91*1.7771 
8(61*3.3496 
8(71*2.9043 
CdMO.O 
C(2l*C.O 
C(SI».S486 
C (41 *1.5999 
C (6) *.1599 
C(6>*.1619 
C(7I*.1931 
0 ( 11 * 0.0 
0(2)*C.O 
0(31 **31. 7050 
0(4I*>34.5268 
0(5t*«1.8504 
0(61*1.3139 
0(7) *-1.4976 
E(l)«t.C 
E(2I*0.0 
E(3)«E.3657 
E(4) *30.9104 
E(5)*2.5«21 
E(6)«4.3679 
E (7) *2.6093 
6(11*404.5564 
6(21*29.1774 
6(3>«-26.9024 
6(4)s-8.006 
6(5) *-.522 
6(6) *3.4213 
6<7)»*.5961 
Z(l>**063 
Z(2)«1.0 
Z(3)*t.l3 
Z(4)«.126 
7 - )*2.0 
Z(6)*1.063 
Z(7)«1.7S 
8SSSS*P8ESSZ 
KASE*1PUEL 


POOR aiAun 
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IFIKASl.CQ.2l PKESSl«PRCSSI«.3S 

IFtKASC.CC.SI PQCSSI«.2«PQFSSI 

RMOI«RHQI*PRCSSI/PSSSS 

KTEST=0 

CLQsl.O 

OLTIsQ.O 

CPSs.aOt 

T IHC Q *1 . 3Sr 2SE*S*ELO 

OrsOT/TIPCQ 

PObPACSS1«1.01329C6 

RNOO«PC*1 .92<iA6Sc* 19 

RMO IsRHCI * .S15«/RH0a 

PRESSIsl.fi 

TsTI 

HIsQ.O 

00 65 Isl.7 

IFfT'TlClM 62.61.61 

61 HIslOCII«CIII*TI*AtPHtm.MI 
CO TO 65 

62 tftT'TOtXII 63.63.66 

63 HIc(GUl«eiII*TI«A(.PMI(IMHI 
60 TO 65 

66 HX«<&(II«eai«T«Cai*<ToT0mi««2t«ALPHl(l|6Hl 

65 CONTINUE 

92 CONTINUE 
JUJ s 25 
JJsQ 
T « TI 
TSAVCsT 
KOUNTsfi 
RMOsRNOI 
OElTAsOlTI 
6 AHNAs 0T«0ECTA41. 

PRESS-PRESSI 

H*HI 

SUHVsfi. 

00 11 Isl.7 
AlPNA(IlsALP9I(II 
vmaRN0«ALPHAfI»/2III 
VNdlsfi.O 

11 SUNY«SUHV«V(I) 

0(mi««.6703lC-7*RH09*CLO 

OUN28CUN1«RHOO/16. 

XF(ALPNA(3I.GT.l.E-lfilG0 TO 6 
IFtALPNA<6t.GT.l.E-10l 60 TO 6 
If (ALPNM5I.GT.1.E-10I GO TO 30 
IF(ALPHA(2I.GT.).C*10I 60 TO 30 
F5*fl.05El7*EXP ’«25./TI»»<0UN1«EXP (•29./TI/T) 
05*1 . E16*DUH1*RS V 0/16 . 
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•ll«-CFS*.S»2.*8f •mi l•SUNV 

CCl«Bf «V (11 ••2»SU«f 

COGAMIAviYf 2I*Y(11 *.51 

Cl»F5«CC«SimY«CCi 

«li«0ELTF«811 

OUIt>Cl/«li 

YNIil»-OQI«»IYIil»IHINI*EKP («tl*OTI 
XFIYNCII.IT.E.OI YNtll«8.0 
YNHI«GC-Ylllil*.5 
60 TO 99 

30 tFUtFM«(9I.CT.1.6<*iOI GO TO 6 
IFMLmiCll.61.1.E*ttl GO TO G 

F8«l5.0Ei€«EXF C-3t«3/TII«<0IINl«EXP {•SO.S/Tl/TI 

MbG • EtG*OUMl *0000/16 • 

eil«- (F«* •00*Y <i 1 1 *9009 

cci«eo* Y I 1 I *T « 1 I *sunY 

88«G«PI«6«(YI2I*VI1I«.SI 

C1«F 8*e8«SUNY«CCl 

811«0ELT8*eil 

OUHaC 1/811 

YN<ll»-0UF*(Yfll*009l«EXP «811«0TI 
YNt2|s8eoYNfll*.5 
tF(YNI2I.LT.0«01 YN(2I*0.0 
YN(51*88-YN(21*.5 
GO TO 99 
6 CONTINUE 
K0UNT«1 

IFIKASE.EQ.2I T«l./fl.l087/T>. 096971 
IFCKASE.E0.3I Tst./«.786/T*.2S811 
Fle3.E16«EXP <-8.81/Tl«0UNl 
F2*3«E16*EXP (oA.eS/TI^OUNl 
F3«3.E16«EXP 1-3.02/Tl^tlUNl 
F6«F3 

01«2.68E13*EXP C-,66/TI *0001 
82«i.3E16*EXP 1-2.69/Tl^OUNl 
B3«1.33E15*EXP I-10.95/T1«OUN1 
86*3. 12E1 ••EXP C-12.S1/TI«0UM1 
T«TS8»E 
TS8VE«T 

IF(KASE.E0.2I T>1.261*.e5926*T 

F6*9.€6E18*EXP (-62.2/Tl/T*0UNl 

F7«8.00E1E«EXP <>S2.S/TI/T*0UN1 

86«1.E17*0U02 

87*1.E16«0UN2 

T«TSAYE 

FS*1.8SE17*EXP f-S6./TI/T*0UNl 
F8«S.80E1«*EXP (•60.6/Tl/T*DUNl 
8S*1.E^.6*0UN2 
08«6.E16«0UN2 
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0UNlstVC2t»Vf6MVt3ll«.S 
0UH23(V(1I*VI6I !••$»¥ C3) 

0UNSsY(l»«*S«Vf6l«V(3l 

OUH%sr l*v (1 1 «DUI*1»B1* Y 12 1 «Y ( 6 1 

0UNSsF2*Yt2l«0UN2*a2«VU»«Y(6l 

0UN6sF3«Y(6l*0UN2*e3*V(ll«Y(3l 

0UN7sFli«YfSI«YI6l«a4«Y(21*YC3l 

ounss (F2«.s>a7«suHvt • Y(2i ♦aa^Y («i 

OUK9S IFt«.S«a7*SUaYt*Ylll *ai*Y (61 

OOaiO«(F2*«5-ail*Y(2M(82-Fl».5l*Y(ll 

OUHll**Fl»,5-a3l*Y(l»-F3*Y(6l 

0UM12=F1»0UH1-83»Y(3I -F3*.S*Y(6I 

OOHl 3a ( F8 •SUHY^Fl *Y ( I M *. 5 

0UNl«iaB6«Y(ll*S(mY>F3*0UH3 

OOat5«2.»F6*Y(6» 

0UN16aSUHY«V(tl 

OUNl 7ae6«SUNY«Y (61 

ai2a 0UN9'-’F2*0Ua2 

B2la0U»*8-FX*0Wl 

019a (F6-F5I asuNYoFZaY (21 ♦OUNli 

029a (F2-86I «Y (2l>0Uai3 

e91a0UM12*a21*CUN8»0Uni7 

827aSuHY« (F7>F8/2. I »OUN10*OUH1S 

879».F6*SUMY-0UH11* (2. •04-F2 1 • Y ( 2 1 

erra-(CU814»S0aY*Fr*(Fl*.5*B2»*Y(ll*(ei*F2*.5l*Y(2M2.»0UM15l 

892a*86«Y(3t 

822a-SU0Y«(2«*B8*Y(2l »8r«Y(ll I •81«Y(6I»F2«0UN2-0UN13»B92 
Blla0Uai2-F5*SuaY*.5-(F2*.5»e7*SU«YI*Y(2l-82*Y(6l-0Uai7-2.*85*0Wl 
16 

a9YaOUN16»0UN19 

899aOUNll*(Fl«. 5*Y(1I ♦F6*SUHY»8«.*Y(2I » 
aria- (ouai2*ouH8»ou0i7i 
8T 2a 2 . *8<i • Y ( 3 1 -0UH9-F2*0O«2 
8l7aSUHY*(F7oF$/2.l>OUNlO«OUai<»*2.«F3*OUN3 
CCla0UNS*0UN6«0(fft6»(e6*Y(6Me5*Y(l)>87aV(2ll«0UH16 
CC2s 0UM4>DUN5«DUF7F (8 F^Y ( 1 » ♦Bi* Y ( 2 » » •SUHY* Y (2 1 
CC7aOUN6«OUNS*OU06»2.*OUN7« (06*V(6I>87*Y(2I )*0UH16 
CC9«0tM6-0UN7-B6*Y(6l*0UNl6 
16 Bea6«n0A*(Y(5l«(Y(2MY(6l»Y(3ll«.5l 
CCa6AHNA«YY(6l»Y(3IF(V(lt*V(6tl«.SI 
AO ( 1 t 1 I «eii »oci T AoFi*ea 
A0(l«2laei2»F2«CC 
A0(l*3>«ei7»F3*CC 
A0(l«6)«ei9 
A0(2*ll«e21FFl«8B 
AO (2 • 2) «e22»0eLT A«F2*CC 
A0(2.3la627 
A0(2*61ae29 
A0(3«i»«e7iFFi*ea 
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AOI3«2lBe72*P2*CC 

«0l3.3l«e77*0etTA-F3*CC 

A0(3*<i»*e79 

A0l4«&l«f>91 

A0f4«2l>-892 

*0l%.3lse97»F3*CC 

A0l%*«l«e99« DELTA 

cmi»cci*F5*sot»y*cc 

CI«2»*CC2*F8*SOFT*ea 

CK3I«CC7 

CKVlsCCS 

SCALEsfl.C 

00 SO Isl,(» 

00 SO 

50 SCALE»AFAXlfSCALE*AOS(AO(I*Jtll 
00 SI I«1*L 

00 S2 

S2 AOII*JI«AO(I*Jt/SCALE 

51 CI(II*CI(1I/SCALE 

CALL HCRFAN(VN.OT«AOtV«CI*ee«CC«SCALEI 
90 00 90 

IFCTNfJI.GE.O.O) CO TO 90 
OTsOT/tO. 

KTEST«KT€ST*1 
IF(KTCST-3» 92«27«27 
90 CONTINUE 
Dun«o.o 
00 t J«1.C 

1 0UN*0UM»VN(JI«2UI 
RH0Na0UH/«l.-ALFHAt7) ) 
YN(7I«RM0N*ALPHA<7>/Z(7t 
SUHTN>0.« 

00 2 J«L*7 

2 SUHYN«SUNYN*rN(JI 
TT«PRE$S/SUHYN 
DO A J*l«6 

h ALPNAIJ|sYNtJt«2IJI/RH0N 
AM«0.0 
8H«0«C 
CH«0«0 

00 SOS 1«1*7 

IFIft-TUin S02«S01*S01 
sot eH«qM-Em*ALPHA<It*.S 

CH«CH »0 < I I • ALPH A ( I I 

GO to SOS 

S02 IF(rT*TO(Itl S03*S03tSQA 
SOS eH«aN«e(n*ALPHA(ii«.s 
CH«CM*G<It«ALPHA(II 
GO TO SOS 


ORIGINAi; PA(aTB 

OP POOR QUALmj 
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sell 

BM»8H«AlPHA(II*(Ctll«T0III>e(I)*.SI 
CH*CH*«LPHA(II*ie<II»CfII*T0m«*2t 
505 CONTI MIC 
CHaCHHI 

IFIAHI 5er.5ec.507 
50e T«CH/BN/2. 

60 TO 50e 

507 T«CeH»SQRT tBN*BM-AH«CHII/AH 
50B CONTINUE 

IPIJJI81. 31.22 
31 ERR1>TT*T 

IF(ABSITT/T-1.0I.LC.EPSI GO TO 27 
GANlsCANNA 
GANNAs.se •GANNA 
ISe 6AN23GANNA 

DELTAS (CANNAol • I 70T 

IF a4«a4l•l^ 

a% If IIC0UNT.EQ»1I GO TO IH 
TsTSAVE 
GO TO S 
Em2«TT*T 

IFUeSCTT/T-i«OiaC«£PSI GO TO 27 
as GANHA«GaMl-fRRl*CGAH2«GAMll/fCRR2-£mil 
GAN1«GAII2 
ERRi«ERRa 
GO TO 130 

12 NRITEIGtiai 

13 F0RNAmN0t23N JJ IS GREATER THIN JJJI 
27 TN*T 

00 2 a B|9lt7 

2a ALPHNIJIsALRHAUI 
OtsOT^tlHEO 
PRESSI«PSSSS 
RETURN 
ENO 



